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ABSTRACT 


When  first  introduced  to  orbital  mechanics,  students  often  experience 
difficulty  in  visualizing  a  satellite's  actual  path  through  space.  The  Satellite 
Maneuver  Evaluation  Tool  (SMET)  seeks  to  alleviate  that  problem.  SMET  is  a 
three-dimensional  color  graphics  simulation  of  satellites  in  flight.  It  allows  a 
student  to  interactively  modify  a  satellite's  orbital  parameters  and  see  the  effects 
as  the  satellites'  positions  are  updated  continuously  with  respect  to  time.  The 
user  can  change  parameters  by  defining  a  maneuver  or  by  directly  entering  a 
change  through  keyboard  or  dial  inputs.  SMET  offers  the  user  the  capability  to 
demonstrate  difficult  concepts,  as  well  as  a  method  to  simulate  actual  satellite 
maneuvers.  Instructors  can  videotape  SMET  sessions  for  classroom 
demonstrations.  Instantaneous  images  can  also  be  saved  for  redisplay  or  for 
printout. 

This  thesis  provides  a  background  on  the  mathematical  formulas  modeled  by 
SMET.  It  also  includes  a  detailed  user’s  guide  for  SMET. 
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I.  INTRODUCTION 


Students  in  the  Space  Systems  Operations  and  Space  Systems  Engineering 
curricula  at  the  Naval  Postgraduate  School  are  required  to  develop  a  thorough 
understanding  of  orbital  mechanics.  They  use  equations,  tables  of  data  and  two- 
dimensional  static  illustrations  to  see  the  relationship  between  a  satellite's  orbital 
parameters  and  its  path  in  space.  These  students,  most  of  whom  arrive  with  little 
or  no  background  in  the  field,  often  have  difficulty  in  visualizing  the  actual  path 
of  a  satellite's  flight.  There  is  no  readily  available  tool  assist  in  learning  orbital 
mechanics;  a  tool  which  would  allow  a  student  to  interactively  modify  a  satellite’s 
orbital  parameters  and  see  the  effects. 

A  mathematical  relationship  can  be  represented  as  two-dimensional  graph.  A 
graph  conveys  information  more  readily  to  human  beings  than  does  a  table  of  the 
underlying  numbers.  But  if  more  than  two  variables  can  be  manipulated  such 
representations  quickly  loose  there  utility.  Adding  a  third  dimension  to  a  graph 
and  varying  it  with  time  can  be  used  to  portray  more  information.  But  that 
format  may  not  be  intuitive  to  the  user.  Visual  simulation  is  another  approach  to 
be  considered.  Visual  simulation  is  the  creation,  by  computer,  of  a  realistic, 
easily  modified,  moving  image  from  the  mathematical  model  of  a  phenomenon. 
Realism  implies  high-resolution,  color  graphics.  [Ref.  l:p.  3-5] 

The  Satellite  Maneuver  Evaluation  Tool  (SMET)  is  a  three-dimensional  color 
graphics  simulation  of  satellites  in  flight.  It  offers  the  user  three  ways  of 
observing  a  satellite's  position  with  respect  to  the  Earth.  In  one  window,  the 
user  is  placed  in  space  at  an  arbitrary  point  of  his  choice.  From  this  position,  he 
sees  the  motion  of  the  satellite  around  a  rotating  model  of  the  Earth.  Another 
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window  allows  the  user  to  see  that  portion  of  the  Earth  which  is  visible  to  the 
satellite.  A  third  window  shows  a  Mercator  projection  map  of  the  Earth  with  a 
plot  of  the  satellite's  ground  track.  After  a  user  enters  a  satellite’s  parameters, 
the  satellite's  position  is  continuously  updated  with  respect  to  time  in  all  three 
windows. 

Although  SMET  allows  the  user  to  change  parameters  through  the  use  of 
dials  or  keyboard  input,  its  greatest  utility  is  in  its  ability  to  let  the  user  execute  a 
satellite  maneuver  at  any  given  moment.  A  maneuver  is  entered  as  a  velocity 
change  (magnitude  and  direction)  or  as  a  combination  of  three  maneuver 
components  (flight  path  angle  change,  heading  angle  change  and  velocity 
magnitude  change).  The  user  can  also  enter  the  burnout  parameters  of  a  launch 
vehicle  in  order  to  place  a  satellite  into  its  orbit. 

The  user  is  also  able  to  read  a  satellite’s  position  and  velocity  continuously  in 
different  coordinate  systems  and  unit  types.  By  allowing  the  user  to  enter  as 
many  satellites  as  he  desires,  SMET  also  lets  him  visualize  the  relationship 
between  several  satellites.  As  an  example,  the  user  could  enter  several  identical 
satellites  but  vary  one  of  the  parameters  to  visually  compare  the  effect  of  that 
variation.  SMET  also  allows  students  to  create  illustrations  and  to  program 
simulations. 

SMET  is  written  for  any  Silicon  Graphics  IRIS  4D/GT,  GTX  or  VGX 
workstation.  A  personal  computer  implementation  was  considered  and  rejected 
for  two  reasons.  One,  there  is  no  commonly  available  set  of  three-dimensional 
color  graphics  routines  for  personal  computers  that  offers  the  functionality  to 
generate  the  images  desired  for  SMET.  Second,  graphics  workstations  are 
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becoming  more  affordable  all  the  time.  An  IRIS  is  not  an  extravagant  piece  of 
hardware.  Currently,  there  are  at  least  10  IRIS  workstations  around  the  campus. 

The  next  three  chapters  provide  a  background  on  the  mathematical  equations 
used  in  SMET’s  simulation  of  satellite  motion,  maneuver  and  launch.  The 
equations  are  not  derived,  in  most  cases,  because  that  was  considered  outside  of 
the  scope  of  this  document.  Those  derivations  can  be  found  in  the  refe^nces  or 
other  publications  listed  in  the  bibliography.  All  angles  in  trigonometric 
functions  are  in  radians.  Inverse  trigonometric  functions  may  require  additional 
quadrant  checking  that  is  not  shown.  Although  distances  and  velocities  can  be 
displayed  in  metric  or  english  units,  SMET  actually  stores  values  and  performs 
all  computations  in  canonical  units.  This,  however,  is  transparent  to  the  user. 
Chapter  V  provides  detailed  instructions  on  how  to  use  SMET.  Chapter  VI 
discusses  ways  to  modify  the  source  code  for  specific  simulations  not  currently 
modeled  in  SMET.  Appendix  A  is  a  consolidated  list  of  common  variables  in 
orbital  mechanics.  Appendix  B  has  procedures  for  recording  SMET  sessions  on 
video  tape.  It  also  provides  directions  on  how  to  print  out  Postscript  files 
generated  by  SMET.  Appendix  C  is  the  source  code  listing  for  SMET. 


II.  SATELLITE  MOTION 


A.  CLASSICAL  ORBITAL  PARAMETERS 

A  satellite's  path  in  space  can  be  described  by  five  independent  quantities. 
These  five  quantities,  called  "orbital  elements"  or  parameters,  are  sufficient  to 
completely  describe  the  size,  shape  and  orientation  of  the  orbital  ellipse.  A  sixth 
element  is  required  to  describe  the  satellite's  position  along  the  path  at  a 
particular  point  in  time.  [Ref.  2:p.  58] 

There  are  several  sets  of  variables  used  to  describe  a  satellites  position. 
SMET  uses  the  classical  set  of  six  orbital  elements.  The  semi-major  axis  (a)  and 
eccentricity  (e)  describe  the  size  and  shape  of  the  ellipse.  The  inclination  (i), 
longitude  of  the  ascending  node  (h)  and  argument  of  perigee  (g)  describe  the 
orientation  of  the  ellipse  are  depicted  in  Figure  2-1.  The  epoch  (tc)  is  the  time  of 
perigee  passage.  It  is  used  as  the  reference  from  which  to  compute  the  satellite's 
position  with  respect  to  time. 


B.  OTHER  STATIC  PARAMETERS 


Based  on  the  orbital  elements,  several  other  parameters  are  computed  to  help 
describe  the  orbit.  These  elements  do  not  change  with  respect  to  time.  The 
following  equations  show  how  to  compute  the  semi-minor  axis  (b),  period  (T), 
perigee  range  (rp)  and  apogee  range(ra). 


b 


t=2,c  VI" 

rp  =  a  (1  -  e) 
ra  =  a  (1  +  e) 


[Ref.  2:p.  33] 

[Ref.  2:p.  33] 

[Ref.  2:p.  25] 
[Ref.  2:p.  25] 
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Figure  2-1.  Orbital  Elements 
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A  transformation  matrix  R  is  also  computed  using  the  longitude  of  the 
ascending  node  (h),  argument  of  perigee  (g)  and  inclination  (i).  This  matrix  is 
used  in  drawing  the  satellite  and  its  orbital  ellipse.  It  is  also  used  to  convert 
perifocal  coordinates  into  geocentric-equatorial  coordinates. 

cos(h)cos(g)  -  sin(h)sin(g)cos(i)  -cos(h)sin(g)  -  sin(h)cos(g)cos(i)  sin(h)sin(i) 
sin(h)cos(g)  +  cos(h)sin(g)cos(i)  -sin(h)sin(g)  +  cos(h)cos(g)cos(i)  -cos(h)sin(i) 

sin(g)sin(i)  cos(g)sin(i)  cos(i) 

[Ref.  2:p.  83] 


C.  KEPLER  PROBLEM 

Given  a  satellite's  position  with  respect  to  its  classical  elements,  it  is  fairly 
simple  to  compute  the  elapsed  time  since  epoch.  This  is  described  at  the  end  of 
Chapter  III.  It  is,  however,  more  difficult  to  compute  a  position  given  the 
elapsed  time.  First  the  mean  anomaly  (M)  must  be  computed  using  the  following 
equation. 


M 


(t  - 10)  -  2k7t  +  Mo 


[Ref.  2:p.  220] 


SMET  defines  epoch  (t0)  at  perigee.  Since  the  mean  anomaly  and  true 
anomaly  are  both  zero  at  perigee,  the  last  term  (M0)  will  drop  out  from  the 
previous  equation.  The  next  step,  computing  the  eccentric  anomaly  (E),  is  more 
difficult.  To  do  this,  the  following  equation  must  be  solved  for  E.  It  is,  in 
classical  terms,  referred  to  as  the  Kepler  problem  [Ref.  2:p.  220]. 


M  =  E  -  e  sin  E 


[Ref.  2:p.  220] 
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This  equation  cannot  be  inverted  in  closed  form  to  provide  a  function  for  E 
in  terms  of  M.  Therefore,  another  approach  is  required.  The  following  Newton 
iteration  scheme  is  used  to  solve  for  E. 

E0  =  M 

Mn  =  Eji  -  e  sin  En  [Ref.  2:p.  221] 

M  -  M 

E-*'  =  E°+  l-e  cos  En  [Ref.  2:p.  2221 

The  second  and  third  equations  are  repeated  until  M  -  Mn  becomes 
sufficiently  small.  The  true  anomaly  (0)  can  be  then  computed  from  E. 

9-cos-'(e?cosCE:El)  [Ref.  2:p.  187] 

The  true  anomaly  is  depicted  in  Figure  2-1. 

D.  OTHER  CHANGING  PARAMETERS 

Once  the  position  of  the  satellite  (i.e.,  true  or  eccentric  anomaly)  is 
computed,  other  descriptions  of  its  position  can  be  computed.  First  the  range  (r) 
of  the  satellite  is  computed  using  the  eccentric  anomaly: 

r  =  a  (1  -  e  cos  E)  [Ref.  3:p.  C-7] 

Next,  the  coordinates  in  the  perifocal  coordinate  system  are  computed: 

P  =  r  cos  6,  [Ref.  2:p.  72] 

Q  =  r  sin0,  [Ref.2:p.72] 

W  =  0  by  definition. 
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The  coordinates  in  the  geocentric-equatorial  coordinate  system  (X,  Y  and  Z) 
are  then  computed  with  the  transformation  matrix  R. 


x* 

-p- 

Y 

-Z. 

=  R 

Q 

LwJ 

[Ref.  2:p.  83] 


The  last  set  of  coordinates  to  be  computed  are  the  latitude  and  longitude  or 
the  sub-satellite  point.  Using  simple  geometry,  we  get 


lat  =  sin-1 


Ion  =  tan-1 

These  coordinates  are  fixed  in  space.  Note  that  the  earth  rotates  at  an 
angular  velocity  of  (0®.  In  order  to  compute  the  longitude  with  respect  to  the 
rotated  earth,  subtract  the  earth's  current  rotation  and  allow  for  an  arbitrary 
number  of  rotations  to  correct  the  previous  equation.  The  corrected  equation 
for  longitude  is 


Ion  =  tan1 


(0®  t  +  2k K. 


Although  the  magnitude  of  the  velocity  (v)  is  not  a  description  of  the 
position,  it  does  provide  some  insight  as  to  the  motion  of  the  satellite.  It  is  also 
instrumental  in  computing  orbital  maneuvers.  The  following  equation  provides 
the  instantaneous  velocity  of  the  satellite. 


v  = 


[Ref.  3:p.  C-ll] 
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E.  PERTURBATIONS 


A  perturbation  is  a  variation  in  the  path  of  the  satellite  not  explained  by  the 
Keplerian  model  for  two-body  motion  used  up  to  this  point.  Perturbations  result 
from  the  gravitational  influence  of  other  bodies,  atmospheric  drag  and  the  oblate 
shape  of  the  earth.  The  two  principle  effects  are  the  regression  of  the  line  of 
nodes  and  the  rotation  of  the  line  of  apsides  [Ref  2:p.  156].  These  two  effects 
change  the  longitude  of  the  ascending  node  and  argument  of  perigee, 
respectively.  When  turned  on,  SMET  uses  the  following  equations  to  model 
these  two  principle  perturbations. 


h  = 


-3  VrJjl.r|  cosi 

2  a3-5  (1  -  e2) 


[Ref.  3:p.  4-4] 


•  _  Vii  J2  Re 

8"  a3-5  (1  -  e2)2 


[Ref.  3:p.  4-5] 


J2  is  the  second  order  coefficient  for  perturbation  accelerations  due  to  the 
earth’s  nonsphericity.  Other  perturbations  are  more  computationally  in  ensive 
and  provide  changes  to  small  to  be  visible  with  SMET.  Therefore,  SMET  makes 
no  attempt  to  model  them. 


F.  FIELD  OF  VIEW  CALCULATION 

The  Satellite  View  of  SMET  demonstrates  how  the  altitude  of  a  satellite 
effects  the  amount  of  the  earth's  surface  visible  to  the  satellite.  This  amount  is 
referred  to  as  the  field  of  view  (FOV)  of  the  satellite.  Note  that  the  angle  a  in 
Figure  2-2  represents  half  of  the  FOV.  If  the  range  (r)  is  expressed  in  terms  of 
earth  radii  (i.e.,  canonical  units),  the  problem  is  simplified  to: 

1 

sin  a  = 

r 
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Thus 
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III.  SATELLITE  LAUNCH 


A.  GENERAL  DESCRIPTION 

When  a  satellite  is  launched,  it  is  normally  attached  to  the  top  of  a  rocket. 
While  the  rocket  is  burning,  the  mass  and  velocity  are  constantly  changing. 
Additionally,  the  effects  of  gravity  and  atmospheric  drag  change  continuously. 
At  the  end  of  the  rocket  bum,  the  satellite  will  separate  from  the  final  stage  of  its 
rocket.  The  end  effect  of  the  rocket  burn  is  to  place  the  satellite  at  a  particular 
point  in  space  with  a  particular  velocity.  At  this  burnout  point,  the  velocity 
vector  of  the  satellite  corresponds  to  a  unique  orbital  path.  SMET  allows  the 
user  to  enter  parameters  associatied  this  burnout  point  to  start  a  satellite's 
motion.  Modeling  the  complex  interaction  of  the  large  number  of  variables  that 
define  the  actual  bum  sequence  is  beyond  the  scope  of  SMET. 

B.  PARAMETERS 

The  values  of  the  following  variables  at  the  time  of  burnout  are  used  by 
SMET  to  compute  the  orbital  parameters  of  the  satellite.  This  is  by  no  means  the 
only  way  to  define  the  burnout  of  the  rocket.  They  were  chosen  because  they  are 
easy  to  visualize. 

v0  Velocity  (with  respect  to  a  rotating  earth) 

<{>  Latitude  in  earth-fixed  coordinates 
A  Longitude  in  earth-fixed  coordinates 

Yo  Flight  Path  Angle  (measured  upward  from  the  local  horizon  plane) 
\\f0  Launch  Azimuth  (measured  clockwise  from  North) 
alt  Altitude  above  the  surface  of  the  earth 
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C.  EARTH  SPIN  CORRECTION 

The  earth  rotates  about  its  axis  in  23  hours  and  56  minutes.  This  produces  a 
surface  velocity  at  the  equator  of  approximately  465.10  m/sec  (1524  ft/sec) 
[Ref.  2:p.  306].  The  surface  velocity  is  a  function  of  latitude,  decreasing  as  the 
poles  are  approached.  This  explains  why  eastward  launches  or  equatorial  launch 
points  are  more  fuel  efficient.  SMET  adds  this  induced  velocity  to  the  satellite 
velocity  defined  by  the  burnout  parameters. 

The  velocity  is  broken  up  into  components  in  a  satellite  fixed  coordinate 
system.  This  is  depicted  in  Figure  3-1.  This  system  has  a  radial  component,  a 
North/South  component  and  an  East/West  component.  Using  simple  geometry, 
the  following  equations  provide  the  solution  for  the  velocity  components. 

vR  =  v0  sin  y0 
vn/s  =  v0  cos  Yo  cos  \\fQ 
ve/w0  =  v0  cos  Yo  sin  y0 
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The  earth  rotation  component  (vSpin)  is  then  computed  from  the  latitude.  R® 
is  the  radius  of  the  earth  and  0)©  is  the  angular  velocity  of  the  earth. 

vspin  =  R©co©  cos  <J)  [Ref.  2:p.  307] 

SMET  adds  vSpiN  to  ve/w0  and  computes  a  new  horizontal  velocity 

component.  The  following  equations,  based  on  Figure  3-2,  show  the  corrected 
values  of  velocity  magnitude,  flight  path  angle  and  azimuth. 


Ve/W  =  VE/W0  +  VSPIN 


D.  ORBIT  PARAMETER  CALCULATIONS 

SMET  uses  the  corrected  burnout  parameters  to  compute  the  orbital 
parameters.  The  inclination  is  a  function  of  the  azimuth  and  latitude. 


i  =  cos-!  (sin  y  cos  <}>) 

The  range  is  a  function  of  the  altitude. 

r"  1  +Ri 

The  semi-major  axis  is  a  function  of  the  range  and  velocity, 
r 


[Ref.  2:p.  142] 


a  = 


o  . 


r  v 


[Ref.  3:p.  C-5] 


The  semi-latus  rectum  is  a  function  of  the  range  and  horizontal  velocity. 

P  =  [Ref.  3:p.  C-7] 

The  eccentricity  is  a  function  of  the  semi-major  axis  and  semi-latus  rectum. 


■V 


[Ref.  3:p.  C-6] 


The  true  anomaly  is  a  function  of  the  range  and  semi-latus  rectum. 


[Ref.  3:p.  C-9] 


The  arguments  of  latitude  and  longitude  are  functions  of  the  latitude  and 
azimuth. 


sin  0 


u  =  tan-1  , 

(cos  <p  cos 

,  ['sin  <b  sin  ty 

©  =  tan-1  - * - *- 

I  cos  y 


[Ref.  3:p.  A-7] 
[Ref.  3:p.  A-7] 
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The  argument  of  perigee  is  a  function  of  the  argument  of  latitude  and  true 
anomaly. 


g  =  u-0  [Ref.  3:p.  C-9] 

To  compute  the  longitude  of  the  ascending  node,  note  that  X  -  (0  represents 
the  distance  between  satellite's  position  in  a  rotational  coordinate  system  and  its 
position  fixed  coordinate  system.  Also  note  that  0)©t  represents  the  distance 
between  an  ascending  node  of  0°  (fixed  in  space)  and  an  ascending  node  that  has 
rotated  as  much  as  the  earth. 


h  =  X  -  (0  +  (0©  t 

As  discussed  in  Chapter  II,  SMET  computes  the  true  anomaly  as  a  function 
of  the  elapsed  time  since  epoch  (to).  Therefore,  after  the  true  anomaly  has  been 
computed,  the  eccentric  anomaly  (E)  and  mean  anomaly  (M)  are  computed  using 
the  eccentricity  (e)  and  true  anomaly. 


E  =  cos'1 


f  e  +  cos  6  ^ 
l  +  e  cos  9 J 


M  =  E  -  e  sin  E 


The  last  parameter,  the  epoch,  is  then  computed  using  the  mean  anomoly,  the 
semi-major  axis  (a)  and  current  time  (t). 
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IV.  SATELLITE  MANEUVERS 


A.  GENERAL  DESCRIPTION 

A  satellite  maneuver  is  executed  by  changing  the  velocity  of  the  satellite. 
This  change  to  the  velocity  vector  can  be  in  its  magnitude,  direction  or  both.  A 
maneuver  results  in  a  new  set  of  orbital  parameters.  A  velocity  change  (Av)  is 
modeled  as  an  instantaneous  impulse.  Therefore,  the  satellite’s  position  or  range 
does  not  change  while  the  velocity  does.  The  Av  vector  can  be  broken  up  into 
component  vectors  (Avy,  Avv  and  Ava)  which  cause  different  effects 
simultaneously.  These  three  vectors  cause  a  flight  path  angle  change,  a  heading 
angle  change  and  a  change  in  the  magnitude  of  the  velocity  vector,  respectively. 
Both  A v y  and  Avv  change  the  direction  of  the  velocity  vector  but  do  not  change 
its  magnitude.  SMET  computes  the  effects  of  each  of  the  maneuvers 
independently  and  then  updates  the  time  of  epoch.  The  time  of  epoch  is 
computed  such  that  the  elapsed  time  will  place  the  satellite  in  the  same  position 
based  on  the  newly  computed  parameters. 

B.  FLIGHT  PATH  ANGLE  CHANGE 

A  Avy  will  cause  a  change  in  the  flight  path  angle  (y).  A  flight  path  angle 
change  is  a  maneuver  in  the  orbital  plane.  The  flight  path  angle  is  measured 
upward  from  the  local  horizon  plane  to  the  velocity  vector  [Ref.  3:p.  2-22].  The 
local  horizon  plane  passes  through  the  satellite  and  is  perpendicular  to  the  line 
passing  through  the  center  of  the  earth  and  the  satellite.  Avy  is  chosen  such  that 
the  magnitude  of  the  velocity  vector  is  unchanged  (i.e.,  Vi  =  V2).  These  vectors 
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form  an  isosceles  triangle  as  can  be  seen  in  Figure  4-1.  Given  Ay,  values  for  the 
eccentricity  (e),  argument  of  perigee  (g),  and  true  anomaly  (0)  can  be  computed. 


Figure  4-1.  Flight  Path  Angle  Change  Vectors 


The  eccentricity  changes  as  the  impulse  changes  the  shape  of  the  orbit.  The 
argument  of  perigee  and  true  anomaly  change  as  a  result  of  the  rotation  of  the 
line  of  apsides.  This  can  be  seen  in  Figure  4-2.  Since  AvY  has  no  out-of-plane 
component,  there  is  also  no  change  in  the  inclination  (i)  or  the  longitude  of  the 
ascending  node  (h).  The  following  equation  shows  that  the  semi-major  axis  (a) 
does  not  change  since  the  range  and  the  magnitude  of  the  velocity  do  not  change. 

a  = - X—2  [Ref.  3:p.  C-5] 

2  -  — 

M- 

The  eccentricity  can  be  computed  given  Ay  by  combining  the  following  two 
equations  for  the  cosine  of  the  flight  path: 


Vu p 

cosy  =— ■ 

1  rv 


[Ref.  3:p.  3-15] 


cos  y2  =  cos  yi  (cos  Ay  -  tan  yi  sin  Ay). 
This  results  in 


VupT  _  Vup7 

r2v2  HV! 


(cos  Ay  -  tan  yi  sin  Ay). 


[Ref.  3:p.  3-15] 


Recalling  that  the  range  and  magnitude  of  the  velocity  are  unchanged,  they 
can  be  removed  from  the  previous  equation. 

=  (cos  Ay  -  tan  7i  sin  Ay) 

The  semi-latus  rectum  is  computed  with 

p  =  a  (1  -  e2).  [Ref.  2:p.  24] 

Substituting  this  into  the  previous  equation  while  recalling  that  the  semi¬ 
major  axis  remains  unchanged,  results  in 

"\j\  '  p  =  (cos  Ay  -  tan  y}  sin  Ay), 

1  -  =  (1  -  e2)(cos  Ay  -  tan  yi  sin  Ay)2. 

The  tangent  of  the  flight  path  angle  can  be  computed  from  the  eccentricity 
and  eccentric  anomaly. 

e,  sin  E, 

tan  =-Jp=dr  [Ref.  3:p.  C-10] 

Vl-e2 

Thus  a  new  eccentricity  can  be  computed  with  Ay,  and  Ej. 


e2  = 


V 


1  •  (1  -  e?) 


cos  Ay 


e2  sin  Ej  sin  AyYi2 


vr 


e2 

el 


JJ 


Another  equation  relating  the  flight  path  angle  to  the  satellite's  position 
provides  the  means  for  computing  the  new  true  anomaly. 


0  =  tan 


.1  ^E_tan_yj 


[Ref.  3:p.  C-9] 

Using  tan  ylt  which  was  computed  earlier,  we  can  compute  the  tangent  of  the 
new  flight  path  angle. 


tan  y,  +  tan  Ay 

tan  y2  =  - u - — - 

1  -  tan  yj  tan  Ay 


[Ref.  3:p.  3-19] 
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Using  p2  =  a  (1  -  e2),  the  new  anomaly  can  thus  be  computed. 


Since  the  satellite  has  not  moved  and  the  line  of  nodes  has  not  changed,  the 
argument  of  latitude  (u)  remains  unchanged.  Knowing  that  g  =  u  -  9,  we  can 
compute  the  new  argument  of  perigee  using  the  old  argument  of  perigee  and  old 
true  anomaly  with  the  new  true  anomaly  [Ref  3:p.  3-8]. 


g2  =  gi  +  9i-e2 

C.  HEADING  ANGLE  CHANGE 

A  Avv  will  cause  a  change  in  the  heading  angle  (y).  A  heading  angle  change 
is  a  maneuver  out  of  the  orbital  plane.  The  heading  angle  is  the  angle  in  the  local 
horizon  plane  measured  from  the  line  through  the  satellite  and  the  north  point 
and  the  velocity  vector  [Ref.  3:p.  2-22].  Again,  Avv  is  chosen  such  that  the 
magnitude  of  the  velocity  vector  is  unchanged  (i.e.,  Vi  =  v2).  These  vectors 
form  an  isosceles  triangle  as  can  be  seen  in  Figure  4-3.  With  Aij/,  values  for  the 
inclination  (i),  argument  of  perigee  (g)  and  longitude  of  the  ascending  node  (h) 
can  be  computed. 
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The  inclination  changes  as  the  orbital  plane  rotates  about  the  line  connecting 
the  center  of  the  earth  and  the  current  satellite  position.  This  rotation  can  be 
seen  in  Figure  4-4.  This  rotation  also  changes  the  point  where  the  orbital  plane 
intersects  with  the  equatorial  plane,  causing  the  ascending  node  to  change. 
Finally,  as  the  plane  rotates,  the  distance  between  the  ascending  node  and  the 
current  position  as  measured  along  the  orbital  path  (the  argument  of  latitude  (u)) 
changes,  causing  the  argument  of  perigee  to  change.  Since  A\|f  has  no  in-plane 
component,  there  is  no  change  to  the  shape  or  in-plane  orientation  of  the  orbit. 
Thus,  the  eccentricity  (e),  semi-major  axis  (a)  and  true  anomaly  (0)  remain 
unchanged. 

Given  A\|/,  the  new  inclination  and  argument  of  latitude  can  be  computed 
using  the  following  equations: 


cos  h  =  cos  ii  cos  A\y  +  sin  ii  cos  Ui  sin  A\j/,  [Ref.  3:p.  3-6] 
sin  \2  cos  U2  =  sin  ii  cos  Ux  cos  A\|/  -  cos  ij  sin  A\j f,  [Ref.  3:p.  3-6] 
ui  =  gi  +  0i.  [Ref.  3:p.  3-8] 

The  new  inclination  and  argument  of  latitude  are  found  by  rearranging  the 
previous  equations. 


i2  =  cos-1  (cos  ii  cos  A\[/  +  sin  ij  cos(gi+0i)  sin  A\y  ) 
j^sin  ii  cos(gi+9i)  cos  A\y  -  cos  ii  sin  Ayr 

sin  i2 


e 


■) 


U2  =  COS’1 

Since  the  true  anomaly  is  unchanged  (i.e.,  0j  =  82),  the  new  argument  of 
perigee  is  simply: 


g2  =  u2  -  81. 
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Figure  4-4.  Heading  Angle  Change 
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The  new  longitude  of  the  ascending  node  is  computed  by  noting  that  the 
change  in  the  longitude  of  the  ascending  node  is  equal  to  the  change  in  the 
argument  of  longitude  [Ref.  3:p.  3-11].  The  argument  of  longitude  is  computed 
as  follows: 


(0  =  tan-1 


/sin  u  cos  h 
^  cos  u  / 


[Ref.  3:p.  A-8] 


_  ,  rsin  ui  cos  in  ,  ^sin  U2  cos 

Thus  h2  =  hi  +  tan-1  — -7— - L  -  tan-1  — 

L  1  V.  cos  U1  )  ^  COS  U2  / 


D.  VELOCITY  MAGNITUDE  CHANGE 

A  Ava  will  change  the  magnitude  of  the  velocity  (V2  =  Vi  +  Ava)  without 
changing  the  direction  (Figure  4-5).  Based  on  the  new  velocity  magnitude, 
values  for  the  semi-major  axis  (a)  and  eccentricity  (e)  can  be  computed.  These 
changes  are  reflected  as  a  change  in  the  shape  of  the  orbit  (Figure  4-6).  Since 
the  heading  angle  and  flight  path  angle  are  unchanged,  the  parameters  defining 
the  orientation  of  the  orbital  plane  ( i.e.,  i  and  h)  are  unchanged  [Ref  2.:p.  3-20]. 
If  the  orbit  was  circular  prior  to  the  change,  the  argument  of  perigee  (g)  and 
true  anomaly  (0)  must  also  be  set,  since  those  two  parameters  have  little  meaning 
in  a  circular  orbit.  The  velocity  change  at  a  position  other  than  perigee  or 
apogee  will  also  cause  a  rotation  of  the  line  of  apsides.  Therefore,  the  argument 
of  perigee  and  true  anomaly  must  be  computed  as  they  were  for  a  flight  path 
angle  change. 


v 

1 


> 


AVa 


> 


Figure  4-5.  Velocity  Magnitude  Change  Vectors 
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Figure  4-6.  Velocity  Magnitude  Change 
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The  semi-major  axis  and  the  eccentricity  are  computed  with  the  following 
equations: 

[Ref.  3:p.  3-20] 


a2  = 


r  vl* 


e2 


2  v2 
VH2 


M-  a2 


where  vH2  =  v2  cos  y.  [Ref.  3:p.  3-20] 


The  flight  path  angle  (y)  does  not  change  since  the  velocity  vector  was  not 
rotated.  Therefore,  the  old  values  of  e  and  E  can  be  used  to  compute  y. 


e,  sin  E, 
tan  y  =-* 1 - L 


[Ref.  3:p.  C-10] 


V 1  -  ei 

If  the  original  orbit  was  circular  (i.e.,  ei  =  0),  the  fact  that  the  argument  of 
latitude  remains  unchanged  is  used  to  compute  a  value  for  the  argument  of 
perigee  and  true  anomaly.  The  current  position  will  be  perigee  if  there  was  an 
increase  in  velocity. 


g2  =  U  =  gi  +01 
02  =  U  -  g2  =  0° 

If  the  velocity  decreased,  the  satellite  will  be  at  apogee  of  the  new  orbit. 

g2=  u-  180°=  gi  +  0i-  180° 

02  =  u  -  g2  =  1 80° 

E.  UPDATING  EPOCH 

Changes  in  the  magnitude  of  the  velocity  and  flight  path  angle  cause  a  change 
to  the  true  anomaly  (0).  As  in  the  satellite  launch  equations,  the  ;poch  is 
computed  as  a  function  of  the  true  anomaly  and  the  current  time  with  respect  to 
the  eccentricity  and  semi-major  axis  of  the  new  orbit. 
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F.  IMPULSE  VECTOR  DIVISION 

Given  a  velocity  change  vector  (Av)  with  an  associated  yaw  and  pitch,  the  in¬ 
plane  and  out-of-plane  components  (Avp  and  Avop)  can  be  computed.  Using  Avp 
and  Avop,  the  desired  vector  components  (AvY,  Av¥  and  Ava)  can  be  computed. 
Unfortunately,  AvY,  Avv  and  Ava  are  not  orthogonal  vector  and,  therefore,  are 
not  easily  computed.  Figure  4-7  depicts  the  relationship  between  the  two  sets  of 
vectors.  The  same  set  of  vectors  is  depicted  again  as  part  of  Figure  4-10. 

Earlier  sections  noted  that  AvY  and  Avv  formed  isosceles  triangles  when 
added  to  the  velocity  vector.  These  isosceles  triangles  had  associated  angles  (Ay 
and  A\| /)  which  were  used  to  compute  the  flight  path  change  maneuver  or  heading 
change  maneuver.  Similarly  the  new  velocity  (V2),  and  not  Ava,  is  actually  used 
to  compute  the  new  orbital  parameters  for  a  velocity  magnitude  change 
maneuver.  SMET,  therefore,  does  not  use  Avy,  Avv  and  Ava.  SMET  actually 
computes  Ay,  A\j/andv2. 


The  yaw  associated  with  a  maneuver  is  used  to  compute  the  in-plane  and  out- 
of-plane  components  of  the  Av  vector. 

Avp  =  Av  cos(Y  aw) 

Avop  =  Av  sin(Yaw) 

The  pitch  associated  with  the  maneuver  is  used  to  compute  that  portion  of  the 
in-plane  component  (Avp)  which  is  in  the  direction  of  motion  and  that  portion 
which  is  perpendicular  to  the  direction  of  motion.  These  vectors  are  added  to 
the  original  velocity  as  depicted  in  Figure  4-8. 


Figure  4-5.  Computing  Flight  Path  Angle  Change 


This  right  triangle  provides  the  means  with  which  to  compute  the  change  to 
flight  path  angle  (Ay). 

Ay  =  tan-4  yffitch)  J 
^  vi  +  Avp  cos(Pitch)  J 

The  hypotenuse  of  the  right  triangle  is  computed  for  later  use. 

v\  +  Aw 'a  =  V(Avp  sin(Pitch))2  +  (vi  +  Avpcos(Pitch))2 
In  order  to  compute  the  heading  angle  change  (Ay),  the  out-of-plane 
components  of  Av  (Avop)  are  added  to  the  previously  computed  in-plane 
components  (v\  +  Av^).  Figure  4-9  shows  that  a  right  triangle  is  formed  again. 
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Figure  4-9.  Computing  Heading  Angle  Change 


This  right  triangle  provides  the  means  with  which  to  compute  the  change  to 
the  heading  angle  (Ay). 

Ay  =  tan-1  f~r^f . T'l 

Vvi  +  AvaJ 

The  hypotenuse  of  this  right  triangle  (v ’{  +  Ava)  is  V2.  This  can  be  seen  in 
Figure  4-10  which  combines  the  two  previous  figures. 

v2  =  v"  +  Ava  =  VAVop2  +  (v’i  +  Ava)2 
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V.  USER’S  GUIDE 


A.  OVERVIEW 

This  chapter  describe  SMETs  display  and  user  interface.  Figure  5-1  shows 
the  five  windows  used  by  SMET  -  the  Main  Window,  two  Auxiliary  Windows, 
the  Parameter  Window  and  the  Satellite  List  Window.  The  Main  Window  is  a 
large  square  window  at  the  center  of  the  screen.  There  are  two  smaller 
Auxiliary  Windows  below  the  Main  Window.  These  three  windows  are  called 
viewing  windows.  The  Parameter  Window  is  a  long  narrow  window  on  the  left 
side  of  the  screen.  The  Satellite  List  Window  is  a  long  narrow  window  on  the 
right  side  of  the  screen. 

SMET  updates  all  satellite  positions  continuously  with  respect  to  time.  One 
satellite  is  designated  as  the  Current  Satellite.  This  satellite  is  listed  at  the  top  of 
the  Satellite  List  Window  and  has  its  parameters  displayed  in  the  Parameter 
Window.  Additionally  the  satellite  viewpoint,  discussed  later,  belongs  to  the 
Current  Satellite.  All  user  inputs  affect  the  Current  Satellite.  Any  satellite  can 
become  the  Current  Satellite  with  the  use  of  the  Satellite  List  Window. 

B.  VIEWING  WINDOWS 

The  three  viewpoints  discussed  in  this  section  can  appear  in  any  of  the  three 
viewing  windows.  These  images  are  completely  interchangeable.  The  user  can 
switch  an  image  from  an  Auxiliary  Window  into  the  Main  Window  by  simply 
pressing  the  middle  mouse  button  while  the  cursor  is  in  the  desired  Auxiliary 
Window.  The  image  that  was  in  the  Main  Window  will  be  moved  to  the  selected 
Auxiliary  Window. 
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1.  OmniViewer  View 

The  OmniViewer  View  allows  the  user  to  observe  all  satellites  as  they 
orbit  around  a  rotating  earth  model  as  seen  in  Figure  5-2.  Axes  for  the 
geocentric-equatorial  coordinate  system  are  displayed  for  user  orientation.  An 
orbital  path  is  displayed  for  each  satellite.  The  user  can  choose  to  draw  it  as  a 
colored  line  or  as  a  colored  transparent  ellipse.  Normal  and  perigee  vectors  can 
also  be  displayed  for  any  satellite  to  assist  in  user  orientation. 

The  user's  viewing  position  is  referred  to  as  the  Viewer  X,  Y  and  Z  in 
the  Parameter  Window.  Changing  the  user’s  viewing  position  is  discussed  in 
SMET  System  Parameters.  It  can  be  set  automatically  or  adjusted  by  the  user. 
When  Automatic  OmniViewer  is  turned  on,  the  viewer  will  be  positioned  far 
enough  away  to  be  able  to  see  all  orbital  paths.  The  viewer  position  will  be 
adjusted  automatically  whenever  a  satellite  is  added  or  any  orbital  path  is 
changed  as  a  result  of  a  maneuver  or  user  input.  During  satellite  maneuvers  or 
launches,  Automatic  OmniViewer  can  be  turned  off.  This  allows  the  user  to  see 
the  event  on  a  single  scale. 

2.  Satellite  View 

The  Satellite  View  allows  the  user  to  see  that  portion  of  the  earth  visible 
to  the  Current  Satellite.  The  image  is  a  "Fish-eye  lense"  type  of  view.  The  user 
sees  a  circular  image  of  the  earth  encircled  by  a  ring  colored  in  the  Current 
Satellite's  color.  The  center  of  the  image  corresponds  with  the  Current 
Satellite’s  sub-satellite  point.  The  outer  edge  represents  the  horizon  as  seen  from 
the  altitude  of  the  Current  Satellite.  As  the  Current  Satellite's  altitude  changes, 
the  field  of  view  is  adjusted  to  maximize  the  size  of  the  image.  Therefore,  the 
size  of  the  image  does  not  generally  change.  Figure  5-3  shows  a  typical  view 
from  a  satellite  in  a  Molniya  orbit. 
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Figure  5-2.  OmniViewer  View 
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Figure  5-3.  Satellite  View 
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3.  Mercator  View 

The  Mercator  View  shows  a  map  of  the  world  as  seen  in  Figure  5-4. 
The  sun’s  relative  plot  is  depicted  as  a  yellow  transparent  ellipse  on  the  map. 
This  ellipse  does  not  attempt  to  represent  day/night  transitions.  Each  satellite 
will  have  its  sub-satellite  point  plotted  in  its  color  every  time  the  satellite’s 
position  is  updated.  Satellite  Parameters  discusses  how  to  turn  plotting  on. 

C.  PARAMETER  WINDOW 

The  Parameter  Window  displays  names  and  values  of  parameters  for  the 
Current  Satellite  and  the  SMET  system.  The  user  can  change  the  value  of  many 
parameters  with  mouse  and  keyboard  inputs.  Distance  and  velocity  parameters 
are  displayed  in  one  of  three  types  of  units.  The  three  types  are  canonical  (DU, 
DU/TU),  metric  (km,  km/sec)  and  english  (NM,  ft/sec). 

1 .  Satellite  Parameters 

There  are  five  sets  of  Satellite  Parameters  -  Classical,  Secondary, 
Fosition,  Options  and  Combined.  The  Classical  set  includes  the  parameters  that 
describe  the  shape  and  orientation  of  the  orbit.  Most  of  these  parameters  can  be 
changed  by  the  user.  The  Secondary  set  also  includes  parameters  that  describe 
the  shape  of  the  orbit.  Most  of  these  parameters  are  actually  computed  by  SMET 
for  display  and  cannot  be  changed  by  the  user.  The  Position  set  includes 
parameters  that  describe  the  satellites  position  in  various  coordinate  systems. 
These  parameters  are  also  computed  and  cannot  be  changed  by  the  user.  The 
Options  set  is  a  set  of  parameters  that  SMET  uses.  The  Combined  set  includes  all 
the  members  of  the  ether  sets  listed  above.  All  parameter  sets  display  the 
Current  Satellite's  name  at  the  top.  Figure  5-5  shows  the  Combined  set  of 
satellite  parameters  with  the  system  parameters  listed  at  the  bottom. 
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Figure  5-4.  Mercator  View 
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Satellite  Parameters; 


Name  Molnlva  Orbit 
Semi-Major  $656178  It 
Semi-Minor  17568.97  Is 
Eccentricity  0,75 
Inclination  63.43 
Ascending  Node  59.55 
Arg  <rf  Periapsis  270  00 
True  Anomaly 18097 
Epoch  25;54;00 
Period  11:58:02 
Apogee  46483 12  km 
Perigee  6640.45  km 
Range  46463  18  km 
Velocity'  1 47  km/sec 
X -18308-79  km 
Y  9851.25  km 
Z  41551.97  km 
P  -46456:53  km 
Q  786,01  km 
Latitude  63  42  N 
longitude  15009  W 
Color  Red 
Fill  Off 
Vectors  On 
Plot  On 


System  Parameters: 
Clock  On 

Current  Time  44:30:00 
Qock  Speed  00:30.00 


Political  Borders  On 
Perturbations  Off 
AutoOmniViewerOn 
Viewer  X  46483.12  km 
Y  637815  km 
Z  637815  km 
UnitsMetric 


Figure  5-5.  Parameter  Window 
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a.  Classical  Set 

The  Classical  set  of  satellite  parameters  includes  the  semi-major 
axis,  semi-minor  axis,  eccentricity,  inclination,  longitude  of  the  ascending  node, 
argument  of  periapsis,  true  anomaly,  and  time  of  epoch.  The  true  anomaly  is 
computed  by  SMET  as  a  function  of  the  elapsed  time  since  epoch.  It  cannot  be 
changed  by  the  user.  All  other  parameters  can  be  adjusted  by  the  user  as 
described  below. 

b.  Secondary  Set 

The  Secondary  set  of  satellite  parameters  includes  the  perigee 
distance,  apogee  distance,  time  of  epoch,  orbital  period,  current  range  and 
velocity.  Except  for  the  time  of  epoch,  all  of  these  parameters  are  computed  by 
SMET  and  cannot  be  changed  by  the  user. 

c.  Position  Set 

The  Position  set  of  satellite  parameters  displays  the  current  position 
of  the  satellite  in  three  different  sets  of  coordinate  systems.  The  X,  Y  and  Z 
coordinates  correspond  with  the  geocentric-equatorial  coordinate  system 
displayed  in  the  OmniViewer  window.  The  P  and  Q  coordinates  correspond 
with  the  perifocal  (P-Q-W)  coordinate  system.  Both  of  these  sets  of  coordinates 
are  displayed  in  the  current  unit  type  (ex.  metric).  The  Latitude  and  Longitude 
coordinates  correspond  with  the  earth-fixed  coordinate  system  seen  in  the 
Mercator  View  window.  None  of  these  parameters  can  be  change  by  the  user. 

d.  Options  Set 

The  Options  set  include  the  color  of  the  satellite  and  its  ellipse,  if 
the  orbital  ellipse  should  be  filled  or  empty,  if  the  normal  and  perigee  vectors 
should  be  displayed  and  if  the  satellite  should  be  plotted  on  the  Mercator  map. 
All  of  these  parameters  can  be  changed  by  the  user. 
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e.  Combined  Set 

The  Combined  set  includes  all  the  members  of  the  other  sets  listed 

above. 

2.  SMET  System  Parameters 

SMET  system  parameters  are  listed  below  the  set  of  satellite 
parameters.  This  set  of  parameters  includes  the  clock,  current  time,  clock  speed, 
detailed  maps,  political  borders,  perturbations,  Automatic  OmniViewer, 
OmniViewer  X,  Y,  and  Z,  and  units.  The  clock  parameter  is  used  to  start  and 
stop  the  clock.  Current  time  displays  the  elapsed  time  in  hours:minutes:seconds 
format.  It  uses  00:00:00  GMT  on  the  day  of  Vernal  Equinox  as  the  reference 
point.  Clock  speed  allows  the  user  to  set  the  number  of  minutes  between  position 
updates.  It  is  initially  set  to  30  minutes.  The  detailed  maps  parameter  allows  the 
user  to  switch  between  detailed  and  coarse  globe  and  map  data.  The  political 
borders  parameter  allows  the  user  to  turn  off  the  display  of  political  borders  on 
the  globe  and  map.  Both  of  these  options  allow  a  user  with  a  large  number  of 
satellites  to  speed  up  SMET  by  drawing  less  detailed  images.  They  are  both 
initially  set  to  on.  The  perturbations  parameter  turns  on  the  regression  of  the 
line  of  nodes  and  the  rotation  of  the  line  of  apsides.  This  is  initially  set  to  off. 
Automatic  OmniViewer  positioning  can  be  turned  on  and  off.  It  is  initially  set  to 
on.  If  any  of  the  OmniViewer  coordinates  (X,  Y  or  Z)  are  changed,  Automatic 
OmniViewer  positioning  is  turned  off.  The  OmniViewer  coordinates  are 
geocentric-equatorial  coordinates.  They  are  displayed  in  the  current  unit  type. 
The  units  parameter  displays  the  current  unit  type  (canonical,  metric  or  english). 
It  is  initially  set  to  canonical.  All  system  parameters  can  be  changed  by  the  user. 
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3 .  Parameter  Input 

A  parameter  in  the  Parameter  Window  is  selected  by  placing  the  cursor 
over  it  and  pressing  the  middle  mouse  button.  If  it  can  be  changed,  a  dialog  box 
or  menu  will  appear.  Dialog  boxes  are  answered  with  keyboard  entries.  All 
keyboard  entries  must  be  entered  in  the  current  unit  type  (e.g.,  metric).  A  menu 
is  answered  by  selecting  an  entry  and  pressing  the  right  button.  Parameters  that 
toggle  between  "On"  and  "Off"  will  change  with  no  other  input  required. 
Similarly,  the  units  parameter  will  cycle  through  "Canonical",  "Metric"  and 
"English"  each  time  it  is  selected.  Certain  parameters  are  computed  by  SMET. 
If  one  of  these  parameters  is  selected,  the  user  is  told  that  it  cannot  be  changed. 

SMET  also  allows  the  user  to  adjust  certain  parameters  with  five  dials 
on  the  dial  box.  These  include  the  clock  speed  and  four  of  the  classical 
parameters.  Figure  5-6  shows  how  the  dials  are  arranged.  The  dials  are  used 
for  user  convenience  and  simple  demonstrations  of  the  manipulation  of 
parameters.  If  no  dial  box  is  available,  these  variables  can  still  be  input  by  using 
the  Parameter  Window.  Whenever  a  dial  is  turned,  the  appropriate  parameter 
will  be  changed.  The  clock  speed  dial  is  initially  calibrated  to  operate  between 
zero  and  six  hours.  If  a  clock  speed  greater  than  six  hours  is  desired,  the  user 
can  enter  it  by  using  the  Parameter  Window. 

A  change  to  a  satellite's  classical  elements  will  radically  alter  its 
position.  Therefore,  SMET  erases  all  displayed  plots  for  the  Current  Satellite 
and  begins  plotting  again.  Similarly,  if  the  current  time  is  changed,  all  plots  for 
all  satellites  are  erased. 
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D.  SATELLITE  LIST  WINDOW 

The  Satellite  List  Window,  seen  in  Figure  5-7,  lists  all  of  the  satellites  in  the 
system.  The  Current  Satellite  is  listed  at  the  top  of  the  list.  The  user  changes  the 
Current  Satellite  by  placing  the  cursor  over  the  desired  satellite  and  pressing  the 
middle  mouse  button.  The  other  satellites  will  shift  down  one  position  as 
necessary  to  reorganize  the  list. 

E.  MENU  COMMANDS 

When  the  user  presses  the  right  hand  button  on  the  mouse,  a  menu  will 
appear  near  the  position  of  the  cursor.  A  menu  item  is  selected  by  pressing  the 
right  hand  mouse  button  again  while  the  cursor  is  over  a  particular  command  or 
option.  Certain  menu  items  will  have  an  arrow  to  the  right.  These  items 
indicate  that  another  menu  (a  rollover  menu)  will  appear  if  the  cursor  is  moved 
to  the  right  of  that  menu  item.  Pressing  the  right  hand  mouse  button  while  it  is 
over  the  menu  title  or  at  any  time  that  it  is  not  over  a  menu  item  will  abort  the 
process  of  making  a  menu  selection. 

The  main  menu  allows  the  user  to  add  a  satellite,  maneuver  the  current 
satellite,  delete  the  Current  Satellite,  save  satellites’  parameters  into  a  file,  take  a 
picture,  choose  the  set  of  satellite  parameters  to  display  in  the  Parameter 
Window,  list  the  user  instructions  or  exit  the  program. 

1.  Add  Satellite 

The  Add  Satellite  command  has  three  options  that  appear  in  a  rollover 
menu.  The  user  can  define  an  orbit's  parameters,  define  the  parameters  for  a 
launch,  or  read  a  file  of  previously  saved  satellites. 
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Figure  5-6.  Dial  Box 
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SMET  Satellite  List 

MoJniya  Orbit 
Low  Earth  Ott*bit 
Geostationary  Orbit 
Geosynchronous  Orbit 
Another  Orbit 
Yet  Another  Orbit 


Figure  5-7.  Satellite  List  Window 
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a.  Define  Orbit 

Defining  an  orbit's  parameters  also  has  a  rollover  menu.  Here  the 
user  is  given  the  option  of  selecting  from  three  pre-defined  orbits  (a  Low  Earth 
orbit  (300  NM),  a  Molniya  orbit,  and  a  Geosynchronous  orbit).  The  user  also 
can  select  "Input  Parameters".  This  allows  the  user  to  enter  the  classical 
parameters  of  the  desired  orbit  through  a  series  of  dialog  boxes.  Parameters 
must  be  entered  in  the  current  unit  type  (e.g.  km,  if  in  metric).  The  user  is 
asked  to  enter  a  name  for  the  new  satellite  in  all  cases.  It  will  be  given  a 
randomly  selected  color  and  have  plotting,  fill,  and  vectors  turned  off.  The  new 
satellite  becomes  the  Current  Satellite. 

b.  Define  Launch 

The  Define  Launch  command  also  initiates  a  series  of  dialog  boxes. 
The  user  is  asked  to  enter  a  name  for  the  new  satellite  and  burnout  parameters 
related  to  the  launch.  These  parameters  include  the  time  at  which  burnout  occurs 
(not  tc  be  confused  with  the  length  of  the  rocket  bum),  the  altitude  above  the 
surface  of  the  earth  (in  the  current  unit  type,  e.g.  km),  the  velocity  (in  the 
current  unit  type,  e.g.  km/sec),  the  earth-fixed  burnout  latitude  and  longitude 
(±degrees),  the  flight  path  angle  and  the  launch  azimuth.  If  the  time  of  burnout 
is  earlier  than  the  current  time,  the  satellite  will  be  placed  at  a  position  where  it 
would  have  been  had  the  launch  been  executed  at  the  appropriate  time.  If  not,  a 
launch  record  is  saved  until  the  appropriate  time  of  execution.  When  executing  a 
launch,  SMET  will  add  in  the  appropriate  earth  rotation  correction.  The  new 
satellite  will  be  given  a  randomly  selected  color  and  have  plotting,  fill,  and 
vectors  turned  off.  It  will  become  the  Current  Satellite.  See  Chapter  III  for 
more  information. 
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c .  Read  File 

The  Read  File  command  will  read  a  satellite  file  previously  saved 
using  the  Save  Satellites  command  discussed  below. 

2 .  Maneuver  Satellite 

The  Maneuver  Satellite  command  initiates  a  sequence  of  dialog  boxes. 
The  user  is  asked  to  enter  parameters  for  a  maneuver  of  the  Current  Satellite. 
The  first  parameter  is  the  true  anomaly  at  which  the  maneuver  should  occur.  If 
there  are  any  maneuvers  pending  for  that  satellite,  the  user  will  be  asked  to  enter 
the  number  of  degrees  of  travel  following  a  previous  maneuver.  When  a 
previous  maneuver  is  executed,  this  quantity  will  be  added  to  the  true  anomaly  at 
the  end  of  the  previous  maneuver.  For  example,  if  900  is  entered  the  new 
maneuver  will  occur  two  and  a  half  revolutions  after  the  previous  maneuver 
occurs.  There  is  no  limit  to  the  number  of  maneuvers  which  can  be  added  to  a 
satellite's  queue.  Maneuvers  will  be  executed  in  the  order  in  which  they  are 
entered.  A  satellite  does  not  have  to  be  the  Current  Satellite  at  the  time  of  the 
maneuver. 

After  the  first  entry  the  user  is  asked  to  choose  between  two  methods  of 
entering  the  parameters  of  a  maneuver.  The  first  option  is  to  enter  the  magnitude 
of  the  velocity  change  along  with  an  associated  pitch  and  yaw.  The  user  enters 
the  magnitude  in  the  current  units  (e.g.,  km/sec).  The  pitch  and  yaw  are  entered 
in  degrees.  Pitch  is  a  right-hand  rotation  about  the  satellite’s  z-axis  (thumb 
pointing  towards  the  earth).  Yaw  is  a  right-hand  rotation  about  the  satellite’s 
y-axis  (thumb  pointing  opposite  to  the  orbital  plane's  normal  vector).  SMET 
uses  the  magnitude,  yaw  and  pitch  of  the  velocity  change  to  compute  the  changes 
to  the  magnitude  of  the  velocity,  flight  path  angle  and  heading  angle  of  the 
satellite  at  the  time  of  the  maneuver.  The  second  option  for  entering  a  maneuver 
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is  to  enter  the  changes  to  the  magnitude  of  the  velocity,  flight  path  angle  and 
heading  angle  directly.  This  option  is  useful  to  users  who  are  more  interested  in 
the  affects  of  a  maneuver  than  in  the  total  Av  required.  See  Chapter  IV  for  more 
information. 

3.  Delete  Satellite 

The  Delete  Satellite  command  will  remove  the  Current  Satellite  from 
the  system  permanently.  The  second  satellite  listed  in  the  Satellite  List  Window 
becomes  the  Current  Satellite. 

4.  Save  Satellites 

The  Save  Satellites  command  will  first  ask  the  user  for  a  file  name.  The 
user  may  enter  any  Unix  path  and  file  name.  If  the  file  name  is  not  unique,  the 
old  file  will  be  overwritten.  The  user  will  then  be  asked  for  a  number  of 
satellites  to  be  saved.  Satellites  are  saved  from  in  the  order  that  they  appear  in 
the  Satellite  List  Window.  If  the  user  only  desires  a  specific  subset  of  the 
satellites  in  the  system,  each  satellite  desired  should  first  be  selected  in  sequence 
until  they  are  at  the  top  of  the  list.  Any  file  saved  using  this  procedure  can  be 
read  by  Read  File  command  discussed  above.  Saving  satellites  does  not  remove 
them  from  the  system. 

5.  Take  Picture 

The  Take  Picture  command  has  four  options  that  will  appear  in  a 
rollover  menu.  All  options  request  a  file  name.  The  first  two  options  will  save 
an  sgi-rgb  formatted  file.  These  files  can  be  shown  and  manipulated  with  several 
image  tools  available  on  the  IRIS.  These  image  tools  include  showsgi,  shrinksgi 
and  others.  They  are  available  in  the  imagetool  subdirectory  of  SMETs 
directory.  The  first  option  saves  the  entire  screen  as  it  appears.  The  second 
saves  the  entire  screen  after  the  background  has  been  changed  to  white  and  white 
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items  have  been  changed  to  gray.  The  third  and  fourth  options  will  save  the 
entire  screen  or  the  Main  Window  in  Encapsulated  PostScript  (EPS)  format. 
These  files  can  be  printed  on  a  PostScript  laserwriter  or  read  by  any  software 
capable  of  reading  EPS  images.  SMET  will  remain  frozen  for  30  to  90  seconds 
while  the  image  is  being  saved. 

6.  Satellite  Parameters 

The  Satellite  Parameters  command  has  a  rollover  menu.  This  menu 
allows  the  user  to  select  which  of  the  five  sets  of  satellite  parameters  will  be 
displays  1  in  the  Parameter  Window. 

7.  Instructions 

The  Instructions  command  will  temporarily  halt  SMET  and  display  the 
user  instructions. 

8.  Exit 

The  Exit  command  will  quit  SMET.  The  user  is  offered  an  opportunity 
to  save  any  satellites  desired.  Any  satellites  in  the  system  that  have  not  been 
saved  will  be  lost. 
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VI.  PROGRAM  MODIFICATIONS 


A.  GENERAL 

There  are  several  dummy  function  calls  entered  into  SMETs  source  code. 
These  dummy  procedures  can  be  used  to  make  minor  modifications  to  SMET 
without  making  major  changes  to  the  code.  This  capability  is  directed  primarily 
at  thesis  students  requiring  a  specific  simulation.  All  of  the  functions  described 
below  are  included  in  the  file  smet_user.c  which  is  included  in  Appendix  C.  This 
file  can  be  modified  by  novice  C  programmers  for  desired  modifications  to 
SMET.  Other  smet_xxx.c  files  should  not  be  changed  by  inexperienced 
programmers.  Modified  versions  of  SMET  should  have  the  name  altered  to 
prevent  confusion  with  the  original  SMET. 

B.  USER  PROGRAM  INITIALIZATION  ROUTINE 

The  procedure  user_init  can  be  used  to  read  in  satellite  files  and  initialize 
program  variables  to  values  different  than  the  default  value. 

C.  USER  SATELLITE  INITIALIZATION  ROUTINE 

The  procedure  user_new_sat  can  be  used  to  modify  any  satellite  variables 
to  values  different  than  the  default  value.  The  SMET  satellite  structure  provides 
four  unused  variables  (user_floatl,  user_float2,  user_intl  and  user_int2)  which 
can  be  used  in  a  modified  program  for  such  things  as  orbital  perturbations  or 
ASAT  simulations. 
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D.  USER  PROGRAM  UPDATE  ROUTINES 

SMET  provides  two  procedures  that  can  be  used  to  check  for  an  exit 
condition  or  a  condition  for  executing  a  launch  or  maneuver.  Procedures 
userjpre_update_satellites  and  user__post_update_satellites  can  also  be 
used  to  check  the  elapsed  time  or  to  compare  the  positions  of  the  satellites.  One 
routine  is  called  before  all  of  the  satellites  are  updated  and  the  other  is  called 
after. 

E.  USER  SATELLITE  UPDATE  ROUTINES 

SMET  provides  two  procedures  that  can  be  used  to  perturb  an  orbit  or 
delete  a  satellite  based  on  a  given  set  of  parameters.  The  procedure 
user_pre_orbital_parameters  is  called  before  the  satellite's  changing  orbital 
parameters  are  updated.  The  procedure  user_post_orbital_parameters  is 
called  after  the  changing  parameters  have  been  updated.  Both  procedures  take  a 
satellite  and  the  time  as  arguments. 
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VII.  SUMMARY,  RECOMMENDATIONS  AND  FUTURE 

DIRECTIONS 


SMET  provides  beginning  students  with  the  ability  to  learn  the  difficult 
concepts  of  orbital  mechanics  with  greater  ease.  Satellite  maneuvers,  which  are 
not  covered  in  introductory  courses  in  great  detail,  can  now  be  taught  with 
simple  presentations.  Additionally,  a  means  of  generating  simulations  is 
available  for  more  advanced  students. 

In  order  to  achieve  the  maximum  benefit  of  this  project,  the  Space  Systems 
Academic  Group  should  purchase  an  IRIS  graphics  workstation  capable  of 
running  SMET.  Several  students  in  the  Space  Systems  Operations  Curriculum 
have  already  received  instruction  in  programming  the  IRIS.  Already,  one 
student  is  creating  a  smet_user  module  for  a  SMET  based  ASAT  simulation.  As 
a  start,  modules  could  be  written  for  standard  transfer  orbit  maneuvers  or 
perturbations  not  currently  modeled.  Space  Systems  Engineering  students  could 
also  benefit  from  the  use  of  a  CAD  system  or  other  engineering  software 
available  for  the  IRIS.  As  an  interim  measure,  the  Space  Systems  Academic 
Group  should  negotiate  limited  access  to  an  IRIS  belonging  to  another 
department  for  the  development  of  class  presentations. 

It  is  hoped  that  the  software  will  be  used  elsewhere  and  continue  to  evolve. 
Future  enhancements  not  currently  provided  include  the  ability  to  read 
NAVSPASUR  Charlie  Element  formatted  files  for  entering  satellites,  an 
illumination  model  to  represent  the  sun,  a  more  accurate  representation  of 
perturbations  and  a  capability  to  automatically  compute  orbit  transfer  or 
intercept  requirements. 
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APPENDIX  A  ■  COMMON  VARIABLES 


Variable 

Name 

Formula 

Name 

Program 

Name 

Semi-Major 

Axis 

a 

major 

Description  of  conic  section 

imsmm 

e 

ecc 

Description  of  conic  section 

Semi-Latus 

Rectum 

p 

P 

Description  of  conic  section 

Range 

r 

range 

Distance  from  the  center  of  the  earth  to  the 
satellite 

Velocity 

V 

velocity 

Magnitude  of  the  velocity  of  the  satellite 

Inclination 

i 

incl 

Angle  between  the  plane  of  the  equator  and 
the  plane  of  the  orbit 

Longitude  of 
the  Ascending 
Node 

h 

asc 

Angle  in  the  equatorial  plane  between  the 
right  ascension  of  the  sun  at  Vernal  Equinox 
and  the  right  ascension  of  the  satellite 

uSollMm 

g 

peri 

Angle  in  the  orbit  plane  from  the  line  of 
nodes  to  perigee 

mm 

0 

anom 

Angle  in  the  orbit  plane  between  perigee  and 
the  satellite 

Eccentric 

Anomaly 

E 

E 

Angle  in  the  orbit  plane  between  perigee  and 
a  "projection"  of  the  satellite’s  position  onto  a 
circumscribed  auxiliary  circle 

mm 

M 

M 

The  product  of  the  mean  motion  and  elapsed 
tirr 

Epoch 

to 

Time  of  perigee  passage 

Heading 

Angle 

V 

Angle  in  the  local  horizon  plane  between  the 
line  through  the  satellite  and  the  North  point 
and  the  line  of  motion  measured 

Flight  Path 

Y 

f_path 

Angle  in  the  orbit  plane  measured  upward 
from  the  local  horizon  plane  to  the  line  of 
motion 

(0 

Angle  in  the  equatorial  plane  between  the  line 
of  nodes  and  the  meridian  of  the  satellite 

Argument  of 
Latitude 

u 

u 

Angle  in  the  orbit  plane  between  the  line  of 
nodes  and  the  satellite 

Latitude 

iat 

Latitude  of  the  subsatellite  point 

Longitude 

m 

Ion 

Earth-fixed  Longitude  of  the  subsatellite 
point 

Gravitational 

Parameter 

Constant  for  computation  of  earth  orbits 
(1  DUffi3/TU©2  =  1.407647  x  1016  ft3/Sec2) 

50 


APPENDIX  B  -  VIDEO  RECORDING  AND  PRINTING 


The  following  instructions  apply  to  the  implementation  of  SMET  in  the 
Graphics  and  Video  Laboratory  located  in  the  Computer  Science  Department  of 
the  Naval  Postgraduate  School. 

Video  Recording  -  With  Professor  Zyda's  permission,  an  S-VHS  or  a 
standard  VHS  tape  can  be  recorded  while  running  SMET  on  IRIS1.  This  tape 
can  then  be  edited  and  dubbed  as  necessary  for  presentations.  Note  that  some 
resolution  will  be  lost  because  standard  video  signals  have  a  lower  resolution 
than  the  IRIS.  The  top  and  bottom  edges  of  the  image  will  also  be  clipped 
because  of  the  different  aspect  ratios.  The  YEM  scan  converter  should  be  turned 
on  in  order  to  record.  The  video  recorder  will  record  everything  appearing  on 
the  video  monitor. 

Printing  EPS  Files  -  EPS  files  are  saved  using  one  of  the  EPS  options  of 
the  Take  Picture  command.  These  files  can  be  printed  on  any  PostScript  laser 
printer.  The  file  can  be  printed  with  the  following  command:  "Ipr  -Pps2 
filename".  The  file  can  also  be  imported  into  other  programs  such  as 
Framemaker  for  incorporation  into  other  documents  as  was  done  in  the 
preparation  of  this  document. 
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APPENDIX  C  -  SMET  PROGRAM  LISTING 


SMET  consists  of  the  following  files: 

smet.h  (p.  54  -  59)  contains  all  constants  and  structures  used  by 
SMET.  Also  has  all  forward  function,  procedure  and  global 
variable  declarations. 

smet_vars.c  (p.  60  -  64)  defines  all  global  variables. 

smet_addelete.c  (p.  65  -  68)  includes  procedures  to  add,  delete 
and  switch  satellites.  Additionally,  procedures  to  add  and 
delete  plots  to  satellite  structures  are  included. 

smet_axes.c  (p.  69)  includes  procedure  to  draw  the  axes  in  the 
OmniView  window. 

smet_dialog.c  (p.  70  -  72)  includes  procedures  for  dialog  and 
message  boxes. 

smet_files.c  (p.  73  -  77)  includes  procedures  that  read  and  write 
satellite  files,  as  well  as  reading  SMET  data  files  and  writing 
image  files. 

smet_globe.c  (p.  78  -79)  includes  procedures  to  draw  the  globe  in 
the  OmniView  and  Satellite  View  windows. 

smet_init.c  (p.  80  -  82)  includes  SMET  initialization  procedures. 

smet_kepler.c  (p.  83  -  85)  includes  procedures  that  update  a 
satellite's  orbital  parameters. 

smet_launch.c  (p.  86  -  87)  includes  procedures  to  launch  a 
satellite. 
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smet_main.c  (p.  88  -  89)  is  the  heart  of  SMET.  It  includes  the 
main  execution  loop. 

smet_maneuver.c  (p.  90  -  95)  includes  procedures  to  maneuver  a 
satellite. 

smet_map.c  (p.  96  -97)  includes  procedures  to  draw  the  map  in 
the  Mercator  View  window. 

smet_orbit.c  (p.  98)  includes  procedure  to  draw  a  satellite’s 
orbital  path  in  the  OmniView  window. 

smet_param.c  (p.  99  -  104)  includes  procedures  that  handle 
Parameter  window  events. 

smet_queue.c  (p.  105  -  109)  includes  procedures  read  the  queue 
and  handle  events  not  covered  in  other  files. 

smet_sat.c  (p.  110  -  113)  includes  procedures  to  draw  a  satellite 
in  the  OmniView  window. 

smet_user.c  (p.  114)  includes  dummy  procedures  discussed  in 
Chapter  VI. 

smet_windows.c  (p.  115  -  123)  includes  procedures  that  draw  the 
various  windows. 

makefile  (p.  124)  can  be  used  to  compile  SMET  in  the  NPS 
Computer  Science  Department  Graphics  and  Video 
Laboratory.  Files  in  the  imagetools  sub-directory  will  also  be 
required. 

smet_data.d  (not  included)  is  a  binary  data  file  necessary  to  draw 
the  globe  and  map. 
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it .h  Mon  Sep  3  17:59:00  1990  1 


/**********************************/ 
/*  */ 

/ *  smet . h  *  / 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  global  variables,  constants  and  */ 

/*  structures  used  by  the  SMET  demonstrator.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/**********************************/ 

♦include  <math.h>  /*  math  functions  file  */ 

♦include  <string.h>  /*  string  routines  */ 

♦include  "gl.h"  /*  graphics  library  file  */ 

♦include  "device. h"  /*  device  header  files  */ 

♦include  "stdio.h"  /*  I/O  routines  */ 

♦include  "fmclient.h"  /*  IRIS  font  manager  */ 

♦include  "imagesupport/image_types . h"  /*  for  picture  saving  */ 

♦include  "imagesupport/image_funcs .h"  /*  for  picture  saving  */ 


♦define  FINEPOINTS  12843 
♦define  COARSEPOINTS  3322 
♦define  BORDERPOINTS  1407 

♦define  OMNISATMERC  1  /*  Which  view  is  in  which  window?  */ 

♦define  OMNIMERCSAT  2  /*  These  constants  used  by  switch_windows ( )  */ 

♦define  SATMERCOMNI  3 

♦define  SATOMNIMERC  4 

♦define  MERCOMNISAT  5 

♦define  MERCSATOMNI  6 

♦define  EARTH  9  /*  Color  Codes  */ 

♦define  COAST  9 

♦define  BORDER  10 

♦define  LATLONG  11 

♦define  SUN  12 

♦define  SAT_RED  13 

♦define  SAT_GREEN  14 
♦define  SAT_YELLOW  15 
♦define  SAT_BLUE  16 
♦define  SAT_MAGENTA  17 
♦define  SAT_CYAN  18 
♦define  SAT_WHITE  19 
♦define  SAT_GRAY  20 
♦define  MAX  COLORS  21 


♦define  CANONICAL  0  /*  Unit  Codes  */ 

♦define  METRIC  1 

♦define  ENGLISH  2 


♦define  DISTANCE  1  /* 


♦define  SPEED  2 
♦define  DEGREE  3 
♦define  LAT  4 
♦define  LON  5 
♦define  TIME  6 
♦define  FLAG  7 
♦define  FLOAT  8 
♦define  STRING  9 


♦define  COLORCODE  10 
♦define  UNITCODE  11 


for  parameter  type  field  */ 
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SOMt.h 

Mon  Sep  3 

17:59:00  1990 

2 

♦define 

ROT AT I ON_RAT E 

4. 37526 95 128 6e -3 

/* 

earth  rate 

in  rads/min 

V 

♦define 

TU 

13.44686457 

/* 

minutes  per  TU 

*/ 

♦define 

DU_METRIC 

6378.1459997 

/* 

KM  per  DU 

*/ 

♦define 

DU  ENGLISH 

3443.9227864 

/* 

NM  per  DU 

V 

♦define 

DU_PER  TU_METRIC 

7.9053682830 

/* 

KM/ SEC  per 

DU/TU 

*/ 

♦define 

DU_PER  TU_ENGLISH 

25936.24764 

/* 

FT/SEC  per 

DU/TU 

*/ 

♦define 

PI 

3.14159265359 

/* 

Not  mom' s  apple  . . . 

*/ 

♦define 

TO_DEG 

57.2957795132 

/* 

Convert  degrees  to  radians 

*/ 

♦define 

TO  RAD 

1. 74532925199e-2 

/^Convert  degrees  to  radians 

V 

♦define 

DEG30 

0.523598775598 

/* 

30  degrees  in  radians 

V 

♦define 

DEG60 

1.04719755120 

/* 

60  degrees  in  radians 

*/ 

♦define 

DEG90 

1.57079632679 

/* 

90  degrees  in  radians 

*/ 

♦define 

DEG120 

2.09439510239 

/* 

120  degrees  in  radians 

V 

♦define 

DEG180 

3.14159265359 

/* 

180  degrees  in  radians 

*/ 

♦define 

DEG270 

4.71238898038 

/* 

270  degrees  in  radians 

*/ 

♦define 

DEG360 

6.28318530718 

/* 

360  degrees  in  radians 

*/ 

♦define 

J2 

1 . 08264e-3 

/* 

perturbation 

constant 

*/ 

♦define 

SUN  INCL 

-23 . 454 *TO  RAD 

/* 

pseudo-incl 

for  plots 

V 

♦define 

SUN  TU 

83711.26224 

/* 

minutes  per 

Solar  TU 

*/ 

♦define  CLASSICAL  1  /*  parameter  display  options  */ 

♦define  SECONDARY  2 
♦define  POSITIONS  3 
♦define  OPTIONS  4 
♦define  COMBINED  5 

♦define  LETTER  18  /*  font  size  and  line  spacing  */ 

♦define  LINE  22 

♦define  ERROR  -9999 

♦define  max (x,  y)  (x  >  y)  ?  x  :  y  /*  return  maximum  value  */ 

♦define  nonzerop  (x)  fabs  (x)  >  0.00001  /*  ~  float  !  =  0.0  */ 

♦define  zerop(x)  f abs (x)  <  0.00001  /*  ~  float  ==0.0  */ 
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•Mt.h  Mon  Sep  3  17:59:00  1990  3 

typedef  struct  { 

float  x [ FINEPOINTS ] ; 
float  y [FINEPOINTS] ; 
float  z [ FINEFOINTS ] ; 
int  pen [FINEPOINTS] ; 

)  fine_data; 


typedef  struct  { 

float  x [COARSEPOINTS] ; 
float  y [COARSEPOINTS] ; 
float  z [COARSEPOINTS] ; 
int  pen [COARSEPOINTS ] ? 
)  coarse  data; 


typedef  struct  { 

float  x [BORDERPOINTS] ; 
float  y [ BORDERPOINTS ] ; 
float  z [BORDERPOINTS] ; 
int  pen [ BORDERPOINTS ] ; 
)  border_data; 


typedef  struct  ( 

char  ‘name; 
int  type; 
char  * value; 
I  ear; 


/*  Displayable  parameters  */ 

/*  for  dialog  boxes  */ 

/*  DISTANCE,  SPEED,  TIME  */ 
/*  cast  this  to  whatever  */ 


typedef  struct  { 

int  numparams; 
Par  *params; 

)  ParBlock; 


/*  Block  of  parameters  */ 

/*  number  of  parameters  */ 
/*  list  of  parameters  */ 


typedef  struct  (  /* 

int  orgx,orgy, 
sizex, sizey, 
gid; 

]  window_struct; 


Window 

structure 

V 

/* 

lower  left  corner 

*/ 

/* 

window  size 

*/ 

/* 

window  grahics  ID 

V 

typedef 


struct 

launch_st  ruct 

{  /* 

Launch  structure 

*/ 

char 

name  [30] ; 

/* 

satellite  name 

*/ 

float 

time. 

/* 

when  burnout  will  occur 

*/ 

v. 

/* 

burnout  velocity 

V 

alt. 

/* 

burnout  altitude 

*/ 

lat. 

/* 

earth-fixed  latitude 

*/ 

Ion, 

/* 

earth-fixed  longitude 

*/ 

elev. 

/* 

elev  angle  -  above  horizon 

*/ 

azim; 

/* 

azimuth  angle 

*/ 

struct 

launch  struct 

*next 

;/*  the  next  launch 

V 

)  launch_struct  ; 
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Mon  Sap  3  17:59:00  1990 


typedef  struct  plot_struct  {  /*  Plot  structure  */ 
float  lat,lon;  /*  a  plot 

struct  plot_struct  *next;  /*  the  next  plot 
)  plot_struct; 


struct 

man_struct  {  /* 

Maneuver  structure 

*/ 

int 

no  rma l_ma ne u v e r 

;/* 

flag  for  deltav, yaw, pitch 

*/ 

float 

time. 

/* 

when  maneuver  will  occur 

*/ 

anom. 

/* 

true  anomaly  of  occuran ce 

*/ 

deltav. 

/* 

Change  in  velocity 

*/ 

yaw. 

/* 

yaw  or  heading  change 

*/ 

pitch; 

/* 

pitch  or  flight  path  change 

*/ 

struct 

man  struct  *next 

;/* 

the  next  maneuver 

*/ 

)  man_ 

struct ; 

struct 

satellite  { 

/* 

Satellite  structure 

*/ 

char 

name [30] ; 

/* 

satellite  name 

*/ 

float 

major. 

/* 

semi-major  axis  of  orbital 

ellipse 

*/ 

minor. 

/* 

semi-minor  axis  of  orbital 

ellipse 

*/ 

ecc. 

/* 

eccentricity  of  orbital  ellipse 

*/ 

incl , 

/* 

inclination  of  orbital  ellipse 

*/ 

asc. 

/* 

longitude  of  the  ascending 

node 

*/ 

peri. 

/* 

argument  of  periapsis 

*/ 

anom. 

/* 

true  anomaly  at  epoch 

*/ 

epoch. 

/* 

time  of  epoch  +ve  (if  after 

0000Z) 

*/ 

E, 

/* 

eccentric  anomaly  at  epoch 

*/ 

P- 

/* 

X-coordinate  in  orbital  plane 

*/ 

q. 

/* 

Y-coordinate  in  orbital  plane 

*/ 

X, 

/* 

X-coordinate  in  IJK  system 

*/ 

/* 

Y-coordinate  in  IJK  system 

*/ 

2, 

/* 

Z-coordinate  in  IJK  system 

*/ 

range. 

/* 

Sat  distance  from  center  of 

earth 

*/ 

apogee  range. 

/* 

Sat  distance  at  apogee 

*/ 

perigee  range. 

/* 

Sat  distance  at  perigee 

*/ 

lat , 

/* 

Satellite  latitude 

*/ 

Ion, 

/* 

Satellite  longitude 

*/ 

period. 

/* 

Satellite  period  (one-revolution) 

*/ 

velocity. 

/* 

current  satellite  velocity 

*/ 

asc  change. 

/* 

pertrubation  constant 

*/ 

peri_change. 

/* 

pertrubation  constant 

*/ 

user_f loatl. 

/* 

for  program  modifications 

*/ 

user_f loat2. 

/* 

for  program  modifications 

*/ 

matrix [4] [4] ; 

/* 

orbit's  rotation  matrix 

*/ 

int 

plot. 

/* 

are  we  plotting  ? 

*7 

fill. 

/* 

do  we  fill  in  the  orbit  ellipse? 

*/ 

vectors. 

/* 

do  we  draw  vectors  on  the  e 

llipse? 

*  / 

color. 

/* 

display  color 

*/ 

user  inti. 

/* 

for  program  modifications 

*/ 

user  int2; 

/* 

for  program  modifications 

*/ 

plot_struct  *plots;  /* 
man_struct  ‘maneuver;  /* 
struct  satellite  *next;/* 
)  satellite; 


plots  since  plotting  started 

next  maneuver 

pointer  to  next  satellite 
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SflMt  .h 


Mon  Sep  3  17:59:00  1990 


5 


/*  global  variable 

declarations 

*/ 

extern 

satellite  header  sat. 

/* 

Header  to  sat  list 

*/ 

*current_sat 

;  /* 

Current  satellite 

*/ 

extern 

launch_struct  Spending 

_launches; /* 

List  of  upcoming  launches 

V 

extern 

float  clock, 

/* 

Time  elapsed 

*/ 

clock_speed. 

/* 

time  per  minutes 

*/ 

viewer^x. 

/* 

Viewer  X  pos  in  Omniview 

*/ 

viewer_y. 

/* 

Viewer  Y  pos  in  Omniview 

*/ 

viewer^z, 

/* 

Viewer  Z  pos  in  Omniview 

*/ 

viewer^range; 

/* 

Viewer  Range  in  Onniview 

*/ 

extern 

int  top  menu. 

/* 

Menu  ID 

*/ 

addmenu. 

/* 

Menu  ID 

*/ 

orbit  menu. 

/* 

Menu  ID 

V 

color  menu. 

/* 

Menu  ID 

*/ 

picture  menu. 

/* 

Menu  ID 

*/ 

param  menu. 

/* 

Menu  ID 

*/ 

window  order. 

/* 

which-where  ex.  OMNISATMERC 

*/ 

omni  view  gid. 

/* 

gid  of  OMNI  view  window 

*/ 

sa t_view_gid. 

/* 

gid  of  SAT  view  window 

*/ 

mere  view_gid. 

/* 

gid  of  MERC  view  window 

*/ 

unit_type. 

/* 

CANONICAL,  METRIC,  AMERICAN 

*/ 

param  option. 

/* 

Which  satellite  parameters? 

*/ 

clock_running. 

/* 

is  the  clock  running? 

*/ 

f ineflag. 

/* 

Detailed  Earth  drawing? 

*/ 

borderilag. 

/* 

Draw  political  borders? 

*/ 

auto  viewer. 

/* 

Auto  pos  for  Omni  view? 

*/ 

perturbations ; 

/* 

Do  we  do  perturbations? 

V 

/*  Globe  and  map  data  records  */ 

extern  fine_data  fine2d_data,  fine3d_data; 
extern  coarse_data  coarse2d__data,  coarse3d_data; 

extern  border_data  border2d_data,  border3d_data; 

/*  parameter  block  containing  all  displayable  parameters  */ 

extern 

ParBlock  classicalparamblock, 
secondaryparamblock, 
posit ionparamblock, 
optionsparamblock, 
combinedparamblock, 
sysparamblock; 
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s»et  .h 


Mon  Sep  3  17:59:00  1990 


6 


extern 

char  ^instructions [28] , 

*color_names lMAX_COLORS] 
*unit_names [3] , 
*speed_unit [3] , 
*distance_unit [3] ; 

extern 

float  speed_f actor [3] , 

distancef actor [3 ] ; 

extern  /* 

window_struct  backwin, 
mainwin, 
auxwinl, 
auxwin2, 
satlistwin, 
dialogwin, 
paramwin; 


/* 

User 

instructions 

*/ 

,  /* 

Nai.ies 

of 

all 

the  colors 

*/ 

/* 

Names 

of 

the 

Units 

*/ 

/* 

Names 

of 

the 

speed  units 

*/ 

/* 

Names 

of 

the 

distance  units 

*/ 

/*  Speed  conversion  factor 

*/ 

/*  Distance 

conversion  factor 

»/ 

up  all  our 

window  parameters 

*/ 

extern 

float  identjmatrix[4]  [4]  , 

Color [MAX  COLORS ]  [ 4 ]  ; 


/*  Nice  constant  to  have  laying  around  */ 
/*  Array  of  colors  for  c3f  or  c4f  */ 


/*  forward  function  declarations  */ 

draw_instructions () ; 
add  launch () ; 
ask_and_read_sat_f ile () ; 
ask_and_write_sat_f ile () ; 
add_maneuver ( ) ; 
add_satellite (int) ; 
custom_3atellite ()  ; 
delete_3atellite ()  ; 
take_picture {)  ; 
program_exit ()  ; 
smet_param_option (int) ; 
smet_sat_color (int) ; 

float  ask  time(); 
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/********************************* 

/* 

/*  smet_vars.c 

/*  by  Carlos  I.  Noriega 

/* 

/*  This  file  contains  the  global  variables,  constants  and 

/*  structures  used  by  the  SMET  demonstrator. 

/* 

/*  Written  -  September  3,  1990 

/*  Modified  -  September  3,  1990 

/* 

/********************************* 


*/ 

V 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 

*/ 


♦include  "smet.h" 


satellite  header  sat. 

/* 

Header  to  sat  list 

*/ 

*current_sat  —  NULL; 

/* 

Pointer  to  list  of  sats 

*/ 

launch  struct  'pending  launches 

~  NULL; /*  Upcoming  launches 

*/ 

float 

clock  «  0.0, 

/* 

Time  elapsed  in  minutes 

*  / 

clock  speed  -  30.0, 

/* 

Time  per  cycle  in  minutes 

*/ 

viewer  x, 

/* 

Viewer  X  pos  in  Omniview 

*/ 

viewer  y. 

/* 

Viewer  Y  pos  in  Omniview 

*/ 

viewer_z. 

/* 

Viewer  Z  pos  in  Omniview 

*/ 

viewer  range; 

/* 

Viewer  Range  in  Omniview 

*/ 

int 

top  menu. 

/* 

Menu  ID 

*/ 

add  menu. 

/* 

Menu  ID 

*/ 

orbit  menu. 

/* 

Menu  ID 

*/ 

picture_menu. 

/* 

Menu  ID 

*/ 

color  menu. 

!* 

Menu  ID 

*  / 

param  menu. 

/* 

Menu  ID 

*  / 

window  order. 

/* 

which-where  ex.  OMNISATMERC 

*/ 

omni  view  gid. 

/* 

gid  of  OMNI  view  window 

*/ 

sat  view  gid. 

/* 

gid  of  SAT  view  window 

*/ 

mere  view  gid. 

/* 

gid  of  MERC  view  window 

*/ 

unit  type  **  CANONICAL, 

/* 

CANONICAL,  METRIC,  ENGLISH 

*/ 

param  option  *  CLASSICAL, 

/* 

Which  satellite  parameters? 

*/ 

clock  running  =  TRUE, 

/* 

is  the  clock  running? 

*/ 

fineflag  =  TRUE, 

/* 

Detailed  Earth  drawing? 

*/ 

border flag  ”  TRUE, 

/* 

Draw  political  borders? 

*/ 

auto  viewer  =  TRUE, 

/* 

Auto  pos  for  Omni  view? 

*/ 

perturbations  “  FALSE; 

/* 

Do  we  do  perturbations? 

*/ 

/*  Globe  and  map  data  records 

finedata  fine2d_data,  fine3d_data; 
coarse_data  coar se2d_data,  coar se3d_dat a ; 
border  data  border2d  data,  border3d_data; 
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2 


Par  classical_list [ ]  =  ( 

(  "Name  ",  STRING,  (char  *) header_sat . name  }, 

(  "Semi -Major  ",  DISTANCE,  (char  *)& (header_sat .major)  }, 

(  "Semi -Minor  ",  DISTANCE,  (char  *) & (header_sat .minor)  ), 

(  "Eccentricity  ",  FLOAT,  (char  *) & (header_sat . ecc)  ), 

(  "Inclination  ",  DEGREE,  (char  * ) fi (header_sat . incl)  ), 

{  "Ascending  Node  ",  DEGREE,  (char  *) & (header_sat .asc)  ), 

(  "Arg  of  Periapsis  ",  DEGREE,  (char  *) 6 (header_sat . peri)  ), 

(  "True  Anomaly  ",  DEGREE,  (char  *) & <header_sat . anom)  ), 

{  "Epoch  ",  TIME,  (char  *)& (header_sat . epoch)  )  ), 

secondary_list (]  =  ( 

(  "Name  ",  STRING,  (char  *) header_sat . name  ), 

(  "Perigee  ",  DISTANCE,  (char  *) & (header_sat . perigee_range)  ), 
(  "Apogee  ",  DISTANCE,  (char  *) & (header_sat . apogee_range)  ), 

(  "Epoch  ",  TIME,  (char  *) & (header_sat . epoch)  ), 

{  "Period  ",  TIME,  (char  *) & (header_sat . period)  ), 

(  "Range  ",  DISTANCE,  (char  *)& (header_sat . range )  ), 

{  "Velocity  ",  SPEED,  (char  *) & (header_sat . velocity)  )  ), 


position_list  ( ]  *=  ( 


{ 

"Name 

",  STRING, 

(char  *) header 

sat . name 

} 

< 

"X 

n 

t 

DISTANCE, 

(char  *) & (header 

_sat . x) 

}, 

{ 

"  Y 

fi 

t 

DISTANCE, 

(char  *)«,  (header 

_sat .y) 

), 

( 

"Z 

ft 

r 

DISTANCE, 

(char  *)& (header 

_sat . z) 

), 

< 

"P 

r 

r 

DISTANCE, 

(char  *)& (header 

sat .p) 

), 

{ 

"Q 

n 

/ 

DISTANCE, 

(char  *)& (header 

_sat .q) 

)  r 

(  "Latitude  ",  LAT,  (char  *) S (header_sat . lat)  ), 

(  "Longitude  ",  LON,  (char  * ) & (header_sat . Ion)  )  ), 


optionlist [ )  -  ( 

(  "Name  ",  STRING,  (char  * ) header_sat . name  }, 

(  "Color  ",  COLORCODE,  (char  *)& (header_sat . color)  ), 
{  "Fill  ",  FLAG,  (char  * ) & (header_sat . f ill )  ), 

(  "Vectors  ",  FLAG,  (char  *) S (header_sat . vectors)  ), 

(  "Plot  ",  FLAG,  (char  *) & (header_sat .plot)  )  ), 


sysparaml ist [ ]  ”  ( 


< 

{ 

I 

( 

< 

< 

( 

( 

( 

< 

< 


"Clock  ",  FLAG,  (char  *) 4clock_running  ), 

"Current  Time  ",  TIME,  (char  *)&clock  ), 

"Clock  Speed  ",  TIME,  (char  *) Sclock_speed  ), 
"Detailed  Maps  ",  FLAG,  (char  *) Sfineflag  ), 
"Political  Borders  ",  FLAG,  (char  *) Sborderf lag  ), 
"Perturbations  ",  FLAG,  (char  * ) ^perturbations  ), 
"Auto  Omniviewer  ",  FLAG,  (char  * ) &auto_viewer  ), 
"Viewer  X  ”,  DISTANCE,  (char  *) &viewer_x  ), 

"  Y  ",  DISTANCE,  (char  *)4viewer_y  ), 

"  Z  ",  DISTANCE,  (char  *)&viewer_z  ), 

"Units  ",  UNITCODE,  (char  *)&unit_type  ), 
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combined_list f ]  —  ( 

(  "Name  ",  STRING,  (char  *) header_sat . name  }, 

(  "Semi -Major  ",  DISTANCE,  (char  *) 6 (header_sat .ma jor)  }, 

(  "Semi -Minor  ",  DISTANCE,  (char  *) & (header_sat .minor)  >, 

{  "Eccentricity  ",  FLOAT,  (char  *) & (header_sat .ecc)  }, 

(  "Inclination  ",  DEGREE,  (char  *) & (header_sat . incl)  }, 

{  "Ascending  Node  ",  DEGREE,  (char  *) & (header_sat . asc)  ), 

{  "Arg  of  Periapsis  ",  DEGREE,  (char  *) & (header_sat .peri)  ), 

{  "True  Anomaly  ",  DEGREE,  (char  *) & (header_sat .anom)  }, 

(  "Epoch  ",  TIME,  (char  *) t (header_sat .epoch)  ), 

(  "Period  ",  TIME,  (char  *) & (header_sat .period)  }, 

(  "Apogee  ",  DISTANCE,  (char  *) & (header_aat .apogee_range)  }, 

(  "Perigee  ",  DISTANCE,  (char  *) & (header_sat .perigee_range)  }, 
{  "Range  ",  DISTANCE,  (char  *) & (header_sat . range)  ), 

(  "Velocity  ",  SPEED,  (char  *) 4 (header_sat . velocity)  >, 


( 

"X 

*1 

t 

DISTANCE, 

(char 

*> & (header_ 

sat .x) 

) 

( 

"Y 

If 

t 

DISTANCE, 

(char 

*) 4 (header_ 

sat -y) 

) 

< 

"Z 

n 

r 

DISTANCE, 

(char 

*) & (header_ 

sat . z) 

) 

{ 

"P 

n 

f 

DISTANCE, 

(char 

*) S  (header_ 

sat . p) 

) 

l 

"Q 

it 

r 

DISTANCE, 

(char 

*) & (header_ 

sat . q) 

) 

(  "Latitude  ",  LAT,  (char  *) 6 (header_aat . lat)  }, 

{  "Longitude  ",  LON,  (char  * ) & <header_sat . Ion)  ), 

(  "Color  ",  COLORCODE,  (char  *) S (header_sat .color)  ), 
(  "Fill  ",  FLAG,  (char  *) & (header_sat . fill)  ), 

(  "Vectors  ",  FLAG,  (char  *)& <header_sat -vectors)  ), 

(  "Plot  ",  FLAG,  (char  * )  &  (header__sat .  plot )  )  )  ); 


FarBlock  classicalparamblock  -  { 

9, 

classical_list 

>, 

secondaryparamblock  -  { 

7, 

secondary_list 

), 

positionparamblock  —  ( 

8, 

posit ion_li st 

>, 

optionsparamblock  *  ( 

5, 

option_list 

combinedparamblock  —  ( 

25, 

combined_list 

), 

sysparamblock  -  ( 

11, 

sysparamlist 

>; 


/* 

numparams 

*/ 

/* 

Par  list 

*/ 

/* 

numparams 

*/ 

/* 

Par  list 

*/ 

/* 

numparams 

*/ 

/* 

Par  list 

*/ 

/* 

numparams 

*/ 

/* 

Par  list 

*/ 

/* 

numparams 

*/ 

/* 

Par  list 

*/ 

/* 

numparams 

*  / 

/* 

Par  list 

*/ 
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char  "instructions [28 ]  «  /*  User  instructions  */ 

("  SATELLITE  MANEUVER  EVALUATION  TOOL","","","", 

"This  program  is  designed  to  let  the  user  view  how  the  manipulation", 
"of  different  orbital  parameters  affect  the  path  of  a  satellite.", 

"You  will  see  the  orbiting  satellite  from  an  omniscient  point  of", 
"point  of  view,  the  earth  from  the  satellite's  point  of  view  and  a", 
"mercator  projection  map  with  the  track  of  the  satellite  marked.", 

"The  user  can  also  launch  or  maneuver  satellites  by  use  of  a  menu", 
"command. ", 

nn  n  w 

r  t 

" INSTRUCTIONS 

"Rightmouse:  Used  to  pop  up  menus  to  execute  the  desired  command.", 

m  if 

t 

if  n 

9 

"Middlemouse :  Used  to  select  a  parameter,  a  satellite  or  an  ", 

"  Auxiliary  Window.", 

If  If 

f 

"Dials:  The  six  upper  dials  will  allow  the  user  to  adjust  following", 

"  parameters.  The  order  is:  SEMI -MAJOR  AXIS  ECCENTRICITY", 

"  INCLINATION  ASCENDING  NODE", 

"  PERIGEE  CLOCK  SPEED", 

n  n  nn 

f  t 

"Press  the  Rightmouse  button  to  continue."}, 

"colornames [MAX_COLORS]  “  {  "Black",  "Red",  "Green",  "Yellow",  "Blue", 

"Magenta",  "Cyan",  "White", "Earth",  "Coast", 
"Border",  "LatLong",  "Sun", "Red",  "Green", 
"Yellow",  "Blue",  "Magenta",  "Cyan",  "White", 
"Gray" } , 

*unit_names [ 3 ]  =  {  "Canonical",  "Metric",  "English"  }, 

*speec*_unit  [ 3 ]  -  {  "DU/TU",  "km/sec",  "ft/sec"  }, 

*distance_unit [3]  -  {  "DU",  "km",  "NM" } ; 

float  speed__f actor [ 3 ]  «  {  1.0,  DU_PER_TU_METRIC,  DU_PER_TU_ENGLISH  }, 
distancef actor [3]  -  {  1.0,  DU_METRIC,  DU  ENGLISH  ); 
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/*  Set  up  all  our  window  parameters  allowing  for  transfer  */ 

/*  to  a  different  size  screen  system  */ 

window  struct  /*  lower-left  corner,  window-size  and  dummy  G1D  */ 


backwin  =  {  0,  0,  XMAXSCREEN,  YMAXSCREEN,  0  ), 


mainwin 

-  ( 

XMAXSCREEN 

* 

265/1279, 

YMAXSCREEN 

* 

273/1023, 

XMAXSCREEN 

* 

750/1279, 

XMAXSCREEN 

* 

750/1279, 

0 

), 

auxwinl 

-  ( 

XMAXSCREEN 

* 

315/1279, 

YMAXSCREEN 

* 

23/1023, 

XMAXSCREEN 

★ 

250/1279, 

XMAXSCREEN 

* 

250/1279, 

0 

}, 

auxwin2 

-  ( 

XMAXSCREEN 

★ 

715/1279, 

YMAXSCREEN 

* 

23/1023, 

XMAXSCREEN 

it 

250/1279, 

XMAXSCREEN 

* 

250/1279, 

0 

), 

satlistwin  -  (  XMAXSCREEN  *  1015/1279,  YMAXSCREEN  *  0/1023, 

XMAXSCREEN  *  264/1279,  YMAXSCREEN  *  1023/1023,  0  ), 

dialogwin  -  (  XMAXSCREEN  *  270/1279,  YMAXSCREEN  *  620/1023, 

XMAXSCREEN  *  740/1279,  YMAXSCREEN  *  45/1023,  0  ), 

paramwin  =  (  XMAXSCREEN  *  0/1279,  YMAXSCREEN  *  0/1023, 

XMAXSCREEN  *  265/1279,  YMAXSCREEN  *  1023/1023,  0  }; 


float  ident  matrix[4](4]  -  (  (1.0, 

0.0, 

0.0, 

0.0) 

/ 

(0.0, 

1.0, 

0.0, 

0.0) 

f 

(0.0, 

0.0, 

1.0, 

0.0) 

r 

(0.0, 

0.0, 

0.0, 

1.0) 

) 

Color  [MAX_COLOF.3J  [4]  -  {(0.000,  0.000,  0.000,  0.000), 

{1.000,  0.000,  0.000,  0.000), 
(0.000,  1.000,  0.000,  0.000), 
(1.000,  1.000,  0.000,  0.000), 
{0.000,  0.000,  1.000,  0.000), 
(1.000,  0.000,  1.000,  0.000), 
(0.000,  1.000,  1.000,  0.000), 
(1.000,  1.000,  1.000,  0.000), 
(0.250,  0.800,  0.500,  0.000), 
(0.750,  0.200,  0.000,  0.000), 
(0.500,  0.250,  0.000,  0.000), 
(0.750,  0.750,  0.750,  0.000), 
(1.000,  0.800,  0.200,  0.400), 
(0.750,  0.000,  0.000,  0.200), 
(0.000,  1.000,  0.000,  0.200), 
(1.000,  1.000,  0.000,  0.200), 
(0.000,  0.500,  1.000,  0.200), 
(1.000,  0.000,  1.000,  0.200), 
(0.000,  1.000,  ■'.000,  0.200), 
(1.000,  1.000,  1.000,  0.200), 
(0.600,  0.600,  0.600,  0.200)  ) ; 


/* 

BLACK 

*/ 

/* 

RED 

*/ 

/* 

GREEN 

*/ 

/* 

YELLOW 

*/ 

/* 

BLUE 

*/ 

/* 

MAGENTA 

*/ 

/* 

CYAN 

*/ 

/* 

WHITE 

*/ 

/* 

EARTH 

*/ 

/* 

COAST 

V 

/* 

BORDER 

V 

/* 

LATLONG 

*/ 

/* 

SUN 

*/ 

/* 

SAT_RED 

V 

/* 

SAT_GREEN 

*/ 

/* 

SAT_YELLOW 

*/ 

/* 

SAT  BLUE 

V 

/* 

SAT  MAGENTA 

*/ 

/* 

SAT  CYAN 

*/ 

/* 

SAT_WHITE 

V 

/* 

SAT  GRAY 

*/ 
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j ★  A*******************************  ★  / 

/*  *  / 

/*  smet_addelete . c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

f*  This  file  contains  the  routines  used  by  SMET  to  */ 

/*  add  and  delete  satellites,  add  and  delete  plots  and  */ 

/*  switch  and  crash  satellites.  */ 

/*  V 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/**********************************/ 

finclude  "smet.h" 

/*  MAKE_NEW_SATELLITE  */ 

/*  Creates  a  new  satellite  and  sets  all  default  settings  */ 

/*  Asks  user  for  name  if  none  is  provided.  */ 

makenewcurrentsat (name) 
char  ‘name; 

( 

satellite  "sat; 

sat  —  (satellite  *) malloc (sizeof (satellite) ) ; 
sat->epoch  -  clock; 

if  (strlen (name) )  strcpy (sat ->name, name) ; 

else  dialog ("Enter  name  of  new  satellite.",  sat->name) ; 

sat->color  «  SAT_P.ED  +  rand()  %  6; 

sat -> fill  -  FALSE; 

sat->vectors  -  FALSE; 

sat->plot  «  FALSE; 

sat->plots  *  NULL; 

sat->maneuver  *  NULL; 

sat->next  current_sat; 

header_sat  next  —  current_sat  —  sat; 

user_new_sat () ; 

» 
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/*  ADD_SATELLITE  */ 

/*  Set  the  five  basic  parameters  of  the  orbit  and  init-  */ 

/*  ialize  the  dials  as  appropriate.  */ 

add_satellite (n) 

int  n;  /*  orbit  type  selected  by  pop  up  menu  */ 

< 

make_new_current_sat ("") ; 
switch  (n) { 


case  1:  /*  300  nm  circular  LEO  */ 

current_sat->ma jor  —  1  +  300/distance_factor [ENGLISH] ; 
current_sat->ecc  “  0.0; 
current_sat->incl  -  DEG30; 
current_sat->asc  ”  DEG30; 
current_sat->peri  ■  0.0; 
break; 

case  2:  /*  Molniya  Orbit  */ 

current_sat->ma jor  —  pow (ROTATION_RATE*TU*2,  -2. 0/3.0); 
current_sat->ecc  ■  0.75; 
cur rent_sat->incl  —  asin (sqrt (0 . 8) ) ; 
current_sat->asc  -  DEG60; 
cuirent_sat->peri  —  DEG270; 
break; 

case  3:  /*  Geosynchronous  Orbit  */ 

current_5at->ma jor  ■«  pow (ROTATION_RATE*TU,  -2. 0/3.0); 
current_sat->ecc  -  0.0; 
current_sat->incl  ”  0.0; 
current_sat->asc  «  0.0; 
current_sat->peri  -  0.0; 
break ; 

case  4:  /*  Custom  Orbit  */ 

ask_ma  jor ( ) ; 
ask_ecc ( ) ; 
ask_incl ( )  ; 
ask  asc ( )  ; 
ask_peri ()  ; 
break ; 

) 

update_static_parameters (current_sat ) ; 

update_sat_matrix (current_sat) ; 

reset  dials ();  /*  set  dials  to  current  satellite 


V 
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/*  DELETE_SATELLITE 

/*  Removes  the  current  satellite  from  SMET 

delete_satellite () 

( 

header_sat .next  -  current_sat->next ; 
delete_plots (current_sat) ; 
free((char  *)  current_sat) ; 
current_sat  -  header_sat.next; 
if  (current  sat  !<■  NULL)  reset_dials  () ; 


*/ 

*/ 


/*  ADD_PLOT  &  DELETE_PLOTS  */ 

/*  Add  a  plot  to  the  linked  list  or  purge  the  list.  */ 

addplot (sat) 
satellite  *sat; 

I 

plot  struct  *new_plot;  /*  link  for  new  lat/lon  */ 

new_plot  -  (plot_struct  *)malloc (sizeof (plot_struct) ) ; 

new_plot->lat  ”  sat->lat; 

new_plot->lon  —  sat->lon; 

rew_plot->next  **  sat->plots; 

sat->plots  -  new_plot; 


deleteplots  (sat) 
satellite  *sat; 

( 

plot_struct  *temp; 

for  (temp  -  sat->plots;  sat->plots  !«  NULL;  temp  -  sat->plots) 

( 

sat->plots  -  sat->plots->next ; 
free ((char  *)  temp); 

) 

) 
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/*  SWITCH_CURRENT_SAT  */ 

/*  Figures  out  which  satellite  was  chosen  from  mouse  y  */ 

/ *  position.  Walk  down  list  until  you  got  it  and  switch' em.  */ 

switch_current_sat (y) 
short  y; 

( 

int  sat_number  —  (float) (satlistwin . sizey  -  y)  /  (float) LINE  -  3, 
count;  / *  dummy  loop  counter  */ 

satellite  ‘walker; 

if  (sat_number  >  1) 

( 

walker  —  current_sat; 

for  (count  <■  1;  walker  !—  NULL  ti  count  +  1  <  s_~_number;  count++) 
walker  «•  walker->next ; 
set_next_to_current (walker) ; 

) 

) 


/*  CRASHSATELLITE  ‘/ 

/*  Set  the  desired  sat  to  the  current  position  so  that  */ 

/*  you  can  do  a  delete  satellite.  Called  when  range  >—  1 .  */ 

crash_satellite (sat,  previous_sat) 
satellite  *sat,  *previous_sat; 

( 

char  st r [ 50 ) ; 

sprintf(str,  "%s  crashed  into  the  earth  and  was  deleted. ", sat->name) ; 
message (str )  ; 

set_next_to_current (previous_sat)  ; 
delete_satellite ()  ; 
sat  -  previous_sat ; 


/*  SET_NEXT_TO_CURRENT  */ 

/*  Used  by  two  previous  procedures  to  switch  the  linked  */ 

/*  list  of  satellites.  */ 

set_next_to_current (previous) 
satellite  ‘previous; 

( 

if  (previous  !-  NULL  &&  previous->next  !-  NULL  &&  previous  !»  4header_sat) 
( 

header_sat . next  -  previous->next ; 
previous->next  «  previous->next->next ; 
header_sat . next- ’next  -  current_sat; 
current_sat  -  header_sat . next ; 
reset_dials ( ) ; 

) 

) 


68 


mnmt  uut.c  Mon  Sap  3  17:58:57  1990  1 

/A****************************'*****/ 
/*  */ 
/*  smet_axes.c  */ 
/*  by  Carlos  1.  Noriega  */ 
/*  V 
/*  This  file  contains  the  routines  to  draw  the  axes  */ 
/*  used  by  the  SMET  demonstrator.  */ 
/*  */ 
/*  Written  -  September  3,  1990  */ 
/*  Modified  -  September  3,  1990  */ 
/*  */ 
/**********************************/ 


♦include  "smet.h" 


/*  DRAW-AXES  */ 
/*  draw  the  IJK  axes  at  the  center  of  the  */ 
/ *  world  coordinate  system.  Remember  that  */ 
/*  IJK  is  not  aligned  the  same  as  the  systems  */ 
/*  XY2  coordinate  system.  */ 


draw_axos() 

( 

/*  qo  to  white  */ 
c3f  (Color (WHITE) )  ; 

/*  set  the  linewidth  -  what  a  comment  ! !  */ 

linewidth (3) ; 

/*  label  the  I  axis  */ 
cmov (0.0,  0.2,  1.2); 

fmpistr ("I") ? 

/*  label  the  J  axis  */ 
cmov (1.2,  0.2,  0.0); 
fmprstr ( " J" ) ; 

/*  label  the  K  axis  */ 
cmov (0.2,  1.2,  0.0); 

fmprstr ( "K") ; 

/*  draw  the  J  axis  */ 
move  (1.0,  0.0,  0.0); 
draw  (12. 5,  0.0,  0.0); 

/*  draw  the  K  axis  */ 
move  (0.0,  1.0,  0.0); 

draw(0.0,  12.5,  0.0); 

/*  draw  the  I  axis  */ 
move  (0.0,  0.0,  1.0); 
draw(0.0,  0.0,  12.5); 
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/*****«****************************/ 
/*  */ 

/*  smet_dialog. c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  V 

/*  This  file  contains  the  routines  to  give  the  user  messages  */ 

/*  and  ask  him  questions.  It  is  used  by  the  SMET  */ 

/*  demonstrator.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  ‘/ 

/********************************* 

♦include  "smet.h" 

/*  DIALOG  &  MESSAGE  */ 

/*  Use  the  dialogbox  routines  to  send  a  message  or  ask  a  ‘/ 

/*  question.  */ 


dialog (prompt, response) 
char  ‘prompt,  ‘response; 

( 

opendialogbox () ; 

get_user_response (prompt,  response)  ; 
close_dialogbox ( ) ; 

I 


message (prompt) 
char  ‘prompt; 

( 

opendialogbox () ; 
clearprompt (prompt) ; 
ringbell () ; 
sleep  (3) ; 

close_dialogbox ( ) ; 
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/*  DIALOGBOX  Routines  */ 

/*  Uses  the  hidden  dialogwin  for  all  user  messages.  */ 

open_dialogbox()  /*  "POP"  the  box.  */ 

( 

winset (dialogwin . gid) ; 

ortho2 (0 . 0, dialogwin. sizex,  0.0,  dialogwin. sizey) ; 

loadmatrix (ident_matrix) ; 

singlebuffer () ; 

gconfig () ; 

winpopO  ; 

zciear  ()  ; 

qdevice (KEYBD) ; 


clearprompt (prompt)  /*  Clear  the  box  and  type  the  message.  */ 

char  "prompt; 

( 

c3f (Color [CYAN] ) ; 
clear ( ) ; 

c3f (Color [BLACK] ) ; 
linewidth (2)  ; 

recti(7,  1,  dialogwin . sizex-7,  dialogwin.sizey-7); 
cmov2i(12,  20); 

charstr  (prompt )  ;  charstrC  "); 

» 


close_dialogbox ( )  /*  "PUSH"  the  box.  */ 

I 

zciear  ()  ; 
winpush  ( ) ; 
unqdevice (KEYBD) ; 


71 


smt  dialog. c 


Mon  Sap  3  17:56:57  1990 


3 


/*  GET_USER_RESPONSE  */ 

/*  a  simple  parser  for  user  interaction.  Does  not  currently  */ 

/*  check  if  response  string  is  shorter  than  fixed  length  of  40.*/ 

/*  Ignores  events  other  than  the  keyboard.  */ 

get_user_response (prompt, response) 
char  *prompt,  * response ; 

( 

int  cur_str_len, 

donetyping  -  FALSE; 
short  c; 
char  *str; 

Device  dev; 

response fO]  —  '\0'; 
cur_str_len  -  strlen (response) ; 

clearprompt (prompt) ;  /*  display  prompt  */ 

while (! donetyping)  /*  read  till  eof  ('\n'  or  '\r')  */ 

if (dev  —  qread(ic)) 

{ 

if (dev  !-  KEYBD)  continue;  /*  don't  care  */ 

switch (c) 

( 

case  '\027' :  /*  AW  sets  cursor  back  to  start  */ 

responsefO]  »•  '  \0' ; 
cur_str_len  ■  strlen (response) ; 
clearprompt (prompt) ; 
break ; 

case  ' \n' : 

case  '\r':  /*  carriage  return  or  enter  */ 

donetyping  —  TRUE; 
break ; 

case  127 : 

case  ' \b' :  /*  backspace  or  delete  */ 

if (cur_str_len)  ( 

response ( — cur_str_len]  —  ' \0' ; 
clearprompt (prompt) ; 
charstr (response) ;  ) 

else  ringbellO; 
break ; 

default:  /*  add  to  string  */ 

if (cur_str_len  <  41)  ( 

str  -  iresponse (cur_ str_len] ; 
response [cur_str_len++]  —  c; 
response (cur_str_len]  -  '\0'; 
charstr  ( str) ;  ) 

else  ringbellO; 
break ; 

) 

I 

I 
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/*  ********************************  * j 
!*  */ 
!*  smet_files.c  */ 
/*  by  Carlos  I.  Noriega  */ 
/*  */ 
/*  This  file  contains  the  routines  used  by  SMET  to  */ 
/*  read  and  write  files.  */ 
l*  */ 
/*  Written  -  September  3,  1990  */ 
/*  Modified  -  September  3,  1990  */ 
/*  */ 
/**********************************/ 


♦include  "smet.h" 


/ *  READ_WORLD_F I LE 

*/ 

/*  Reads  earth  data  for  map  and  globe  from  file  in  the 

V 

/*  current  directory. 

*/ 

read  world  file() 

TILE  *fd; 

fd  *  fopen ("smet  data . d", "r") ; 

if ( f d  MULL) 

( 

message ( "Unable  to  open  file  smet_data.d,  quiting  SMET"); 
program_exit ( ) ; 

) 

else 

( 

lseek (fd, 0L, 0)  ; 

f read (if ine3d_data,  sizeof (fine3d  data) ,  1,  fd)  ; 
f read ( 4 f ine2d_dat a,  sizeof  ( f ine2d_data) ,  1,  fd) ; 
f read (4coarse3d_data,  sizeof (coarse3d_data)  ,  1,  fd); 
f read ( icoa r se2d_data,  sizeof (coarse2d_data) ,  1,  fd) ; 
f read (iborder3d_data,  sizeof (border3d_data) ,  1,  fd) ; 
f read ( 4border2d_data,  sizeof (border2d_data) ,  1,  fd) ; 
f close (fd) ; 

) 

) 
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/*  ASK_AND_READ_SAT_FILE  &  RE AD_S AT_F I LE  */ 

/*  Reads  satellite  data  previously  saved.  */ 

ask_and_read_sat_f ile () 

( 

char  name[40J; 

dialog ("Enter  the  name  of  the  satellite  f ile name) ; 
read  sat  file (name); 


read_sat_f ile (f ile_name) 
char  f ilename ( ] ; 

( 

FILE  *fd; 
satellite  *sat; 
int  magic; 
char  str [ 50  J ; 

fd  *  fopen (file_name, "r") ; 
if (fd  -«  NULL) 

( 

sprintf (str, "Unable  to  open  file  -  %s . " , f ile_name) ; 
message (str) ; 

) 

else 

( 

lseek ( fd, 0L, 0) ; 

fread (&magic,  sizeof(int),  1,  fd)  ; 

if  (magic  !-  4567)  /*  ensure  properly  formatted  file  V 

( 

sprintf (str, "%s  is  not  a  satellite  f ile . ",  f ile_name) ; 
message (str) ; 

) 

else 

(  /*  read  satel.1  ites  until  you  get  a  read  error  */ 

sat  =  (satellite  *)malloc (sizeof (satellite)  )  ; 
while  (fread(sat,  sizeof  (satellite) ,  1,  fd) ) 

I 

sat->plots  —  NULL; 
sat->maneuver  —  NULL; 
sat->next  ”  current_sat; 
current_sat  -  sat; 

sat  -  (satellite  *) malloc (sizeof (satellite)  )  ; 

) 

f ree ( (char  *)  sat); 

header_sat . next  «*  current_sat;  /*  make  last  one  current  */ 
reset_dials ( ) ; 

I 

f close (fd) ; 

) 

) 


74 


omet_fil«9 . c  Mon  S«p  3  17:58:58  1990  3 

/*  ASK_AND_WRITE_SAT_FILE  &  WRITE_SAT_FILE  */ 

/*  Write  satellite  data  into  a  file.  */ 

ask_and_write_sat_f ile () 

( 

char  name [40],  number [5]; 

dialog ("Enter  the  name  of  the  satellite  file name) ; 
dialog ("Enter  the  number  of  the  satellites  saved. ",  number) ; 
write  sat  file (name, atoi (number) ) ; 


write_sat_f ile ( f ile_name, numsats) 
char  f ilename ( ] ; 
int  numsats; 

( 

FILE  *fd; 
satellite  *sat; 
int  c,  magic  =  4567; 
char  str l 50 ] ; 

fd  =  f open (f ile_name, "w" ) ; 
if (fd  ==  NULL) 

( 

sprint f  (str, "Unable  to  open  file  -  %s . ",  f ile_name) ; 
message  (str) ; 

) 

else 

( 

lseek (fd, 0L, 0) ; 

fwrite (Smagic,  sizeof (int) ,  1,  f d)  ;  /*  ensure  properly  formatted  file  */ 
sat  -  current_sat; 

for  (c  »•  0;  (c  <  numsats  &&  sat  !—  NULL);  C++) 

< 

fwrite  (sat,  sizeof  (satellite) ,  1,  fd) ; 
sat  “  sat->next; 

) 

f close (fd) ; 

) 

) 
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/*  TAKE-PICTURE  */ 
/*  Use  image  library  routines  to  write  an  image  */ 
/*  in  the  selected  format.  Black  and  White  are  */ 
/*  inverte--’  *ir  use  in  EPS  format.  */ 


take_picture (num) 
int  num; 

( 

NPS IMAGE  *imgl,*img2,  *img3;  /*  temp  ptr  to  images  */ 
char  filename [50] , 
str (75] ; 

switch  (num) 

( 

case  1: 

dialog ("Enter  the  name  for  the  RGB  file  filename); 

message ("Saving  the  screen  will  take  about  90  seconds."); 
ref resh_main ( ) ; 

save_window_as_sgi_rgbimage (filename,  backwin. gid) ; 
break; 
case  2: 

dialog ("Enter  the  name  for  the  RGB  file  filename); 

message ( "Saving  the  screen  will  take  about  90  seconds."); 
invert_bw ( ) ; 

save_window_as_sgi_rgbimage (filename,  backwin . gid)  ; 
revert_bw ( ) ; 
break; 
case  3 : 

dialog ("Enter  the  name  for  the  EPS  file  -",  filename); 
message ( "Saving  the  screen  will  take  about  60  seconds."); 

invert_bw () ; 

imgl  =  read_window_fb_as_rgb_npsimage (backwin. gid, filename) ; 
img2  *  shrink  npsimage (imgl,  2)  ; 
delete_an_image (imgl) ; 
img3  -  rotate_npsimage (img2) ; 
delete_an_image (img2) ; 
write_eps_image (filename,  img3,  125); 
delete_an_image (img3)  ; 
revert_bw ( ) ; 
break ; 
case  4: 

dialog ("Enter  the  name  for  the  EPS  file  filename); 
message ( "Saving  the  main  window  will  take  about  30  seconds.") 
invert_bw ( ) ; 

imgl  -  read_window_fb_as_rgb_npsimage (mainwin. gid, filename) ; 

img2  -  shrink_npsimage (imgl,  2); 

delete_an_image (imgl)  ; 

write_eps_image (filename,  img2,  125); 

delete_an_image (img2) ; 

revert_bw ( ) ; 

break ; 

) 

I 
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*/ 
*/ 


/ *  Support  routines  for  TAKE_PICTURE 
/* 


refresh_main  ') 

{ 

switch (window_order)  /*  refresh  main  window  after  message  */ 

{ 

case  OMNISATMERC: 

case  OMHIMERCSAT:  draw_omni_view ( )  ;  break; 
case  SATOMNIMERC: 

case  SATMERCOMNI :  draw_sat_view () ;  break; 
case  MERCOMNISAT: 

case  MERCSATOMNI:  draw_merc_view ( ) ;  break; 

} 

> 


invert_bw ( ) 

( 

Color [BLACK] [0]  -  Color [ BLACK] [ 1 ]  -  Color [BLACK] [2]  -  Color [BLACK] [3]  -  1.0; 
Color [WHITE] [0]  -  Color [WHITE] [ 1 }  -  Color [WHITE] [2]  -  Color [WHITE] [3]  - 
Color [SATWHITE] [0]  -  Color [SAT_WHITE ][ 1 ]  -  Color [SAT_WHITE] [2]  - 

Color [SAT_WHITE] [3]  -  0.25 

winset (backwin . gid)  ; 
c3f  (Color [BLACK]) ; 
clear ()  ; 
swapbuf fers ( )  ; 
draw_windows ()  ; 

draw  windows));  /*  draw  twice  to  ensure  that  proper  buffers  are  showing  */ 


revert_bw ( ) 

( 

Color [BLACK] [0]  -  Color [BLACK] [1]  -  Color [BLACK] [2]  -  Color [BLACK] [3 ]  -  0.0; 
Color [WHITE] [0]  -  Color [WHITE] [ 1 ]  -  Color [WHITE] [2]  -  Color [WHITE] [3]  - 
Color [SAT_WHITE] [0]  -  Color [SAT_WHITE ][ 1 ]  -  Color [ SAT_WHITE] [ 2 ]  - 

Color [SAT_WHITE] [3]  -  1.0; 

winset (backwin . gid)  ; 
c3f (Color [BLACK] ) ; 
clear  ( ) ; 
swapbuf fers () ; 
draw  windows () ; 
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/**********************************/ 


/*  */ 

/*  smet_globe.c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/ *  This  file  contains  the  routines  used  by  SMET  to  */ 

/*  draw  the  globe  in  a  viewing  window.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 


/*  **  +  ***■****************■*****■*****  */ 
♦include  "smet.h" 

/*  DRAW-GLOBE  */ 

/*  Assumes  that  a  central  disk  normal  to  viewer  drawn  */ 

/*  by  calling  routine.  */ 

draw_globe (rotation) 

float  rotation;  /*  amount  of  y-rotation  of  the  globe  */ 

( 

pushmatrix () ; 

rot (rotation*TO_DEO,  '  y' ) ;  /*  turn  the  globe  */ 

draw_coasts3d () ; 
draw_borders3d ()  ; 
draw_latlons3d ( )  ; 

popmatrix ()  ; 

I 


/*  DRAW-COASTS 3D  */ 

/*  Draw  the  desired  data  FINE  or  COURSE.  */ 

draw_coasts3d  () 

( 

int  c; 

linewidth  (2)  ; 
c3f (Color [COAST] ) ; 
if  (fineflag) 

for  (c  -  0;  c  <  FINEPOINTS;  C++) 
if  (fine3d_data.pen[cj) 

draw (f ine3d_data .x [c] ,  fine3d_data .y [c] ,  f ine3d_data . z [c] ) ; 
else 

move (fine3d_data.x[c] ,  fine3d_data.y[c] ,  fine3d_data . z [c] ) ; 

else 

for  (c  -  0;  c  <  COARSEPOINTS;  C++) 
if  (coarse3d_data.pen[c] ) 

draw (coarse3d^data.x[c] ,  coarse3d_data .y [c] ,  coarse3d_data . z [ c ] ) 
else 

move (coarse3d_data.x(c] ,  coarse3d_data .y [c] ,  coarse3d_data. z[c] ) 

) 
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*/ 

*/ 


I*  DRAW-BORDERS 3D 

/*  Draw  the  border  data  if  desired. 

draw_borders3d () 

{ 

int  c ; 

linewidth (1)  ; 
c3f (Color [BORDER] )  ; 
if  (borderflag) 

for  (c  -  0;  c  <  BORDERPOINTS;  C++) 
if  (border3d_data.pen[c] ) 

draw (border3d_data. x[c] ,  border 3d_dat a .  y [c] ,  border3d_data . z [c] ) 
else 

move (border3d  data.x[c],  border3d  data.y[c],  border3d  data.zfc]) 


/*  DRAW-LATLONS3D  */ 

l*  Draw  the  latitude  and  longitude  lines  on  the  globe.  */ 

draw_latlons3d () 

( 

int  c; 

linewidth  ( 1 )  ; 

c3f (Color [LATLONG] ) ; 

pushmatrix ( ) ; 

/*  draw  the  longitude  lines  */ 
pushmatrix ()  ; 

for  (c  -  1;  c  <*=  12;  C++)  /*  draw  an  arc  every  30  degrees  */ 

1 

rotate (300, ' y' ) ; 

arc (0.0, 0.0,  1.0, -900,  900)  ; 

) 

popmatrix ( )  ; 


/*  draw  the  latitude  lines  */ 

rotate (900, ' x' ) ;  /*  flip  over  on  its  side  */ 

circ  (0 . 0, 0 . 0, 1 . 0) ;  /*  equator  */ 

translate (0.0, 0.0, -0.5) ; 

circ (0 . 0, 0 . 0, 0 . 866) ;  /*  30  degrees  north  */ 

translate (0.0,0.0,-0.366) ; 

circ  (0 . 0, 0 . 0, 0 . 5) ;  /*  60  degrees  north  */ 

translate (0.0, 0.0,  1.366)  ; 

circ  (0 . 0, 0 . 0, 0 . 866) ;  /*  30  degrees  south  */ 

translate (0.0,0.0,0.366) ; 

circ (0 . 0, 0 . 0, 0 . 5) ;  /*  60  degrees  south  */ 


popmatrix ( ) ; 

) 
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/**********************************/ 
/*  */ 

/*  smet_axes.c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  initialization  routines  */ 

/*  used  by  the  SMET  demonstrator.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/**********************************/ 


tinclude  "smet.h" 


/*  1NITIALIZE_SMET  */ 
/*  Open  windows,  set  up  menus,  initialize  queue,  */ 
/*  draw  the  user  instructions.  */ 


initialize_smet (argc,argv) 
int  argc; 
char  **argv; 

{ 

int  c; 

fmfonthandle  paramfont; 


/*  initialize  windows 


init_window (Sbackwin)  ; 

/* 

a  black  background 

V 

init_window (fidialogwin) ; 

/* 

the  dialog  window 

*/ 

init  window (isatlistwin) ; 

/* 

where  the  satellite  names  are  listed 

*/ 

init_window (iparamwin)  ; 

/" 

where  the  parameters  are  listed 

*/ 

init_window (Smainwin) ; 

/" 

the  centered  large  window 

*/ 

init  window (iauxwinl) ; 

/* 

the  lower  left  small  window 

*/ 

init_window (Sauxwin2) ; 

/* 

the  lower  right  small  window 

*/ 

message (" 

Loading  SMET") ; 

/*  make  the  pop  up  menus 


colorjmenu  —  defpup ("Color  %t  %F  |  Red  %x0  |  Green  %xl  |  Yellow  %x2  \ 

|  Blue  %x3  |  Magenta  %x4  |  Cyan  %x5  |  White  %x6  |  Gray  %x7",  smet_sat_color) ; 


orbit_menu  —  defpup ("Define  Orbit  %t  %F  |  Low  Earth  |  Molniya  \ 
|  Geosynchronous  I  Input  Parameters",  add_satellite) ; 


add_menu  ■  defpup ("New  Satellite  %t  I  Define  Orbit  %m  |  Define  Launch  %f  \ 

I  Read  File  %f",  orbit_menu,  add_launch,  ask_and_read_sat_f ile) ; 

param_menu  -  defpup ("Parameter  Options  %t  %F  |  Classical  |  Secondary  \ 

I  Position  |  Options  |  Combined",  smet _param_option) ; 

picture_menu  ■>  defpup ("Picture  Options  %t  %F  |  Color  RGB  Image  (Regular)  %xl  \ 

|  Color  RGB  Image  (Inverted)  %x2  |  Post  Script  -  Entire  Screen  %x3  I  Post  Script  \ 

-  Main  Window  %x4",  take_picture) ; 


top_menu  ■  defpup ("SMET  menu  %t  |  Add  Satellite  %m  |  Maneuver  Satellite  %f  \ 
|  Delete  Satellite  %f  |  Save  Satellites  %f  |  Take  Picture  %m  |  Satellite  \ 
Parameters  %m  I  Instructions  %f  I  Exit  %f  ",  add_menu,  add_maneuver, 
delete_satellite,  ask_and_write_sat_file,  picture_menu,  param_menu, 
draw  instructions,  program  exit); 
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/*  queue  the  necessary  devices  */ 

qdevice (MIDDLEMOUSE) ; 
qdevice  (RIGHTMOUSE)  ; 
qdevice (D1AL2) ; 
qdevice (DIAL3) ; 
qdevice (DIAL4) ; 
qdevice (DIAL5) ; 
qdevice (DIAL6) ; 
qdevice (DIAL7)  ; 

tie (MIDDLEMOUSE, MOUSEX, MOUSEY) ; 

qreset  () ;  /*  initialize  the  queue  */ 

/*  dial  limits:  one  minute  to  6  hours  */ 
setvaluator (DIAL3,  (short) (clock_speed) ,  1,  360); 


/*  miscellaneous  initializations 


*/ 


/*  read  data  files  for  drawing  globe  and  map  */ 

read_world_f ile () ; 

/*  Initialize  the  satellite  list  */ 

headersat .next  «  NULL; 

/*  read  any  sat  files  listed  by  user  on  startup  */ 

for  (c  «  1;  c  <  argc;  C++)  read_sat_file (argvfcj ) ; 

/*  Initialize  font  for  windows  */ 

fminit  ()  ; 

fmsetfont (fmscalefont (fmfindfont ("Times-Roman") ,  LETTER) ) ; 

/*  User  Initialition  Routine  */ 

user  init ( ) ; 


/*  put  the  different  views  into  the  default  windows  */ 

omni_view_gid  —  mainwin.gid; 

sat_view_gid  —  auzwinl . gid; 

merc_view_gid  —  au::win2  .  gid; 

window  order  -  OMNISATMERC; 


/*  draw  instructions  for  the  user  */ 

draw  instructions (); 
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/*  INIT-WINDOW  */ 
/*  Set  up  a  window  using  the  data  from  the  given  */ 
/*  window  structure  */ 


initwindow (win) 

window_3truct  *win;  /*  window  structure  */ 

( 

/*  set  the  pref  position  according  to  structure  */ 

prefposition (win->orgx,  win->orgx  +  win->sizex, 
win->orgy,  win->orgy  +  win->si*ey) ; 


win->gid  —  winopen(""); 

/* 

create  a  window  */ 

noborder  () ; 

/* 

no  edges  or  window  manager  pop  ups 

*/ 

RGBmode ( ) ; 

/* 

RGB  color  specifications 

V 

doublebuffer () ; 

/* 

no  flickering 

*/ 

zbuffer (TRUE) ; 

/* 

hidden  surface  removal 

*/ 

lsetdepth(0,  0x7fffff) ; 

/* 

get  the  largest  Zbuffer  resolution 

*/ 

gconfig () ; 

/* 

attach  settings  to  the  current  window  * / 

mode  (MVIEWING)  ; 

/* 

put  winod  in  viewing  mode 

*/ 

/*  go  black  and  clear  both  buffers  */ 

c3f (Color [BLACK] ) ; 

clear  ()  ; 

swapbuffers  ()  ; 

clear  ()  ; 
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I *  ********************************  */ 
/*  */ 
/  *  smet_)cepler .  c  */ 
/*  by  Carlos  I.  Noriega  */ 
/*  */ 
/*  This  file  contains  the  routines  to  update  orbital  */ 
/*  parameters.  It  is  used  by  the  3MET  demonstrator.  */ 
/*  */ 
/*  Written  -  September  3,  1990  */ 
/*  Modified  -  September  3,  1990  */ 
/*  */ 
/**********************************/ 


♦include  "smet.h" 


/*  UPDATE-STATIC-PARAMETERS 
/*  See  Chapter  II  of  thesis. 

update_static_parameters (sat) 
satellite  *sat; 

( 

/*  update  the  semi-minor  axis,  period,  perigee,  apogee  when  */ 

/*  ecc  or  major  have  changed  due  to  a  maneuver  or  user  input  */ 
sat->minor  —  sat->ma jor  *  sqrt (1  -  sat->ecc  *  sat->ecc) ; 
sat->period  —  pow (sat->ma jor, 1 . 5)  *  DEG360  *  TU; 
sat->apogee_range  -  sat->ma jor  *  (1  +  sat->ecc) ; 

sat->perigee_range  -  sat->major  *  (1  -  sat->ecc) ; 
update_perturb_constants (sat) ; 

) 


*/ 

*/ 


/*  UPDATE-PEP.TURB_CONSTANTS  */ 

/*  See  Chapter  II  of  thesis.  */ 

update_perturb_constants (sat) 
satellite  *sat; 

( 

float  perturb_constant  —  J2  /  (pow (sat->ma jor,  3.5)  *  pow((l  -  sat->ecc  *  sat->ecc) ,  2)) 
/*  constant  for  regression  of  line  of  nodes  */ 

sat->asc_change  -  -  1.5  *  cos (sat->incl)  *  perturb_constant  /TU; 

/*  constant  for  rotation  of  line  of  apsides  */ 

sat->peri_change  -  (3  -  3.75  *  sin (sat->incl) *sin (sat->incl) )  *  perturb  constant  /TU; 
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/*  UPDATE-CHANGING-PARAMETERS 
/*  See  Chapter  II  of  thesis. 

update_changing_parameters  {sat,  time) 
satellite  *sat; 
float  time; 

( 

float  E,  M;  /*  Eccentric  and  Mean  anomalies.  */ 

if  (clock_running  &&  perturbations) 

( 

/*  regression  of  line  of  nodes  */ 

sat->asc  *  sat->asc  +  sat->asc_change  *  clock_speed; 

/*  rotation  of  line  of  apsides  */ 

sat->peri  —  sat->peri  +  sat->peri_change  *  clock_speed; 
update_sat_matri::  (sat)  ? 

) 

/*  use  M  as  initial  approximation  of  eccentric  anomaly  */ 

sat->E  =  M  *  fmod(((time  -  sat->epoch)  *  DEG360) /sat->period,  DEG360) ; 

if  (M  <  0.0)  sat->E  -  M  -  DEG360  +  M;  /*  ie .  time  <  epoch  */ 

/*  compute  the  Eccentric  anomaly  using  Newtonian  Approximation  */ 

/*  normally  converges  after  3  to  5  iterations  at  this  accuracy  */ 
do 
( 

E  «  sat->E; 

sat->E  »  E  +  (M  -  E  +  sat->ecc  *  sin(E))/(l  -  sat->ecc  *  cos(E)); 

)  while  (nonzerop (sat->E  -  E) ) ;  /*  repeat  till  almost  zero  */ 

/*  compute  the  true  anomaly  from  the  eccentric  anomaly  */ 

sat->anom  -  acos ( (cos (sat->E)  -  sat->ecc) / (1  -  sat->ecc  *  cos ;sat->E) ) ) ; 

if  (sat->E  >  DEG180)  sat->anom  -  DEG360  -  sat->anom;  /*  get  correct  quadrant  */ 

/*  compute  the  distance  to  the  satellite  */ 

sat->range  ■  sat->major  *  (1  -  sat->ecc  *  cos (sat->E) ) ; 

/*  compute  the  velocity  of  the  satellite  */ 

sat->velocity  «  sqrt ( (2/sat->range)  -  (l/sat->ma jor) ) ; 

/*  compute  the  planar  coordinates  of  the  satellite  */ 

sat->p  —  sat->range  *  cos (sat->anom) ; 
sat->q  ■  sat->range  *  sin (sat->anom) ; 

/*  compute  the  XYZ  coordinates  of  the  satellite.  Don't  */ 

/*  forget  that  this  is  for  an  IJK  system  which  is  */ 

/*  different  than  the  display  systems  XYZ  system  */ 

sat->y  -  sat->q*sat->matrix[0)  (0 ) +sat->p*sat->matri:: [2 ]  [ 0 J +sat->matrix [ 3 ]  [0J; 
sat->z  —  sat->q*sat->matrix(0] ( 1 ] +sat->p*sat->matrix [2] (1 J+sat->matrix(3) (1J ; 
sat->x  -  sat->q*sat->matrix(0] [2]+sat->p*sat->matrix[2] [2] +sat->matrix [3] [2] ; 

/*  compute  the  sat's  lat  and  long  for  plotting  adjustinf  for  earth  spin  */ 
sat->lat  ”  asin(sat->z  /  sat->range) ; 

sat->lon  —  atan2 (sat->y,  sat->x)  -  fmod(time  *  ROTATION_RATE, DEG360) ; 
if  (sat->lon  <  -DEG180)  sat->lon  •«  sat->lon  +  DEG360; 

/*  save  the  plot  if  it  is  being  plotted  */ 
if  (sat->plot)  add_plot (sat) ; 


*/ 

*/ 


) 
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/*  UPDATE-SAT_MATRIX  */ 
/*  compute  the  orbital  plane  transformation  matrix  */ 
/*  lets  let  the  graphics  hardware  do  some  of  our  */ 
/*  tough  calculations.  what  a  neat  trick!  */ 


update_sat_matrix (sat) 
satellite  *sat; 

( 


loadmatrix (ident_matrix) ; 
rot (sat->asc*TO_DEG, '  y' ) ; 

/* 

longitude  of  the  ascending 

node 

*/ 

rot (sat->incl*TO_DEG, ' z' ) ; 

/* 

inclination 

of  the  orbital 

plane 

*/ 

rot (sat->peri*TO_DEG,  '  y' )  ; 
getmatrix (sat->matrix) ; 

/* 

argument  of 

periapsis 

*/ 

) 


/*  ALT_UFDATE-SAT_MATRIX  */ 
/*  Use  for  multiprocessing  when  encorporated.  */ 
/ *  Gets  same  matrix  as  update_sat_matrix .  */ 
/*  Slightly  different  than  matrix  in  Chapter  2  of  */ 
/*  thesis  due  to  difference  in  coordinate  systems.  */ 


alt_update_sat_matri::  (sat) 
satellite  *sat; 

{ 

float  cos^asc  =  cos (sat->asc) , 
sin^asc  «■  sin (sat->asc)  , 
cos_incl  -  cos (sat->incl)  , 
sin  incl  =*  sin (sat->incl) , 
cos__peri  =  cos  (sat->peri)  , 
sin__peri  **  sin (sat->peri) ; 


sat->matrix [0] [0] 
sat->matrix[0] [1] 
sat->matrix ( 0 ] (2] 
sat->matrix [ 1 } (0) 
sat->matrix[l] [1) 
sat->matrix ( 1 1 [2] 
sat->matrix[2] [0] 
sat->matrix(2J [1J 
sat->matrix[2] [2] 
sat->matrix[0] [3] 
sat->matrix[3] [0] 
sat->matrix[3] [3] 


-sin_asc*sin_peri  +  cos_asc*cos_peri*cos_incl ; 
cos_peri*sin__incl; 

-cos_asc*sin_peri  -  sin_asc*cos_peri*cos_incl; 

-cos_asc*sin__incl; 

cos_incl; 

sin_asc*sin_incl; 

sin_asc*cos_peri  +  cos_asc*sin_peri*cos_incl; 
sin_peri*sin_incl ; 

cos_asc*cos_peri  -  sin_asc*sin_peri*cos_incl; 
sat->matrix(l]  [3]  -  sat->matrix[2]  [3]  ••  0.0; 
sat->matrix(3] [1]  -  sat->matrix[3] [2]  ■*  0.0; 
1.0; 
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/***************"***"**"*"**»******/ 
/*  */ 

/*  smet_launch. c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  routines  used  by  SMET  to  "/ 

/*  launch  a  satellite.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/**********************************/ 

♦include  "smet.h" 

/*  APD_LAUNCH  */ 

/*  Use  a  series  of  dialogs  to  get  required  info.  */ 

/*  See  Chapter  3  for  more  information.  */ 


add_launch ( ) 

( 

launch_struct  "newlaunch  -  (launch_struct  *)malloc (sizeof (launch_struct) ) 
char  str [ 15 ] ; 

dialog ( "Enter  name  of  new  satellite.",  new_launch->name) ; 

new_launch->time  -  ask_time ("time  at  burnout"); 

dialog  ("Enter  the  altitude  at  burnout",  str)  ; 

new_launch->alt  —  atof (str) /distance_f actor [unit_type] ; 

dialog < "Enter  the  velocity  at  burnout.",  str); 

new_launch->v  -  atof (str) /speed_f actor [unit^type] ; 

dialog ("Enter  the  latitude  at  burnout.",  str); 

new_launch->lat  *=  atof  (str)  *TO_RAD; 

dialog ("Enter  the  longitude  at  burnout.",  str) ; 

new_launch->lon  -  atof (str) *TO_RAD; 

dialog ("Enter  the  angle  above  the  horizon.",  str); 

new_launch->elev  -  atof (str) *TO_RAD; 

dialog ( "Enter  the  heading  angle.",  str); 

new_launch->azim  -  atof (str ) *TO_RAD; 

new_launch->next  -  pending_launches; 

pending_launches  —  new_launch; 


/*  LAUNCH_PENDING_SATS  */ 

/*  Check  all  pending  to  see  if  it  is  time  to  launch.  */ 

launch_pending_sats () 

I 

launch_struct  "launch,  "prev; 

for  (launch  -  pending_launches;  launch  !-  NULL;  launch  -  launch->next) 

( 

if  (launch->time  <-  clock) 

( 

initialize_orbital_parameters (launch)  ; 

if  (launch  -  pending_launches)  pending_launches  -  launch->next; 
else  prev->next  —  launch->next; 
free ( (char  *)  launch) ; 

) 

prev  —  launch; 

) 
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/*  INITIALIZE_ORBITAL_PARAMETERS  */ 

/*  See  Chapter  3  for  more  information.  */ 

initialize_orbital_parameters (launch) 
launch_struct  "launch; 

( 

char  str [70]; 
float  p, 

v_radial  —  launch->v  *  sin (launch->elev) , 

v_ns  —  launch->v  *  cos (launch->elev)  *  cos (launch- >azim) , 
vew  —  launch->v  *  cos (launch->elev)  *  sin (launch->azim) 

+  ROTATIONJRATE  *  TU  *  cos (launch— >!at) , 

v_hor  “  sqrt (v_ns*v_ns  +  v_ew*v_ew) ; 
satellite  "sat; 

make_new_current_sat (launch->name) ; 
sat  —  current_sat; 

/*  adjust  velocity,  elevation  angle  and  azimuth  for  earth  spin  "/ 
launch->v  —  sqrt ( v_radial*v_radial  +  v_hor*v__hor) ; 
launch->elev  -  atan2 (v_radial, v_hor) ; 
launch->azim  *  atan2  (v__ew,v_ns)  ; 

/*  compute  orbital  parameters  as  discussed  in  Chapter  3  of  thesis  */ 
sat->incl  —  acos (sin (launch->azim)  *  cos (launch->lat) ) ; 
sat->range  -  1  +  launch->alt; 

sat->ma jor  -  sat->range/ (2  -  sat->range  *  launch->v  *  launch->v) ; 
p  -  sat->range  *  v_hor  *  sat->range  *  v_hor; 
sat->ecc  -  f sqrt (1.0  -  p/sat->ma jor) ; 

sat->anom  *  atan2 (p  *  tan (launch->elev) ,  p  -  sat->range)  ; 

sat->E  -  acos ( (sat->ecc  +  cos (sat->anom) ) / (1  +  sat->ecc  *  cos (sat->anom) ) ) ; 
if  (sat->anom  >  DEG180)  sat->E  -  DEG360  -  sat->E;  /*  get  correct  quadrant  */ 

sat->epoch  -  launch->time  -  pow (sat->ma jor, 1 . 5)  "  (sac->E  -  sat->ecc*sin (sat->E) )  *  TU 
sat->peri  —  fmod (atan2 (sin (launch->lat) , cos (launch->lat)  *  cos (launch->azim) ) 

-  sat->anom  +  DEG360,  DEG360) ; 

sat->asc  —  fmod (launch->lon  +  clock*ROTATION_RATE  -  atan2 ( sin (launch->lat ) 

*sin (launch->azim) ,  cos  (launch->azim) )  +  DEG360,  DEG360); 

update_3tatic_parameters (sat) ; 
update_sat_matrix (sat) ; 
reset_dials ()  ; 

if  (0.0  >  sat->ecc  II  sat->ecc  >—  1.0) 

( 

sprintf (str, "%s  was  launched  into  an  open  orbit  and  deleted  from  the  system.", 
sat->name) ; 
delete  satellited; 

I 

else  sprintf (str, "%s  was  launched.",  sat->name) ; 
message (str) ; 
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/★  ********************************  *  j 
/*  */ 
/*  smet_main.c  */ 
/*  by  Carlos  I.  Noriega  */ 
/*  */ 
/*  This  file  main  routines  used  by  SMET.  */ 
/*  */ 
/*  Written  -  September  3,  1990  */ 
/*  Modified  -  September  3,  1990  */ 
/*  */ 
/**********************************/ 


# include  "smet.h" 


/*  MAIN  */ 

/*  Initialize  everything  and  enter  a  loop  to  check  the  */ 

/*  queue,  act  on  it  and  update  the  windows.  */ 

main (argc, argv) 
int  argc; 
char  **argv; 

< 

initialize_smet (argc, argv) ;  /*  windows,  menus,  queue,  variables  */ 

while  (TRUE)  /*  loop  until  we  exit  via  a  pop  up  menu  */ 

( 

handle_queue ( ) ; 

user_pre_update_satellites () ;  /*  User  defined  routine  */ 

update_satellites  () ; 

user_post_update_satellites  () ;  /*  User  defined  routine  */ 
draw_windows () ; 

if  (clock_running)  clock  -  clock  +  clock  speed;  /*  increment  clock  if  on  */ 

) 
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/**************  This  is  where  I  want  to  do  the  multi  processing  *********** / 


/*  UPDATE_SATELLITES  */ 
/*  Calling  routine  for  updating  all  the  satellites.  */ 
/*  Also  checks  for  earth  impact  and  sets  Automatic  */ 
/*  Viewing  position  if  auto_viewer  is  on.  */ 


update_satellites () 

( 

satellite  *sat,  *previous_sat ; 
float  max_apogee  —  1.2; 

previous_sat  **  £header_sat; 
launch_jpending_sats  ()  ; 

for (sat  »  header_sat . next ;  sat  ! “  NULL;  sat  —  sat->next) 

( 

if  (sat->maneuver  !-  NULL  £&  clock  >-  sat->maneuver->time)  maneuver_sat (sat) ; 
user_pre_orbital_parameters (sat,  clock) ; 
update_changing_parameters (sat,  clock) ; 
user_post_orbital_parameters (sat,  clock) ; 

if  (sat->range  <-  1.0)  crash_satellite (sat,  previous_sat) ; 
previous  sat  “  sat; 

max_apogee  =  max (sat->apogee_range,  max_apogee) ; 

) 

if  (auto_viewer )  /*  put  our  viewing  position  out  towards  Ares  */ 

( 

viewer_x  «  max_apogee; 
viewer_y  -  1.0; 
viewer_z  “  1.0; 

viewerrange  ”  sqrt (max_apogee*max_apogee  +  2) ; 

) 

) 


/*  DRAW_WINDOWS  */ 

/*  Calling  routine  for  all  the  drawing  routines.  */ 


draw_windows () 

l 


draw_omni_view  ()  ; 

/* 

draw_sat_view ( ) ; 

/* 

draw_merc_view ( ) ; 

/* 

draw_parameters () ; 

/* 

draw_sat_list ()  ; 

/* 

) 


draw  the  omniscient  point  of  view  */ 
draw  the  satellite  point  of  view  */ 
draw  the  mercator  map  with  plot  */ 
show  the  current  satellite's  stats  */ 
draw  the  list  of  satellites  */ 
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/* 

★  * 

*********************** 

********/ 

/* 

*/ 

/* 

smet  maneuver. c 

*/ 

/* 

by  Carlos  I.  Noriega 

*/ 

/* 

*/ 

/* 

This  file  contains  the  routines  used  by  SMET 

to  */ 

/* 

maneuver  a  satellite. 

*/ 

/* 

*/ 

/* 

Written  -  September  3,  1990 

*/ 

/* 

Modified  -  September  3,  1990 

*/ 

/* 

*/ 

/* 

★  * 

*********************** 

********/ 

♦include  "smet.h" 

/* 

ADD 

MANEUVER 

V 

/* 

Use 

a  series  of  dialogs  to  get  required  info. 

*/ 

/* 

See 

Chapter  4  for  more  information. 

*/ 

add 

maneuver  ( ) 

( 

char  str[25]; 

manstruct  "walker,  ‘maneuver; 

if  (current_sat  !«=  NULL) 

{ 

maneuver  **  (manstruct  *)malloc (sizeof (man_struct) ) ; 
if  (current_sat->maneuver  ■“  NULL) 

( 

dialog ("Enter  the  True  Anomoly  at  the  time  of  the  maneuver.",  str) ; 
maneuver->anom  —  atof (Str) *TO_RAD; 
current_sat->maneuver  -  maneuver; 
compute_maneuver_time (current_sat) ; 

) 

else 

( 

dialog ("Enter  the  number  of  degrees  to  travel  after  the  preceding  maneuver.",  str) 
maneuver->anom  *■  atof (str) *T0_RAD; 

for  (walker  -  current_sat->maneuver ;  walker->next  ! —  NULL;  walker  —  walker->next) 
walker->next  —  maneuver; 

) 

dialog ("Enter  (1)  for  Delta-V/Yaw/Pitch  or  (2)  Velocity/Flight  Path/Heading  Changes. 

,  str)  ; 

maneuver->normal_maneuver  -  atoi(str)  %  2; 
if  (maneuver->normal_maneuver) 

( 

dialog ("Enter  the  Delta-V  for  the  maneuver.",  str); 
maneuver->deltav  -  atof (str) /speed_f actor [unit_type] ; 
dialog ("Enter  the  Fitch  for  the  maneuver.",  str); 
maneuver->pitch  —  atof (str) *TO_RAD; 
dialog("Enter  the  Yaw  for  the  maneuver.",  str) ; 
maneuver->yaw  «  atof (str) *TO_RAD; 

) 

else 

< 

dialog ("Enter  the  velocity  Magnitude  Change  for  the  maneuver.",  str); 
maneuver->deltav  -  atof (str) /speed_factor [unit_type] ; 
dialog ("Enter  the  Flight  Path  Change  for  the  maneuver.",  str); 
maneuver->pitch  —  atof (str) *TO_RAD; 

dialog ("Enter  the  Heading  Change  for  the  maneuver.",  str); 
maneuver->yaw  -  atof (str) *T0_RAD; 

>  90 

maneuver->next  -  NULL; 

> 
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/*  MANEUVER_SAT  */ 

/*  See  Chapter  4  for  more  information.  */ 

maneuver_sat (sat) 
satellite  "sat; 

< 

float  delta_f_path,  delta_heading,  new_velocity, 
old_ecc  —  sat->ecc, 
oldjtnajor  -  sat->major, 
old_incl  “  sat->incl, 
old_peri  «■  sat->peri, 
old_a;ic  ■  sat->asc; 

update_changing_parameters (sat,  sat->maneuver->time) ; 

if  (sat->maneuver->normal_maneuver) 

( 

compute_deltas (sat,  Sdelta_f_path,  Sdelta_heading,  4new_velocity) ; 
f light_path_change (sat,  delta_f_path) ; 
heading_change (sat,  delta_heading) ; 
velocity_change (sat,  new_velocity) ; 

) 

else 

( 

f light_path_change (sat,  sat->maneuver->pitch) ; 
heading_change (sat,  sat->maneuver->yaw) ; 

velocity_change (sat,  sat->velocity  +  sat->maneuver->deltav) ; 

> 

if  (0.0  <-  sat->ecc  &4  sat->ecc  <  1.0) 

( 

sat->epoch  -  sat ->maneuver->t ime 

-  pom  (sat->ma  jo*.,  1 . 5)  *  (sat->E  -  sat->ecc  *  sin(sat->E))  *  TU; 
if  (sat  —  current_sat)  reset_dials () ; 
update_static_parameters (sat) ; 
update_sat_matri::  (sat)  ; 

) 

else 

( 

message ("That  maneuver  put  the  satellite  into  an  open  orbit  and  will  be  ignored."); 

sat->ecc  —  old_ecc; 

sat->major  -  old_major; 

sat->incl  —  old_incl; 

sat->peri  -  old_peri; 

sat->asc  —  old_asc; 

) 

delete  maneuver (sat )  ; 
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/*  COMPUTE_DELTAS  */ 

/*  See  Chapter  4  for  more  information.  */ 

compute_deltas (sat,  delta_f_path,  delta_heading,  new_yelocity) 
satellite  *sat; 

float  *delta_f_path,  *delta_heading,  *new_velocity; 

( 

float 

delvp  —  sat->maneuver->deltav  *  cos (sat->maneuver->yaw) , 
delvps  —  delvp  *  sin(sat->maneuver->pitch) , 
delvop  ■  sat->maneuver->deltav  *  sin (sat->maneuver->yaw) , 
vlmodl  ■  sat->velocity  +  delvp  *  cos (sat->maneuver->pitch) , 
vlmod2  ■■  f sqrt (delvps*delvps  +  vlmodl * vlmodl ) ; 

*delta_f_path  «  fatan2 (delvps, vlmodl) ; 

*delta_heading  ■»  fatan2 (delvop,  vlmod2) ; 

*new_velocity  -  fsqrt (delvop*delvop  +  vlmod2*vlmod2) ; 

) 

/*  COMPUTE_MANEUVER_TIME  */ 

/*  Set  the  time  of  maneuver  based  on  the  current  time,  */ 

/*  current  anomaly  and  desired  anomaly.  */ 

compute_maneuver_time (sat) 
satellite  *sat; 

< 

float  c_anom  —  cos (sat->maneuver->anom) , 

E  -  acos ( (sat->ecc  +  c_anom)/(l  +  sat->ecc  *  c_»nom) ) ;  /*  Eccentric  Anomoly  */ 

/*  number  of  revolutions  completed  plus  number  of  more  whole  revolutions  */ 
int  revs  ->  (int)  (clock  -  sat->epoch)  /sat->period 

+  (int) (sat->maneuver->anom/DEG360) ; 

/*  go  to  next  time  around  if  we're  past  it  on  this  revolution*/ 
if  (sat->maneuver->anom  <-  sat->anom)  revs++; 

/*  make  anomoly  positive  */ 

while  (sat->maneuver->anom  <  0)  sat->maneuver->anom  ■  sat->maneuver->anom  +  DEG360; 

/*  correct  quadrant  */ 

if  (fmod (sat->maneuver->anom, DEG360)  >  DEG180)  E  «■  DEG360  -  E; 

sat->maneuver->time  ■  sat->epoch  +  sat->period* (revs  +  (E  -  sat->ecc  *  sin (E) ) /DEG360) ; 

) 

/*  DELETE_MANEUVER  */ 

/*  Remove  the  maneuver  from  the  linked  list.  */ 

deletejmaneuver (sat) 
satellite  *sat; 

( 

man_struct  *temp; 

temp  -  sat->maneuver; 
sat->maneuver  _  sat->maneuver->next; 
free ((char  *)  temp) ; 

/*  set  execution  time  of  next  maneuver  if  there  is  one  */ 
if  (sat->maneuver  !■  HULL) 

( 

sat->maneuver->anom  *■  sat->maneuver->anom  +  sat->anom; 
compute_maneuver_time (sat)  ;  92 
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/*  FLIGHT_PATH_CHANGE  */ 

l*  See  Chapter  4  for  information.  */ 

flight_path_change (sat,  delta)  /*  In-Plane  Maneuver  */ 
satellite  *sat; 
float  delta; 
l 

float  p, 

old_anom, 
tan_f_path; 
char  str[70]; 

if  (nonzerop (delta) ) 

< 

tan_f_path  -  sat->ecc  *  sin <sat->E) /f sqrt (  1  -  aat->ecc  *  sat->ecc) ; 
sat->ecc  -  fsqrt(l  -  (1  -  sat->ecc*sat->ecc)  * 

po*«  ( (cos  (delta)  -  tan_f_path*sin  (delta) ) ,  2) )  ; 

old_anom  “  sat->anom; 

tan_f_path  -  (tan_f_path  +  tan (delta) )/ <  1  -  tan_f_path  *  tan (delta)); 
p  -  sat->ma jor  *  (1  -  sat->ecc  *  sat->ecc) ; 

sat->anom  —  fmod ( fatan2 (p  *  tan_f_path,p  -  sat->range)  +  DEG360,  DEG360) ; 

sat->E  —  acos((sat->ecc  +  cos (sat->anom) ) / (1  +  sat->ecc  *  cos (sit->anom) ) ) ; 

if  (sat->anom  >  DEG180)  sat->E  —  DEG360  -  sat->E;  /*  get  correct  quadrant  */ 

sat->peri  «=  fmod (sat->peri  +  old_anom  -  sat->anom  +  DEG360,  DEG360) ; 
sprintf (str, "%s  executed  a  flight  path  change  of  %4.3f  degrees.", 
sat->name,  delta  *  TO_DEG) ; 

message (str) ; 

) 

) 
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/*  HEADING_CHANGE  V 

/*  See  Chapter  4  for  information.  */ 

heading_change (sat,  delta)  /*  Out-of-Plane  Maneuver  */ 
satellite  *sat; 
float  delta; 

( 

float  incl_old  ”  sat->incl, 

u_old  -  sat->peri  +  sat->anom, 

unew  —  0.0, 

temp; 

char  strl70]; 

if  (nonzerop (delta) ) 

< 

sat->incl  -  facos (cos  <incl_old) *cos (delta)  +  sin (inel_old) *cos (u_old) *sin (delta) ) ; 
if  (nonzerop (sat->incl) ) 

< 

temp  ■■  sin  (incl_old)  *cos  (u__old)  *cos  (delta)  -  cos (incl_old) *sin (delta) ; 
u_new  —  facos (temp/sin (sat->incl) ) ; 

if  (sat->y  <  0.0)  u_new  «  DEG360  -  u_new;  /*  get  correct  quadrant  */ 

) 

sat->peri  *  fmod(u_new  -  sat->anom  +  DEG360,  DEG360) ; 

sat->asc  —  fmod (sat->asc  +  fatan2 (sin (u_old) *cos (incl_old) , cos (u_old) ) 

-  fatan2 (sin (u_new) *cos (sat->incl)  ,  cos (u_new) )  +  DEG360,  DEG360) ; 
sprintf (str, "%s  executed  a  heading  change  of  %4.3f  degrees .", sat->name, delta*TO_DEG)  ; 
me  3  sage (str) ; 

) 

> 


94 


it  maneuver .  c 


Hon  Sep  3  19:17:10  1990 


€ 


/*  VELOCITY_CHANGE  */ 

/*  See  Chapter  4  for  information.  */ 

velocity_change (sat,  new_vel) 
satellite  "sat; 
float  new_vel; 

( 

float  tan_f_path,  cos_f_path,  p,  old_anom; 
char  str [70]; 

if  (nonzerop  (ne**_vel  -  sat->velocity) ) 

( 

if  (sat->ecc  ”  0) 

if  (new_vel  >  sat->velocity) 

( 

sat->peri  —  fmod (sat->peri  +  sat->anom,  DEG360) ; 
sat->E  ”  sat->anom  “  0.0; 

) 

else 

{ 

sat->peri  =  fmod (sat->peri  +  sat->anom  +  DEG180,  DEG360) ; 
sat->E  —  sat->anom  *»  DEG180; 

) 

sat->ma jor  -  sat->range/ (2  -  sat->range  *  new_vel  *  new_vel) ; 

tan_f_path  -  sat->ecc  *  f sin (sat->E) /sqrt (1  -  sat->ecc  *  sat->ecc) ; 
cos_f_path  —  cos (atan (tan_f_rath) ) ; 

p  -  sat->range  *  sat->range  *  new_vel  *  new_vel  *  cos_f_path  *  cos_f_path; 
sat->ecc  «*  fsqrt  (1.0  -  p/sat->ma  jor)  ; 

oldanom  ”  sat->anom; 

sat->anom  —  fmod (fatan2 (p  *  tan_f_path,p  -  sat->range)  +  DEG360,  DEG360) ; 

sat->E  -  acos ( (sat->ecc  +  cos (sat->anflro) ) / (1  +  sat->ecc  *  cos (sat->anom) ) ) ; 

if  (sat->anom  >  DEG180)  sat->E  -  DEG360  -  sat->E;  /*  get  correct  quadrant  */ 

sat->peri  —  fmod (sat->peri  +  old_anom  -  sat->anom  +  DEG360,  DEG360) ; 

sprintf (str, "%s  executed  a  velocity  magnitude  change  of  %4.3f  %s . sat->name, 

(new_vel  -  sat->velocity) *speed_f actor [unit_typej , speed_unit [unit_type] ) 
message  (str) ; 

> 

) 
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/*  ********************************  */ 
/*  */ 

/*  amet_map.c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/*  Thia  file  containa  the  routinea  to  draw  the  globe  */ 

/*  uaed  by  the  SMET  demonatrator .  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/**********************************/ 

♦include  "amet.h" 

/*  DRAW-MERCATOR  */ 

/*  */ 


draw_mercator () 

{ 

c3f (Color [EARTH] ) ; 

rectf (-DEG180, -DEG90, DEG180, DEG90) ;  /*  map  background  */ 

draw_coaata2d () ; 
draw_bordera2d ()  ; 
draw_latlons2d ( )  ; 

) 


/*  DRAW-COASTS2D  */ 

/*  Draw  the  deaired  data  FINE  or  COURSE.  */ 

draw_coaata2d ( ) 

< 

int  c; 

linewidth (2) ; 
c3f (Color [COAST] ) ; 
if  (fineflag) 

for  (c  -  0;  c  <  FINEPOINTS;  C++] 
if  (fine2d_data .pen [c] ) 

draw(fine2d_data.x[c] ,  fine2d_data.y (c] ,  fine2d_data. z [c] ) ; 
else 

move [fine2d_data.x(c] ,  f ine2d_data .y [c] ,  f ine2d_data . z [c] ) ; 

else 

for  (c  -  0;  c  <  COARSEPOINTS;  C++) 
if  (coarae2d_data.pen[cJ ) 

draw (coarse2d_data .x [c] ,  coarse2d_data . y [c] ,  coarse2d_data . z [c] ) 

else 

move  <coarse2d_data.x[c] ,  coarse2d_data.y[cJ ,  coarse2d_data. z [c] ) 

> 
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/*  DRAW- BORDERS 2D  */ 

/*  Draw  the  border  data  if  desired.  */ 

draw_borders2d () 

( 

int  c; 

linewidth (1) ; 
c3f (Color (BORDER] ) ; 
if  (borderflag) 

for  (c  -  0;  c  <  BORDERPOINTS;  C++) 
if  (border2d_data.pen[c] ) 

draw (border 2d_data.x[c] ,  border 2d_data. y [c] ,  border2d_data. z (c) ) ; 
else 

move (border2d_data .x[c] ,  border2d_data . y [c] ,  border2d_data . * [c] ) ; 


/*  DRAW-LATLONS2D  */ 

/*  Draw  the  latitude  and  longitude  lines  on  the  map.  */ 

draw_latlons2d () 

( 

int  c; 

linewidth  (1) ; 

c3f  (Color (LATLONG] )  : 

for  (c  -  -6;  c  <-  6;  C++)  /*  draw  the  longitudes,  every  30  degrees  */ 

( 

move (DEG30  *  c,-DEG90,  0.0); 
draw (DEG30  *  c,  DEG90,  0.0); 

) 

for  (c  —  -3;  c  <-  3;  C++)  /*  draw  the  latitudes,  every  30  degrees  */ 

( 

move (-DEG180,  DEG30  *  c,  0.0); 
draw (  DEG180,  DEG30  *  c,  0.0); 

) 
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/* 

***************** 

******* 

*  *  *  *  * 

*  *  *  *  l 

/* 

*/ 

/* 

smet  orbit . c 

*/ 

/* 

by  Carlos  I.  Noriega 

*/ 

/* 

*/ 

/* 

This  file  contains  the  routine  used  by  SMET 

to 

*/ 

/* 

draw  the  orbit  ellipses  in  a  viewing  window. 

*/ 

/* 

*/ 

/* 

Written  -  September  3,  1990 

*/ 

/* 

Modified  -  September  3,  1990 

*/ 

/* 

*/ 

/* 

***************** 

******* 

*  *  *  *  * 

* 

♦ 

« 

« 

♦include  "smet.h" 

' 

/* 

DRAW-ORBIT 

*/ 

/* 

draw  an  orbital  ellipse  based  on 

the  global 

*/ 

/* 

parameters:  major,  minor  and 

eccentricity. 

*/ 

draworbit (sat) 
satellite  *sat; 

{ 

pushmatrix ( )  ; 

multmatrix (sat->matrix) ;  /*  set  the  orbital  transformation  plane  matrix  */ 

linewidth (3) ;  /*  set  the  linewidth  */ 

c4f (Color [sat->color] ) ;  /*  go  to  satellite's  color  */ 

if  (sat->vectors) 

(  /*  draw  the  normal  vector  */ 
move (0.0,  1.0,  0.0); 
draw  (0.0,  '5.0,  0.0); 

/*  draw  the  perigee  vector  */ 

move (0.0,  0.0,  1.0); 

draw(0.0,  0.0,  sat->perigee_range) ; 

) 

/*  shift  the  focus  to  the  origin  */ 

translate (0 . 0,  0,0,  -sat->ecc*sat->ma jor) ; 

/*  make  an  ellipse  out  of  a  circle  */ 
scale (sat->minor,  1.0,  sat->ma jor) ; 

/*  put  it  on  equatorial  plane  */ 

rotate (900, '  x' )  ; 

if  (sat->fill) 

(  /*  draw  a  simple  transparent  circle  * / 

zbuffer (FALSE) ; 

blendf unction (BF_SA,  BF_MSA); 

circf (0 . 0, 0 . 0, 1 . 0) ; 

blendfunction (BF_ONE,  BF_2ERO) ; 

zbuffer (TRUE) ; 

) 

else  circ (0 . 0,  0 . 0, 1 . 0)  ;  /*  draw  an  unfilled  oval  */ 

popmatrix ( ) ; 

) 
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Z  **********************************/ 


/*  */ 

/*  smet_param.c  */ 

/*  by  Carlos  X.  Noriega  */ 

/*  *  / 

/*  This  file  contains  the  routines  used  by  SMET  to  */ 

/*  handle  parameter  changes.  *f 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 


♦include  "smet.h" 


smet_param_option (n)  /*  used  by  popup  menu  */ 

i  nt  n  ; 

( 

param_option  =  n; 

) 


smet_sat_color (n)  /*  used  by  popup  menu  */ 

int  n; 

{ 

current_sat->color  -  SAT_RED  +  n; 

J 


/* 

CHANGE_PARAMETER 

*/ 

/* 

Compute  the  parameter 

number  base 

on  the  mouse 

*/ 

/* 

position  and  call  the 

appropriate 

change  routine. 

V 

change_parameter (y) 
int  y; 

( 

int  paramnum  =  (float) (satlistwin. sizey  -  y)  /  (float) LINE  -  4; 

if  (0  <—  paramnum  &&  paramnum  <  25) 
switch  (paramoption)  ( 

case  CLASSICAL:  change_classical_param (paramnum) ;  break; 

case  SECONDARY:  change_secondary_param (paramnum) ;  break; 

case  POSITIONS:  change_posit ion_par am (paramnum) ;  break; 

case  OPTIONS:  change_options_param (paramnum) ;  break; 

case  COMBINED:  change_combined_param (paramnum) ;  break; 

) 

else  change_sysparara (paramnum  -  29) ; 

> 
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/*  CHANGE_xxxxxxxx_PARAM  * / 

/*  Change  the  appropriate  parameter.  */ 

change_classical_param(paramnum) 
int  paramnum; 

{ 

switch (paramnum) 

{ 

case  0:  ask_name();  break; 

case  1:  ask_major();  update_static_parameters (current_sat) ;  break; 
case  2:  ask_minor_for_ecc ()  ;  update_static_pararaeters (current_sat) ;  break; 
case  3:  ask_ecc();  update_static_parameters (current_sat) ;  break; 
case  4:  ask_incl();  update_sat_raatrix(current_sat) ; 

updatejperturb_constants (current_sat) ;  break; 
case  5:  ask_asc<);  update_sat_matrix (current_sat) ;  break; 
case  6:  ask_jjeri();  update_sat_matrix (current_sat) ;  break; 

case  7:  bad_choice_message (classicalparamblock.params [paramnum] .name) ;  break; 
case  8:  ask_epoch();  delete_plots (current_sat) ;  break; 

) 

if  (1  <«  paramnum  &&  paramnum  <-  6) 

{ 

reset_dials () ; 
delete_plots (current_sat) ; 

) 

) 


change_secondary_param (paramnum) 
int  paramnum; 

{ 

if  (paramnum  0)  ask_name(); 

else  if  (paramnum  —  3)  ask_epoch(); 
else  if  (paramnum  <-  6) 

bad_choice_message (secondaryparJOTblock.params (paramnum) .name) ; 

) 


change_position__param  (paramnum) 
int  paramnum; 

< 

if  (paramnum  •—  0)  ask_name(); 
else  if  (paramnum  <■■  7) 

bad  choice_message (positionparamblock.params [paramnum] .name) ; 

) 


change_options_param (paramnum) 
int  paramnum; 

( 

char  str[50); 

switch (paramnum) 

{ 

case  0:  ask_name();  break; 
case  1:  dopup (color_menu) ;  break; 

case  2:  current_sat->fill  —  !current_sat->fill;  break; 

case  3:  current_sat->vectors  -  ! current_sat->vectors;  break; 

case  4:  current_sat->plot  ”  !current_sat->plot;  delete_plots (current_sat) ; 
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change_combined_param (paramnum) 
int  paramnum; 

( 

if  (paramnum  <—  24)  switch (paramnum) 

( 

case  0:  ask_name();  break; 

case  1:  ask_major();  update_static_parameters (current_sat) ;  break; 
case  2:  ask_minor_for_ecc () ;  update_static_parameters (current_sat) ;  break; 
case  3:  ask_ecc();  update_static_parameters (current_sat) ;  break; 
case  4:  ask_incl();  update_sat_matrix (current_sat) ; 

update_perturb_constants (current_sat) ;  break; 
case  5:  ask_asc();  update_sat_matrix (current_sat) ;  break; 
case  6:  ask_peri();  update_sat_matrix (current_sat) ;  break; 
case  8:  ask_epoch();  delete_plots (current_sat) ;  break; 
case  21:  dopup (color_menu) ;  break; 

case  22:  current_sat->fill  -  ! current_sat->fill;  break; 
case  23:  current_sat->vectors  -  ! current_sat->vectors;  break; 

case  24:  current_sat->plot  ”  ! current_sat->plot ;  delete_plots (current_sat) ;  break; 
default:  bad_choice_message (combinedparamblock . params [paramnum]  .name) ;  break; 

I 

if  (1  O  paramnum  &&  paramnum  <—  6) 

( 

reset_dials () ; 
delete_plots (current_sat) ; 

) 

) 
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change_sysparam (paramnum) 
int  paramnum; 

( 

char  str (15]; 
satellite  "sat; 

switch (paramnum) 

( 

case  0:  clock_running  —  ! clock_running;  break; 
case  1 : 

clock  -  ask  time ("current  time"); 

for  (sat  —  current_sat;  sat  !-  NULL;  sat  -  sat->next)  delete_plots  (sat) ; 
break; 
case  2: 

dialog ("Enter  a  clock  speed  (in  minutes  per  update) ", str) ; 
clock_speed  —  atof (str); 

set valuator (DIAL3,  (short) (clock_speed) ,  0,  max ( (short) (clock_speed) "1.2,  360)) 
break; 

case  3:  fineflag  —  !fineflag;  break; 
case  4:  borderflag  =  Iborderflag;  break; 
case  5:  perturbations  =  ! perturbations;  break; 
case  6:  auto_viewer  -  ! auto_viewer ;  break; 
case  7: 

dialog ( "Enter  new  viewer  X  position .", str) ; 
viewer_x  =  atof (str) /distance_f actor [unit_type] ; 

viewer_range  -  sqrt  (viewer_x*viewer__x  +  viewer_y"viewer_y  +  viewer_z*viewer_z) ; 
auto_viewer  —  TALSE; 
break; 
case  8: 

dialog ("Enter  new  viewer  Y  position. ", str) ; 
viewer_y  «  atof (str) /distance_f actor [unit_type] ; 

viewer_range  «  sqrt (viewer_x*viewer__x  +  viewer_y*viewer_y  +  viewer_z*viewer_z) ; 
auto_viewer  —  TALSE; 
break; 
case  9: 

dialog ("Enter  new  viewer  Z  position str ) ; 
viewer_z  —  atof (str) /distance_f actor [unit_type] ; 

viewer_range  —  sqrt (viewer_x*viewer__x  +  viewer_y*viewer_y  +  viewer_z*viewer_z) ; 

auto_viewer  —  FALSE; 

break; 

case  10:  unit_type  —  (unit_type  +  1)  %  3; 

) 

) 


*/ 
*/ 

bad_choice_message (name) 
char  "name; 

( 

char  str(50); 

sprintf (str, "%scannot  be  changed  by  the  user .", name) ; 
message (str) ; 


/*  BAD_CHOICE_MESSAGE 
/* 
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/*  AS K_xxxxxx  -  Custom  dialog  routines  used 
/*  to  ensure  the  correct  input  is  made. 

£loat  ask_time (name) 
char  ‘name; 

< 

char  msg[75], 

response [20] ; 
int  hour  —  0, 
min  —  0; 
float  sec  —  0.0, 
time; 

sprint f (msg,  "Enter  the  %s.",name); 
dialog (msg,  response); 
sscanf (response, "%d%*c%d%*c%f ", Shour,  &min,  (see); 
if  (min  >—  0  &&  min  <  60  &£  sec  >-  0.0  &&  sec  <  60.0) 

if  (hour  >—  0  &&  response [0]  !  —  '-')  time  -  hour"60  +  min  +  sec/60; 
else  time  -  hour*60  -  min  -  sec/60; 
else 
( 

sprintf(msg,  "The  %s  must  be  entered  in  00:00:00  format.",  name) ; 

message (msg) ; 

time  —  ask_time (name) ; 

) 

return  (time)  ; 


*/ 

*/ 


askname ( ) 

( 

char  str[45]; 

sprintf(str,  "Enter  new  name  for  satellite  %s  current_sat->name) ; 

dialog (str, current_sat->name) ; 


askjma jor () 

( 

char  str [15]; 

dialog ("Enter  a  Semi-Major  Axis.", str); 

current_sat->ma jor  -  atof (str) /distance_factor [unit_type] ; 
if  (current_sat->ma jor  <“  0.0) 

[ 

message ("Semi -Major  Axis  must  greater  than  0."); 
askjma jor  () ; 

) 


ask_minor_for_ecc  () 

( 

char  str [15]; 

dialog ("Enter  a  Semi -Minor  Axis.", str); 

current_sat->minor  -  atof (str) /distance_factor [unit_type] ; 

current_sat->ecc  -  sqrt (  1  -  pow (current_sat->minor  /  current_sat->ma jor,  2)) 
if  (current_sat->minor  <-0.0  II  current_sat->minor  >  current_sat->ma jor) 

< 

message ("Semi -Minor  Axis  must  be  between  0  and  the  Semi -Major  Axis."); 
ask  minor  for  ecc(); 

}  "  "  "  103 
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ask_ecc() 

< 

char  str [15]; 

dialog ("Enter  an  Eccentricity .", str) ; 
current_sat->ecc  *  atof(str); 

if  (current_sat->ecc  <0.0  ||  current_sat->ecc  >«  1.0) 

1 

message ("Eccentricity  must  be  between  0.0  and  1.0."); 
ask_ecc () ; 

) 


ask_incl () 

( 

char  str [15]; 

dialog ("Enter  an  Inclination. ", str) ; 
current_sat->incl  —  atof  (str)  *  TO_RAD; 

if  (current_sat->incl  <0.0  ||  current_sat->incl  >  DEG180) 

{ 

message ("Inclination  must  be  between  0  and  180  degrees."); 
ask_incl () ; 

) 

) 


askasc () 

( 

char  str [15]; 

dialog ("Enter  an  Ascending  Node.", str); 
current_sat->asc  -  atof (str)  *  TORAD; 

if  (current_sat->asc  <  0.0  ||  current_sat->asc  >  DEG360) 

< 

message ("Ascending  Node  must  be  between  0  and  360  degrees."); 
ask_asc () ; 

) 

) 


ask_peri () 

( 

char  str [15]; 

dialog ("Enter  an  Argument  of  Periapsis .",  str) ; 
current_sat->peri  -  atof (str)  *  T0_RAD; 

if  (current_sat->peri  <0.0  !|  current_sat->peri  >  DEG360) 

I 

message ("Argument  of  Periapsis  must  be  between  0  and  360  degrees.") 
ask_peri () ; 

) 

) 

ask_epoch () 

< 

current_sat->epoch  -  ask_time ("epoch") ; 

I 
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/**********************************/ 
/*  */ 

/ *  smet_queue . c  *  / 

/*  by  Carlos  X.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  routines  used  by  SMET  to  */ 

/*  check  the  queue  and  handle  events  not  covered  */ 

/*  in  other  files.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/*  ********************************  */ 

♦include  "smet.h" 

/*  GET_DIAL_VALUE  */ 

/*  Given  a  dial  and  an  initial  reading  from  the  queue  keep  */ 

/*  reading  it  until  the  user  stops  tuning  it.  */ 

float  get_dial_value (dev, initial_value) 

Device  dev;  /*  which  dial  we  are  working  with  */ 

short  initial_value ;  /*  value  when  first  detected  on  the  queue  */ 

( 

/*  set  the  current  value  of  the  dial  */ 

short  value  «  initial  value; 


/*  read  the  queue  until  something  other  than  the  */ 
/*  same  dial  is  on  top  or  the  queue  is  empty  */ 
while  (qtest  ()  ■*—  dev)  qread  {£ value )  ; 

if  (dev  !-  DIAL3)  delete_plots (current_sat) ; 

/*  return  the  last  value  in  float  format  */ 

return  '  < float) value) ; 


/*  RESETJDIALS  */ 

/*  Used  when  a  value  of  a  classical  parameter  is  changed  */ 

/*  without  the  use  of  the  dials.  Keeps  everyone  happy.  */ 

reset_dials () 

( 

/*  Set  the  dials  allowing  for  degree  of  accuracy  desired  */ 
set valuator (DIAL2,  (short) (current_sat->peri*TO_DEG*10) ,  0,  360*10); 
setvaluator (DXAL4,  (short) (current_sat->incl*TO_DEG*10) ,  0,  180*10); 
setvaluator (DIAL5,  (short) (current_sat->asc*TO_DEG*10) ,  0,  360*10); 
setvaluator (DIAL6,  (short) (current_sat->ma jor*100) ,  100,  1000); 
setvaluator (DIAL7,  (short) (current_sat->ecc*1000) ,  0,  999); 
qresetO;  /*  suppress  "setvaluator"  induced  events.  */ 
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/*  HANDLE_QUEUE  */ 

/*  Clears  the  queue  of  all  pending  events.  */ 


handle_queue () 
{ 


short 

value; 

/* 

value 

from  a  device  on  the  queue 

*/ 

while 

(qtest () ) 

/* 

clear 

queue  before  updating  windows 

*/ 

switch (qread ((value) ) 

( 

case  MIDDLEMOUSE: 

if  (value  —  1  )  handle_middlemouse () ; 
break; 

case  RIGHTMO0SE: 

if  (value  —  1  )  dopup (topjmenu) ; 
break; 

case  DIAL2:  /*  Adjust  "peri"  */ 

if  (currentsat  !■  NULL) 

( 

current_sat->peri  —  get_dial_value (DIAL2, value) *TO_RAD/10; 
update_sat_matrix (current_sat) ; 

) 

break; 

case  D1AL3 :  /*  Adjust  "clock  speed"  */ 

clock_speed  -  get_dial_value (DIAL3, value) ; 
break; 

case  D1AL4:  /*  Adjust  "incl"  */ 

if  (current_sat  !-  NULL) 

{ 

current_sat->incl  ■  get_dial_value (DIAL4, value) *TO_RAD/10; 
update_sat_matrix (current_sat) ; 
update_perturb_constants (current_sat) ; 

) 

break; 

case  DIAL5:  /*  Adjust  "asc"  */ 

if  (current_sat  !—  NULL) 

( 

current_sat->asc  ■  get_dial_value (DIALS, value) *TO_RAD/10; 
update_sat_matrix (current_sat)  ; 

> 

break; 

case  DIAL6:  /*  Adjust  "major"  */ 

if  (current_sat  !■  NULL) 

( 

current_sat->ma jor  -  get_dial_value (DIAL6, value) /100; 
update_static_parameters (current_sat) ; 

) 

break; 

case  DIAL7 :  /*  Adjust  "ecc"  */ 

if  (current_sat  !«  NULL) 

( 

current_sat->ecc  —  get_dial_value (DIAL7, value) /1000; 
update_static_parameters (current_sat) ; 

) 

break; 

default:  break; 

)  /*  end  of  switch  statement  */ 
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/*  PROGRAM-EXIT  */ 

/*  Exit  the  program  cleanly  */ 

program_exit ( ) 

( 

char  strflOJ; 

dialog("Would  you  like  to  save  any  satellites  before  quitting?", str) ; 

if  (str [0]  “  ' y'  ||  str[0]  —  'Y' )  ask_and_write_sat_file () ; 

message ("Thank  you  for  using  SMET."); 

winclose (satlistwin.gid)  ; 

winclose (paramwin . gid) ; 

winclose (auxwinl .gid) ; 

winclose (auxwin2 .gid) ; 

winclose (mainwin.gid) ; 

winclose (dialogwin . gid) ; 

winclose (backwin.gid)  ; 

gexit();  /*  Close  graphics  ports  */ 
exit();  /*  Exit  orbital_demo  */ 

) 


/*  HANDLE -MIDDLEMOUSE  */ 

I*  Call  the  appropriate  routine  based  on  the  */ 

/*  mouse  position.  */ 

handle_middlemouse ( ) 

( 

short  x,y;  /*  mouse  position  */ 

qread (&x) ;  /*  cetMOUSEX  */ 

qread(Sy);  /*  get  MOUSEY  */ 


/*  Are  we  in  the  area  containing  the  aux  windows?  */ 

if  (  x  >—  auxwinl. orgx  &&  x  <—  auxwin2 . orgx+auxwin2 . sizex) 

( 

if  (  y  >“  auxwinl. orgy  SS  y  <—  auxwinl . orgy+auxwinl . sizey  ) 
switch_windows (x) ; 

) 

else  if  (  x  <•=  paramwin. orgx  +  paramwin. sizex)  change_parameter (y) ; 
else  if  (  x  >—  satlistwin.orgx)  switch_current_sat (y) ; 

) 
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*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 

switch_windows  ( x ) 
short  x; 

( 

/*  Are  we  in  the  left  aux  window?  */ 

if  (  x  <■  auxwinl.orgx  +  auxwinl . sizex  ) 

switch  (window_order) ( 

case  OMNISATMERC: 

sat_view_gid  -  mainwin.gid; 
omni_view_gid  “  auxwinl. gid; 
window_order  ”  SATOMNIMERC; 
break ; 

case  OMNIMERCSAT: 

merc_view_gid  -  mainwin.gid; 
omni_view_gid  -  auxwinl. gid; 
window_order  -  MERCOMN1SAT; 
break; 

case  SATMERCOMNI: 

merc_view_gid  -  mainwin.gid; 
sat_view_gid  —  auxwinl. gid; 
window_order  -  MERCSATOMNI; 
break; 

case  SATOMNIMERC: 

omni_view_gid  -  mainwin.gid; 
sat_view_gid  •  auxwinl. gid; 
windoworder  ”  OMNISATMERC; 
break; 

case  MERCOMNISAT: 

omni_view_gid  ■  mainwin.gid; 
merc_view_gid  -  auxwinl. gid; 
window_order  —  OMNIMERCSAT; 
break ; 

case  MERCSATOMNI: 

sat_view_gid  -  mainwin.gid; 
merc_view_gid  »  auxwinl. gid; 
window_order  “  SATMERCOMNI; 
break; 


/*  SWITCH-WINDOWS 

/*  See  if  we  are  in  one  of  the  auxiliaty  windows.  If  we  are, 
/*  then  swap  the  gid  values  for  the  views  associated  with 

/*  the  current  main  window  and  the  auxiliary  window  which 

/*  was  se'scted  by  the  user  by  pressing  the  MIDDLEMOUSE. 

/*  The  window_order  variable  designates  the  which  views 

/*  are  in  which  windows.  For  example,  OMNISATMERC 

/*  means  that  the  main  window  contains  the  Omniscient 

/*  point  of  view,  auxiliary  window  1  contains  the 

/*  Satellite  point  of  view,  and  auxiliary  window  2 

/*  contains  the  meracator  projection. 
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/*  Are  we  in  the  right  aux  window? 
else  if  <x  >-  auxwin2.orgx  ) 

switch  (window_order) { 

case  OMNISATMERC: 

mere  view_gid  “  mainwin.gid; 
omni~view_gid  -  auxwin2.gid; 
window_order  —  MERCSATOMNI; 
break; 

case  OMNIMERCSAT: 

sat  view_gid  -  mainwin.gid; 
omniviewgid  -  auxwin2.gid; 
window_order  »  SATMERCOMNI; 
break ; 

case  SATMERCOMNI: 

omni  view _ gid  *■  mainwin.gid; 

sat  view_gid  -  auxwin2.gid; 
window  order  “  OMNIMERCSAT; 
break ; 

case  SATOMNIMERC: 

mere  view_gid  ~  mainwin.gid; 
sat  view_gid  -  auxwin2.gid; 
window  order  “  MERCOMNISAT; 
break; 

case  MERCOMNISAT: 

sat  view  gid  ”  mainwin.gid; 
mereviewgid  -  auxwin2.gid; 
window_order  ™  SATOMNIMERC ; 
break; 

case  MERCSATOMNI: 

omni  viewgid  =  mainwin.gid; 
merc~view_gid  =  auxwin2.gid; 
window_order  “  OMNISATMERC; 
break; 

> 

) 
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/*  ********************************  */ 


/*  */ 

/*  smet_sat.c  */ 

/*  by  Carlos  X.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  routines  used  by  SMET  to  */ 

/*  draw  the  satellite  in  a  viewing  window.  It  was  borrowed  */ 

/*  from  the  ASAT  program  written  by  Jim  Zanoli  and  Do  Kyeing  */ 

/*  Ok.  It  was  modified  and  converted  to  procedures.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 


/**********************************/ 
♦include  "smet.h" 


/*  DRAW-SATELLITE  */ 
/*  draws  a  satellite  in  the  correct  attitute  in  */ 
/*  the  appropriate  spot  in  the  orbital  plane  */ 
/*  as  set  by  the  global  orbital  parameters  */ 


/*  anom,  r_p  and  r_q.  The  satellite  is  scaled  */ 

/*  with  respect  to  the  major  axis  of  the  orbit.*/ 

draw_satellite (sat) 
satellite  *sat; 

(  /*  scale  size  increases  with  distance  */ 

float  sat_scale  —  0.005  *  viewer_range;  /*  maybe  vr*vr  ??  */ 

short  angle;  /*  dummy  loop  counter  for  making  the  ball  */ 

pushmatrix ( ) ; 

multmatrix (sat->matrix) ;  /*  set  the  orbital  transformation  plane  matrix  */ 

linewidth (2) ;  /*  set  the  linewidth  */ 


translate (sat->q, 0 .0, sat->p) ; 
rot (sat->anom*TO_DEG,  'y' ) ; 
rotate (900, ' x' ) ;  /* 

scale (sat  scale, sat  scale, sat 


/*  move  it  in  the  orbital  plane  */ 

/*  point  towards  the  earth 
lay  it  on  its  side  */ 

scale) ;  /*  make  much  smaller  than  the  earth 


c3f (Color (WHITE] ) ;  /* 

go 

white 

/*  draw  left  wing  */ 
pushmatrix () ; 
scale  (3 . 0,  0 . 7, 0 . 1)  ; 

/* 

make 

translate (-1.2, 0.0, 0.0) ; 

/* 

move 

draw_cube  () ; 
popmatrix  () ; 

/*  draw  right  wing  */ 
pushmatrix ( ) ; 
scale (3. 0,0. 7, 0.1) ; 

/* 

make 

translate (1 . 2, 0 . 0, 0 . 0)  ; 

/* 

move 

draw_cube ( ) ; 

popmatrix ()  ; 

c3f (Color [ sat->color J )  ; 

/* 

*/ 


it  long,  flat,  and  tt.  . n  */ 
it  out  to  the  side  */ 


it  long,  flat,  and  thin  */ 
it  out  to  the  side  */ 


make  the  body  the  right  color  */ 
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/*  draw  satellite  body  */ 
draw_cylinder ( ) ; 

/*  draw  satellite  middle  globe  part  */ 
pushmatrix () ; 

translate (0 . 0, -3 . 5, 0 . 0) ;  /*  move  below  the  main  cylinder  */ 

c3f (Color [YELLOW] ) ;  /*  go  yellow  */ 

/*  draw  a  0.8  radius  ball  */ 

for  ( angle«*0 ;  angle<-3 600 ;  angle^angle+200 ) 

( 

rotate (300, '  x'  )  ; 
circf (0.6, 0.0, 0.8); 
rotate (300,fy'); 
circf (0 . 0, 0 . 0, 0 . 8) ; 
rotate (300, ' z' ) ; 
circf (0.0, 0.0, 0.8) ; 

) 

popmatrix ( ) ; 

c3f (Color [sat->color] ) ; 

/*  draw  front  small  body 
pushmatrix () ; 
translate (0.0, -5.0, 0.0)  ; 
scale (1.0, 0.3, 1.0) ; 
draw_cylinder () ; 
popmatrix  () ; 

popmatrix  () ; 


/*  make  the  body  the  right  color 

*/ 

/*  move  below  the  ball  */ 

/*  make  it  shorter  */ 


V 
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/*  DRAW-CYLINDER  */ 

/*  draw  a  6  unit  long  by  1  unit  radius  cylinder  */ 

/*  centered  at  the  origin  and  sitting  upright  */ 

draw_cylinder () 

( 

short  angle; 

/*  draw  the  sides  of  the  cylinder  */ 

/*  by  drawing  72  filled  rectangles  */ 
pushmatrix () ; 

for (angle—0 ; angle<— 3  600 ; angle«angle+5  0 ) 

{ 

rotate (50, 'y' ) ; 
rectfi (-1,-3, 1,  3) ; 

) 

popmatrix ( )  ; 

/ *  draw  the  top  of  the  cylinder  */ 

pushmatrix ( )  ; 

translate (0 . 0, 3 . 0, 0 . 0) ;  /*  move  it  to  the  top  */ 

rotate (900, 'x' ) ;  /*  lay  it  flat  */ 

circf i  (0, 0,1); 
popmatrix () ; 

/*  draw  the  bottom  of  the  cylinder  */ 
pushmatrix () ; 

translate (0 .0, -3 . 0, 0 .0) ;  /*  move  it  to  the  bottom  */ 

rotate (-900, 'x' ) ;  /*  lay  it  flat  */ 

circf i  (0, 0,1) ; 
popmatrix  () ; 
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l*  DRAW-CUBE  */ 

/*  draw  a  cube,  used  for  reshaping  to  any  box  shape  */ 

draw_cube ( ) 

( 

/*  draw  the  front  */ 
pushmatrix () ; 
translate (0.0, 0.0, 1.0) ; 
rectfi (-1, -1, 1, 1) ; 
popmatrix () ; 

/*  draw  the  back  */ 
pushmatrix () ; 
translate (1 . 0, 0 . 0, 0 . 0) ; 
rotate (900, '  y' )  ; 
rectfi(-l, -1,1,1) ; 
popmatrix () ; 

/ *  draw  the  right  */ 
pushmatrix () ; 
translate (0.0, 0.0, -1.0) ; 
rotate (1800, 'y'); 
rectfi  (-1, -1, 1,1) ; 
popmatrix ( ) ; 

/*  draw  the  left  */ 
pushmatrix () ; 
translate  (-1.0, 0.0, 0.0) ; 
rotate (-900, ' y' ) ; 
rectfi  (-1,-1, 1,1) ; 
popmatrix () ; 

/*  draw  the  top  */ 
pushmatrix () ; 
translate (0 . 0, 1 . 0, 0 . 0) ; 
rotate (-900, ' x' ) ; 
rectfi  (-1,-1, 1,1) ; 
popmatrix () ; 

/*  draw  the  bottom  */ 
pushmatrix ( ) ; 
translate (0.0, -1.0, 0.0) ; 
rotate (900, ' x' ) ; 
rectfi  (-1,  -1, 1, 1) 
popmatrix () ; 
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/**********************************/ 
/*  */ 

/*  smet_user.c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  user  routines  called  by  SMET.  */ 

/*  They  currently  do  nothing.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/**********************************/ 


♦include  "smet.h" 

userinit () 

l 

/*  ex.  read  sat  files,  init  variables  etc.  */ 

) 


user_new_sat () 

( 

/*  ex.  modify  new  current_sat  */ 

) 


user_pre_update_satellites () 

( 

/*  ex.  check  for  exit,  compare  sat  positions,  maneuver  sats,  launch  sats,  etc.  */ 

) 

user_po3t_update_3atellites () 

{ 

/*  ex.  check  ^^.r  exit,  compare  sat  positions,  maneuver  sats,  launch  sats,  etc.  */ 

) 


user_pre_orbital_parameters (sat,  time) 
satellite  *sat; 
float  time; 

( 


/*  ex.  perturb  orbit  etc.  */ 

) 

user_post_orbital_parameters (sat,  time) 
satellite  *sat; 
float  time; 

{ 

/*  ex.  perturb  orbit  etc.  */ 

) 
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/**********************************/ 
/*  */ 

/*  smet_windows . c  */ 

/*  by  Carlos  I.  Noriega  */ 

/*  */ 

/*  This  file  contains  the  routines  used  by  SMET  to  */ 

/*  draw  all  the  windows.  */ 

/*  */ 

/*  Written  -  September  3,  1990  */ 

/*  Modified  -  September  3,  1990  */ 

/*  */ 

/•A********************************/ 


♦include  "smet . h" 


/*  DRAW-OMNI-VIEW  */ 

/*  Draw  the  globe,  orbital  ellipse,  satellite  and  */ 

/*  and  the  IJK  axis  at  the  center  of  the  window  */ 

/*  form  an  omniscient  point  of  view.  */ 

draw_omni_view () 

« 

satellite  *sat; 

/*  select  the  proper  window  and  and  clear  it  */ 

winset (omni_view_gid) ; 

c3f  (Color [BLACK] ) ; 

zclear () ; 

clear ()  ; 

/*  make  the  box  large  enough  to  hold  any  possible  rotation  */ 

/*  of  the  orbital  plane.  Keep  the  globe  in  the  center.  */ 

ortho (-vie wer_range,  viewer_range,  -viewer_range,  viewer_range,  0,  100); 

loadmatrix (ident_matrix) ; 

/ *  draw  earth  background  */ 

pushmatrix  ()  ; 

look=»t(0.0,  0.0,  viewer_range,  0.0,  0.0,  0.0,  0.0); 
c3f (Color [EARTH] )  ; 
circf (0.0, 0.0, 1.0)  ; 
popmatrix () ; 

/*  place  our  viewer  */ 

lookat (viewer_y,  viewer_z,  viewer_x,  0.0,  0.0,  0.0,  0.0); 

/*  draw  the  orbital  ellipses  */ 

for(sat  ■>  current_sat;  sat  !-  NULL;  sat  -  sat->next)  draw_orbit (sat) ; 

/ *  draw  the  IJK  axes  */ 

draw_axes ( ) ; 

/*  draw  the  earth  w/  rotation  */ 

draw_globe (ROTATIOM_RATE  *  clock); 

/*  draw  the  satellites  */ 

for (sat  -  current_sat;  sat  ! ”  NULL;  sat  ■  sat->next)  draw_satellite (sat) 
swapbuf fers () ;  /*  show  it  */ 
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/*  DRAW-SAT— VIEW  */ 
/*  Show  that  portion  of  the  globe  that  is  visible  */ 
/*  to  the  satellite.  This  is  best  represented  */ 


/*  by  a  "fish-eye"  view  given  by  using  perspective.  */ 

draw_sat_view  < ) 

( 

float  de l_z ; 

/*  select  the  proper  window  and  clear  it  */ 

winset (sat_view_gid) ; 

c3f (Color [BLACK] )  ; 

zclear () ; 

clear  () ; 

/*  if  we're  above  the  surface  draw  it  */ 
if  (current_sat  !-  NULL  S&  current  sat->range  >  1.001) 

{ 

/*  compute  the  fov  based  on  the  range  and  make  the  */ 

/*  box  deep  enough  to  handle  the  near  (viewable)  */ 

/*  half  of  the  globe.  */ 

del_z  —  l/current_sat->range; 

perspective (asin (del_z) "3790/PI,  1.0,  (current_sat->range  -1) /2, current_sat->range) 
loadmatrix (ident_matrix) ; 

pushmatrix ( )  ; 

lookat(0.0,  0.0,  current_sat->range,  0.0,  0.0,  0.0,  0.0); 
translate (0 . 0,  0.0,  del_z) ; 

/*  draw  a  ring  around  the  globe  the  color  of  the  current  satellite  */ 

c3f (Color [current_sat->color] )  ; 

circf(0.0,  0.0,  1.02  *  sqrt(l  -  del_z*del_z) ) ; 

/*  draw  earth  background  */ 
c3f (Color [EARTH] ) ; 

circf(0.0,  0.0,  sqrt(l  -  del_z*del_z) ) ; 
popmatrix ( ) ; 


/*  look  from  the  satellite  to  the  center  of  the  earth  */ 

/*  remember  that  IJK  system  not  same  as  IRIS'  XYZ  system  */ 

lookat (current_sat->y,  current_sat->z,  current_sat->x,  0.0,  0.0,  0.0,  0.0); 

/*  draw  the  earth  w/  provided  rotation  */ 
draw_globe ( ROT AT I ON_RAT E  *  clock) ; 

> 

swapbuf fers () ;  /*  show  it  */ 


116 


mt  windows .  c  Mon  Sap  3  18:31:43  1990  3 

/*  DRAW-MERC-VIEW  V 
/*  Draw  a  mercator  projection  map  and  plot  all  */ 
/*  lat/long  */ 


draw_merc_view () 

( 

satellite  *sat; 

plot_struct  ^plotter; 

int  Jt;  /*  dummy  loop  counter  */ 

/*  get  the  right  window  */ 

winset (merc_view_gid) ; 

/ *  clear  out  any  viewing  or  projection  matrices  left  */ 
/*  on  the  window  and  set  up  a  simple  2D  screen  */ 

ortho2 (-DEG180-0 . 001, DEG180+0 . 001,  -DEG120,  DEG120)  ; 
loadmatrix (ident_matrix) ; 

c3f (Color [BLACK] ) ;  /*  go  to  black  and  clear  window  */ 

clear () ; 
zclear  ()  ; 


draw_mercator  () ; 

/* 

draw 

the  map 

*/ 

draw_sun  () ; 

/* 

draw 

golden  circle  for  sun 

*/ 

for (sat  =  current_sat;  sat  !=  NULL;  sat  “  sat->next)  if  (sat->plot) 

( 

c3f (Color [sat->color] ) ; 

for (plotter  **  sat->plots;  plotter  !-  NULL;  plotter  -  plotter->next) 
drawdot (plotter->lon,  plotter->lat) ; 

) 

swapbuf fers () ;  /*  show  it  */ 
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/*  DRAW-DOT  */ 

/*  Used  in  drawing  the  Mercator  map.  */ 

/*  Draw  a  dot  at  the  given  x,y  coordinate.  Since  it  */ 

/*  is  a  small  dot,  it  is  faster  to  draw  an  8  point  */ 

/*  circle  than  the  standard  circf  circle.  */ 

draw_dot (x,y) 

float  x, y;  /*  coordinates  of  center  */ 

< 

move (x,  y,  1.0); 
rpmv(-0. 00873, -0.02107,0.0)  ; 
rpdr<0. 01234, 0.01234, 0.0)  ; 
rpdr (0.0, 0.01745, 0.0) ; 
rpdr (-0.01234, 0.01234, 0.0)  ; 
rpdr (-0.01745,0.0,0.0)  ; 
rpdr(-0. 01234, -0.01234, 0.0)  ; 
rpdr (0.0,-0.01745,0.0)  ; 
rpdr (0.01234,-0.01234,0.0)  ; 
pclos () ; 

) 


/*  DRAW-SUN  */ 
/*  Used  in  drawing  the  Mercator  map.  */ 
/*  Draw  a  transparent  ellipse  to  represent  the  sun's  */ 
/*  position.  */ 


draw_sun () 

( 

float  sun_lat  -  SUN_INCL  *  sin (clock/SUN_TU) , 

sun_lon  -  fmod (clock  *  (1/SUN_TU  -  ROT AT I ON_RATE )  -  DEG180,  DEG360) ; 

c4f (Color (SUN]) ; 

zbuffer (FALSE)  ; 

blendf unction (BF_SA,  BF_MSA) ; 

circf (sun_lon,  sun_lat, 1 . 0) ?  /*  western  hemisphere  */ 

circf (sun_lon  +  DEG360,  sun_lat,  1.0);  /*  eastern  hemisphere  */ 
blendf unction (BF_ONE,  BF_ZERO) ; 
zbuffer (TRUE); 
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/*  DRAW-PARAMETERS  */ 

/*  Go  through  the  appropriate  parameter  block  and  print  */ 

/*  each  parameter  in  the  right  spot  in  the  parameter  window.  */ 

draw_parameters () 

{ 

ParBlock  "satjblock; 
int  k; 

winset (paramwin. gid)  ; 
c3f (Color (BLACK] ) ; 
clear  () ; 

ortho2 (0, paramwin. sizex,0, paramwin. sizey) ; 
loadmatrix (ident  matrix); 


/*  current  satellite  parameter  block  */ 
/*  dummy  loop  counter  */ 

/*  go  to  the  parameter  window  */ 


/*  draw  satellite  parameters  */ 

c3f (Color (WHITE] ) ;  /*  write  a  title  */ 

cmov2i(10,  paramwin. sizey  -  3  *  LINE); 
fmprstr ("Satellite  Parameters:"); 


if  (current_sat  !=  NULL) 
( 


switch  (param_option)  ( 


case  CLASSICAL: 
case  SECONDARY: 
case  POSITIONS: 
case  OPTIONS: 
case  COMBINED: 


sat_block 
sat_block 
sat_block 
sat_block 
sat  block 


Sclassicalparamblock;  break; 
Ssecondaryparamblock;  break; 
Spositionparamblock;  break; 
Soptionsparamblock;  break; 
Scombinedparamblock;  break; 


) 

for  (k  -  5;  k  <  sat_block->numparams  +  5;  k++) 

( 


) 


cmov2i(10,  paramwin . sizey  -  k  *  LINE); 

draw_parameter (sat_block->params [k  -  5] , current_sat->color, 

(int) current  sat- (int) Sheader 


) 


sat) 


/*  draw  system  parameters  */ 

c3f (Color [WHITE] ) ;  /*  write  a  title  V 

cmov2i(10r  paramwin. sizey  -  32  *LINE) ; 
fmprstr ("System  Parameters:"); 

for  (k  «  34;  k  <  sysparamblock . numparams  +  34;  k++) 

( 

cmov2i(10,  paramwin. sizey  -  k  *  LINE); 
draw_parameter (sysparamblock .params (k  -  34] , WHITE, 0) ; 

) 

swapbuffers <) ;  /*  show  it  */ 
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/*  DRAW-PARAMETER  */ 

/*  Draw  the  contents  of  a  parameter  record  in  the  */ 

/*  desired  color.  Offset  is  used  to  get  the  value  from  */ 

/*  the  desired  satellite.  */ 

draw_partmeter  (par am,  color__code, offset) 

Par  param;  /*  a  parameter  record  */ 

int  color_code,  /*  what  color  to  draw  the  data  */ 

offset;  /*  offset  of  the  string  pointer  */ 

[ 

c3f (Color [WHITE] ) ; 
fmpr st r  (param. name) ; 

param. value  -  (char  *)  ( (int) param. value  +  offset); 
c3f (Color [color_code] ) ; 

switch  (param. type)  {  /*  choose  the  correct  output  routine  */ 

case  DISTANCE:  draw_distance (param. value) ;  break; 

case  SPEED:  draw_speed (param. value) ;  break; 

case  FLOAT :  dr aw_float (param. value) ;  break; 

case  DEGREE:  draw_degree (param. value) ;  break; 

case  LAT:  drawlat (param. value) ;  break; 

case  LON:  draw_lon (param. value) ;  break; 

case  TIME:  dr aw_time (param. value) ;  break; 

case  STRING:  fmprstr (param. value) ;  break; 

case  COLORCODE:  fmprstr (color_names{* (int  *) param. value] ) ;  break; 
case  UNITCODE:  fmprstr ( unit_names [ * (int  *) param. value] ) ;  break; 

case  FLAG:  if(*(int  *) param. value)  fmprstr ("On") ;  else  fmprstr ("Off") ;  break; 

) 


/*  DRAW-xxxxsrx  */ 

/*  Customized  formatting  routines  that  coerce  the  value  */ 

/*  into  the  format  we  want.  */ 

draw_distance (val) 
char  *val; 

( 

char  buf[15]; 


sprintf(buf,  "%4.2f  %s", 
fmprstr (buf)  ; 


* (float  *) val*distance_factor [unit_type] , 
distance_unit (unit_typej ) ; 


draw_speed (val) 
char  *val; 

( 

char  buf l 15]; 


sprintf(buf,  "%4.2f  %s", 
fmprstr (buf)  ; 


* (float  *) val*speed_factor [unit_type] , 
speed_unit (unit_type] ) ; 
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draw_f loat (val) 
char  *val; 

( 

char  bu£[10]; 

sprintf (buf ,  "%3.2£",  *  (float  *)val); 
fmprstr (buf) ; 


draw_degree (val) 
char  *val; 

( 

char  buf ( 10 ] ; 

sprint £ (buf,  "%3.2f",  *(£loat  *) val*TO_DEG) ; 
fmprstr (buf) ; 


draw_lat (val) 
char  *val; 

( 

char  buf [10]; 

if  (* (float  *)val  >—  0.0) 

sprintf (but,  "%3.2f  N",  fabs(*(float  *) val) *TO_DEG) ; 
else 

sprintf (buf,  "%3.2f  S",  fabs(* (float  *) val) *T0_DEG) ; 
fmprstr (buf) ; 


draw_lon (val) 
char  *val; 

( 

char  buf [10]; 

if  (* (float  *) val  >-  0.0) 

sprintf (buf,  "%3.2f  E",  fabs(* (float  *) val) *T0_DEG) ; 
else 

sprintf (buf ,  "%3.2f  W",  fabs(* (float  *) val) *TO_DEG) ; 
fmprstr (buf) ; 


draw_time  (val) 
char  *val; 

( 

char  buf [10]; 

int  hour  «■  * (float  *)val  /  60, 

min  -  * (float  *)val  -  hour  *  60, 

sec  -  (* (float  *)val  -  (hour  *  60  +  min))  *  60; 

if  (*(£loat  *)val  >-  0  I |  *(float  *)val  <-  -60) 

sprintf (buf,  "%02d: %02d: %02d",  hour,  abs (min) ,  abs(sec)); 

else  sprintf (buf,  "-0 : %02d: %02d",  abs (min) ,  abs (sec)); 

fmprstr (buf ) ; 


) 
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/*  DRAW-SAT-LIST  */ 

/*  List  out  the  linked  list  of  satellites  */ 

/*  vertically  in  the  satlistwin.  */ 

draw_sat_list  () 

( 

satellite  *sat; 

int  y  ”  satlistwin. sizey  -  3  *  LINE;  /*  first  line  of  text  */ 
winset (satlistwin.gid) ; 

c3f (Color [BLACK] ) ;  /*  clear  the  window  */ 

clear (> ; 

ortho2 (0, satlistwin. sizex, 0, satlistwin. sizey) ; 
loadmatrix (ident_matrix) ; 

cmov2i (20,  y)  ; 


c3f (Color (WHITE])  ; 

/* 

write  a  title 

*/ 

fmprstr ("SMET  Satellite  List: 
y  -  y  -  LINE; 

")  ; 

/* 

blank  line 

*/ 

for (sat  »  current  sat;  sat  !- 
( 

y  -  y  -  LINE; 
cmov2i  (25, y)  ; 
c3f (Color [ sat->color ] )  ; 

NULL;  sat  =  sat->next) 

/* 

current  sat  color 

*/ 

fmprstr (sat->narae)  ; 

) 

swapbuf fers ()  ; 

/* 

write  a  name 

*/ 

/* 

show  it 

*/ 

it  windows. c 


Mon  Sep  3  18:31:43  1990 


9 


/*  DRAW_INSTRUCTIONS  */ 

/*  Loop  through  the  instructions  array  and  draw  it  on  the  */ 

I*  Main  window,  then  wait  for  the  user  to  press  and  */ 

/*  release  the  rightmouse  button.  */ 

draw_inst ructions () 

( 

int  k;  /*  dummy  loop  variable  */ 

winset (mainwin . gid) ;  /*  get  the  main  window  */ 

c3f (Color [BLACK] ) ; 

clear  ()  ; 

zclear  () ; 

c3f (Color [CYAN]) ; 

/*  clear  out  any  viewing  or  projection  matrices  left  */ 

/*  on  the  main  window  and  set  up  a  simple  2D  screen  * / 
ortho2 (0, mainwin. sizex, 0, mainwin. sizey) ; 
loadmatrix (ident_matrix) ; 

for  (k  *  0;  k  <  28;  k++) 

( 

cmov2(75,700  -  k  *  LINE);  /*  move  to  the  left  margin  */ 
fmprstr (instructions [kj ) ;  /*  draw  the  line  of  text  */ 

) 

swapbuf fers 0 ;  /*  show  it  */ 

while  (qread(sk)  !—  RIGHTMOUSE);  /*  wait  for  press  */ 
while  (qread(sk)  ! »  RIGHTMOUSE);  /*  wait  for  release  */ 
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#/*  *********************************/ 
#/*  */ 

#/*  makefile  for  SMET  */ 

#/*  by  Carlos  1.  Noriega  */ 

#/*  */ 

f/*  Written  -  September  3,  1990  */ 

#/*  Modified  -  September  3,  1990  */ 

#/*  */ 


#/*  a********************************/ 

#  compiler  switches 

CC  -  cc 

INCLUDE  -  -I. 

CFLAGS  -  -O  -C  -D$ (LANGUAGE) 

LIBS  -  -lmalloc  -lfm  -lgl_s  -lm  \ 

image support / libnpsimage . a  ~zyda/utah/lib/librle . a  -limage 

SRCS  *  smet  addelete . c\ 

smet_axea . c\ 
smet_dialog .  c\ 
smet_f iles . c\ 
smet_globe . c\ 
smet_init .  . 

smetkepler . c\ 
smet_lciu:  _h.c\ 
amet  main.c\ 
amet_maneuver . c\ 

3met_map . c\ 
smet_orbit . c\ 
smet_param. c\ 
amet_queue . c\ 
ametaat . c\ 
ametuser . c\ 
ametwindows . c\ 
amet_vars . c 

OB JS  —  smet_addelete . o\ 

smet_axes .o\ 
amet_dialog . o\ 
smet_files .o\ 
smet_globe . o\ 
smet_init .o\ 
amet_kepler . o\ 
smet_launch . o\ 
smet_main.o\ 
smet_maneuver . o\ 
amet_map . o\ 
smet_orbit . o\ 
smet_param .  o\ 
smet_queue . o\ 
smet_sat .o\ 
smet_user.o\ 
smet_windows . o\ 
smet_vars .o 

smet:  $(OBJS) 

cc  -O  $ (INCLUDE)  $ (OBJS)  -o  smet  $ (LIBS) 

clean: 

rm  -£  $ (ALL) 

Secho  "make  clean  new  completed" 

neat : 

rm  -f  . [BC] *  .em*  *.o  *.BA*  *.CK*  core 

8echo  "make  neat  new  completed"  124 


i 


t 


.  c .  o : 


$<CC)  $(CrLAGS)  $ (INCLUDE)  $< 
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